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Predit par Baade et Zwicky en 1934, I'existence d'etoiles a neutrons a ete confirmee par 
la decouverte d'un pulsar radio en 1967. Cette decouverte fut un tournant important en 
astrophysique, tant par la confirmation d'objets macroscopiques de densite nucleaire dans 
I'univers, que par I'enjeu scientifique de comprendre la complexity d'une etoile a neutrons 
fortement magnetisee en rotation rapide, et par dela, d'etudier la physique sous des condi- 
tions extremes. Les pulsars sont formes au sein des supernovae ou par le transfert de moment 
angulaire d'une etoile compagnon vers une etoile a neutrons. Par ces deux processus de for- 
mation, on distingue deux types principaux de pulsars : les pulsars normaux (sujet principal 
de cette these) et les pulsars "milliseconde". 

Ces objets sont utilises pour etudier la physique des etoiles a neutrons, la relativite ge- 
nerale, le champ magnetique galactique, la physique des plasmas, la distribution de densite 
d'electrons dans la Galaxie, ainsi que la magnetohydrodynamique relativiste des vents de 
pulsars interagissant avec le milieu ambiant. Cependant, meme apres plus de 40 annees de 
recherches intenses, utilisant des instruments de plus en plus performants, il reste des ques- 
tions fondamentales qui ne sont pas resolues : 

- Combien de pulsars se trouvent dans la Galaxie et quel est leur taux de formation ? 

- Quelle est la relation entre la supernova d'une etoile massive et les proprietes intrin- 
seques des pulsars ? 

- Quelle est la limite superieure de la masse d'une etoile a neutrons ? 

- Est-ce que le vent de particules chargees injecte par le pulsar dans le mileu interstellaire 
peut alimenter une partie du spectre des rayons cosmiques ? 

- Quelle est la composition de la magnetosphere d'un pulsar et comment interagit-elle 
avec le champ magnetique intense ? 

- Ou et comment remission radio et remission a haute energie (de I'infrarouge au rayon- 
nement y) sont-elles produites dans la magnetosphere du pulsar ? 

- Est-ce que les processus physiques sont les memes pour tous les pulsars ? 

- Quel est le rapport entre les pulsars radio emetteurs y tels que le Crabe ou Vela, et les 
pulsars y non-emetteurs en radio tels que Geminga et CTAl ? 

L'etude de ces questions va bien au-dela du sujet de cette these, mais montrent I'etendue de 
la recherche sur les etoiles a neutrons et les pulsars. 

La periode de rotation des pulsars est relativement stable mais elle n'est pas constante. 
Tous ralentissent a cause de leur rayonnement magnetique dipolaire (ondes electromagne- 
tiques a la frequence de rotation du pulsar), et de leur vent de particules relativistes qui peut 
eventuellement creer en contact avec le milieu interstellaire une nebuleuse a vent de pulsar. 
Une fraction de ces particules accelerees le long des lignes de champ magnetique produit un 
faisceau d'ondes electromagnetiques qui balaye le ciel comme un phare, generant pour I'obser- 
vateur une pulsation par periode de rotation. La base de donnees "ATNF Pulsar Catalogue" 
repertorie 1827 pulsars repartis dans notre Galaxie et dans les nuages de Magellan, dont la 
pulsation a ete mise en evidence, principalement en radio (une cinquantaine emet egalement a 
haute energie, de I'optique au rayonnement y). Theoriquement, les ondes electromagnetiques 
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de la radio jusqu'axix hautes energies sont produites dans differentes regions de la magne- 
tosphere du pulsar. Les ondes radio polarisees sont emises depuis les poles magnetiques de 
I'etoile par emission synchrotron. Les photons y quant a eux sont produits principalement 
par le rayonnement de courbure et la diffusion inverse Compton des electrons (et positons) 
acceleres le long des lignes de champ magnetique. Cette propriete est importante puisque les 
particules chargees accelerees rayonnent leur energie tres vite et sont liees aux sites d'acce- 
leration. C'est dans ce cadre que se place ce travail de these, qui ambitionne d'analyser les 
pulsars y detectes par le telescope spatial Fermi Gamma-ray Space Telescope et d'essayer 
de contraindre la geometrie et les mecanismes des zones d'acceleration, en etudiant en de- 
tail les courbes de lumiere et les spectres resolus en phase d'une large population de pulsars y. 

La premiere partie de ce memoire presente les proprietes intrinseques des pulsars, ainsi 
que les modeles d'emission y bases sur les observations de I'experience EGRET. On dressera 
ensuite pour la mission Fermi une liste de pulsars candidats y etablie sur leur energie de 
ralentissement et de leur distance. 

Le Large Area Telescope, instrument principal de Fermi dedie a I'observation du ciel des 
hautes energies (> 20 MeV), sera brievement decrit dans la deuxieme partie de cette these. La 
mesure des parametres spectraux des pulsars (i.e. le flux integre, I'indice spectral, et I'energie 
de coupure), necessaires pour discriminer les modeles d'emission y des pulsars, depend des 
fonctions de reponse de 1 'instrument. Nous presenterons une methode developpee pour la 
validation de la surface efiicace en utilisant le pulsar de Vela, qui compare I'efScacite des 
coupures entre les vraies donnees et les donnees simulees a chaque niveau de la classification 
des evenements du LAT. Nous finirons cette partie par la description des analyses temporelle 
et spectrale des pulsars. 

Enfin, la troisieme et derniere partie presentera les resultats obtenus a partir des methodes 
d'analyses pour trois pulsars individuels, PSRs J0205+6449, J2229+6114, et J1048-5832, 
ainsi que du premier catalogue de pulsars du telescope Fermi. Une interpretation de la 
mesure des parametres spectraux et du profil des pulsations y des pulsars detectes, sera pre- 
sentee. 

Cette these se base sur les papiers suivants : 

- "The First Fermi Large Area Telescope Catalog of Gamma-ray Pulsars", 
Abdo, A. A. et al. 2009, submitted to Astrophysical Journal, 

Contact authors : A. Caliandro, E. Ferrara, D. Parent, R. Romani 

- "Discovery of Pulsations from the Pulsar J0205+6449 with the Fermi Gamma-ray Space 
Telescope", Abdo, A. A. et al. 2009, Astrophysical Journal Letters, 699, L102, 
Contact author : Damien Pcirent 

- "Fermi LAT detection of pulsed gamma-rays from the Vela-like pulsars PSR J 1048-5832 
and PSR J2229-I-6114", Abdo, A. A. et al. 2009, Astrophysical Journal, 706, 1331, 
Contact authors : Damien Pcirent, Alice K. Harding, Massimiliano Razzano 
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- "The Radio Polarization of Six Gamma-ray Pulsars seen with the Fermi LAT", 
Weltevrede, P. et al. 2009, submitted to Astrophysical Journal 

- "PSR J1410-6132 : a young, energetic pulsar associated with the EGRET source 3EG 
J1410-614r, O'Brien ... Parent, D., ... et al. 2008, Monthly Notices of the Royal As- 
tronomical Society, 388, LI 

- Signataire de 6 autres publications Science et ApJ sur les pulsars radio emetteurs y, 
6 autres sont en cours de publication. 
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Les Pulsars 
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1.1 Decouverte 

Depuis la fin de la seconde guerre mondiale, rastronomie radio connait un developpe- 
ment considerable. Les nouvelles techniques d'observations ouvrent de nouveaux champs de 



recherches et menent a la decouverte inattendue des quasars (Schmidt 1963). Suite a cela, 
Antony Hewish et son etudiante Jocelyn Bell de I'universite de Cambridge esperent differen- 
cier les quasars des radio-galaxies en utilisant la scintillation interplanetaire, refraction des 
ondes radio due a I'atmosphere terrestre et au gaz ionise interstellaire. lis construisent pour 
cela un radiotelescope suffisamment sensible aux sources faibles et variables dans le temps. 
Au cours du mois de juillet 1967, J. Bell decouvre une source qui emet de fagon constante 
et reguliere des impulsions radio qui different de la scintillation. EUes apparaissent avec un 



21 



Chapitre 1. Les Pulsars 



decalage de 4 minutes chaque jour, caracteristique d'un objet d'origine celeste. Un enregis- 
treur avec une reponse en temps plus rapide est ensuite utilise. La nouvelle prise de donnees 



confirme un signal extremement stable dont la periodicite est de 1.337 secondes (Figure 1.1). 



L'annonce de cette decouverte est publiee dans la revue Nature en fevrier 68 (Hewish et al 



1968), et soumet I'idee que la rotation ou I'oscillation d'un astre, tel que les naines blanches 



ou les etoiles a neutrons, pent etre a I'origine de ce rayonnement electromagnetique. Notons 
qu'aujourd'hui la source radio d'Hewish et Bell est connue sous le nom de PSR B1919+21 
ou PSR J1921+2153E] 




Hat] t'^H^ij^-rt'ltt 



r \ 



Figure 1.1 - Enregistrement des pulsations du premier pulsar connu, CP 1919 
(PSR B1919+21) , le 28 novembre 1967. Les signaux observes ont une periode de 1.337s. 



C'est seulement a la fin de I'annee 1968, apres la decouverte de signaux periodiques de 88 



ms et de 33 ms, respectivement dans les restes de supernova de Vela (Large et al 1968) et 



du Crabe ( iStaelin and Reifenstein ill 1968] ), que Ton attribue avec certitude I'origine de ce 
phenomene, non pas aux naines blanches relativement connues, mais a une etoile a neutrons en 
rotation. La periode minimale de rotation d'une etoile avec une densite uniforme p est donnee 
par r-min = (37t/Gp)^/^. Pour les naines blanches dont la densite moyenne est de 10 gem ^, la 



periode minimale de rotation est de I'ordre de la seconde. Quant a leurs oscillations, Melzer 



and Thorne (1966) montrent egalement des periodes minimales voisines de la seconde. En 



consequence, les signaux ne pouvaient venir que d'objets encore plus compacts, notamment 
les hypothetiques etoiles a neutrons. Le nom de pulsar est ainsi associe a cet objet, association 



d'autant plus pertinente qu'un ralentissement dans la periode du Crabe est observe (Richards 
and Comella, 1969), ce que prevoyait les mo deles de pulsar de Gold (1968) et Pacini ( 1968[ ). 

Cette remarquable decouverte confirmera les predictions faites en 1934 par [Baade and 
Zwicky (1934) de I'existence d'etoiles a neutrons au sein des supernovae (pour une revue 



Green ( [2009) ), stade ultime de revolution stellaire, et recompensera Hewish du 



recente, voir 

prix Nobel de physique en 1974. 

1. PSR signifie "Pulsating Source of Radio". B et J correspondent respectivement aux equinoxes standards 
utilisees pour les coordonnees du pulsar (B1950 et J2000), tandis que 1919 (ou 1921) et +21 (ou +2153) 
correspondent aux coordonnees de la source (ascension droite de 191il9' et declinaison de +21°). 
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1.2 Evolution de la Rotation et Proprietes Intrinseques 

Le lien entre etoiles a neutrons et pulsars etant etabli par I'observation de leur pulsation, 
interessons-nous a present a la formation et aux proprietes intrinseques des pulsars. 

L'evolution d'une etoile est liee a sa masse initiale. Plus I'etoile est massive, moins elle 
passe de temps dans la sequence principale a synthetiser des noyaux d'helium a partir des 
protons et des electrons. On pourra se reporter au diagramme d'Hertzsprung et Russell 



pour parcourir revolution stellaire (Russell 1921). Quand les noyaux d'hydrogene viennent 



a manquer, la pression radiative qui compense I'attraction gravitationnelle chute. Le coeur 
de I'etoile se contracte et augmente en temperature, permettant la combustion de I'helium. 
Ce phenomene s'accompagne de I'expansion des couches externes. Pour les etoiles de types 
O et B, superieures a 8 - 10 masses solaires, la fusion thermonucleaire d'elements de plus en 
plus lourds va continuer jusqu'au fer. La nucleosynthese de tout autre element plus lourd ne 
produit pas d'energie, au contraire elle en consomme. A ce stade la pression radiative creee par 
I'energie nucleaire s'arrete. Le coeur de I'etoile s'effondre creant une onde de choc qui expulse 
les couches superieures de I'etoile dans le milieu interstellaire. Ce phenomene est appele 

Si 



supernova (|Baade and Zwicky, 1934) de type Ib/c, ou II (Cappellaro and Turatto 



le coeur restant excede la masse critique de Chandrasekhar ( Chandrasekhar 1935 



2000) 



masse a 



laquelle la pression de degenerescence des electrons s'oppose a la force gravitationnelle, les 
electrons et les protons vont se combiner pour former des neutrons. L'effondrement s'arretera 
quand la pression des neutrons degeneres compensera la pression gravitationnelle, formant 
ainsi une etoile a neutrons. 



1.2.1 Etoiles a Neutrons 

Les etoiles a neutrons sont I'un des objets les plus compacts que Ton connaisse dans 
notre univers (apres les trous noirs). Leur structure (composition interne, masse, et rayon) 
est encore mal connue et depend de I'equation d'etat de la matiere. Cette discipline constitue 
un sujet d'etudes important en physique nucleaire. L'analyse theorique de la structure des 
etoiles a neutrons est basee sur la chronometrie des pulsars et les observations optiques et X 
des etoiles elles-memes. 



La figure 1.2 (gauche) presente un des modeles d'etoile a neutrons pour une masse de 
1.4 masse solaire et un rayon de 15 km. Entre la surface de I'etoile et son centre, la densite 
volumique varie de lO'' a lO^^gcm^^. L'etoile est composee d'une croute surfacique rigide 
d'une epaisseur de 1 km formee essentiellement de noyaux de fer et d'un noyau de neutrons 
superfiuides. La separation entre ces deux parties se situe pres de p = 4.3 x 10^^ gcm^^. La 
plupart des modeles de matiere dense predisent un deconfinement des quarks et une possible 
apparition de matiere exotique (quark etrange) au-dessus d'une densite de ~ 5 x 10^'^gcm^^, 
la densite de la matiere nucleaire froide. Cette structure composee de differents etats de la 
matiere pourrait expliquer a la fois I'observation de perturbations continues de la rotation du 
pulsar appelees "timing noise", ainsi que I'acceleration brusque de la periode de rotation des 



pulsars que Ton nomme "glitch" (Anderson and Itoh 1975; Janssen and Stappers 2007). Ces 
evenements sont probablement la consequence du transfert du moment angulaire de la croiite 
solide, qui tourne a la periode mesuree du pulsar, jusqu'a la matiere au centre de I'etoile 
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tournant plus vite. Les glitch sont suivis par un retablissement de la frequence de rotation, 
qui est la consequence d'un nouvel equilibre entre les deux composants. 



La figure 1.2 (droite) presente la masse des etoiles a neutrons mesuree a partir de pulsars 
presents dans des systemes binaires, revelant une masse moyenne autour de 1.4 masse solaire 
(masse de Chandrasekhar). L'intervalle de masse predit par les modeles theoriques d'etoiles 
a neutrons s'etend de 0.2 a 3.0 masses solaires, et depend des modeles d'equation d'etat de 
la matiere nucleaire. Par la suite, nous prendrons une masse d'etoile a neutrons Mns de 1.4 
Mo. 



Pour plus de details sur les etoiles a neutrons, on pourra se reporter a I'article de Lattimer 



and Prakash] ( |2007| ) et au cours en ligne de |Gourgoulhon| ( |2004 - 2005D 




DitHir 



Radius (hm) 



(a) structure possible d'une etoile a neutrons. 
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(b) Masses d'etoiles a neutrons. 



Figure 1.2 - A : Structure possible d'une etoile a neutrons. B : Masses mesurees d'etoiles d 
neutrons localisees dans des systemes binaires. Pour chaque region, la moyenne des masses 
est indiquee en pointille, tandis que la moyenne ponderee est representee par la ligne discon- 



tinue. La figure est extraite de Lattimer and Prakash (2007). 



1.2.2 Rotation 

Pendant la phase de supernova, le systeme 'coeur de I'etoile massive / etoile a neutrons' 
passe en quelques secondes de Rq ~ 10^ km (rayon moyen) a R^s ~ 10 km. Si Ton assimile 
I'etoile a une sphere, son moment d'inertie I est egal a 2/5 x M^sR^s > ~ IC^^gcm^. En 
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1.2. Evolution de la Rotation et Proprietes Intrinseques 



TIS -,2 



(1.1) 



supposant que le moment angulaire lO avec D. = 2n/V et la masse du systeme soient en 
partie conserves pendant I'effondrement, la periode initiale de rotation de I'etoile a neutrons 
Pns sera determinee par la periode de I'etoile initiale Pq, avec : 

P ~PAf — 
f ns — "^01 n 

Ainsi, pour une etoile initiale de periode Pq lO^s (notre soleil tourne sur lui-meme en 
2 X lO^s ~ 27 jours), la periode de rotation de I'etoile a neutrons est de I'ordre de 10ms. Ce 
chiffre est en accord avec la periode observee des jeunes pulsars tels que le Crabe (33 ms), 
PSR J0205+6449 (64ms), ou PSR J1833-1034 (61.9ms). 



D'autre part, comme I'ont observe Richards and Comella (1969) pour le Crabe, la periode 



de rotation des pulsars augmente au cours du temps. En d'autres termes, ils ralentissent, 
indiquant un transfert de leur energie cinetique de rotation. Si Ton considere le pulsar comme 
un rotateur rigide avec un moment d'inertie I, I'energie cinetique disponible est donnee par 
\lO.^, avec CI la vitesse angulaire. En consequence, I'energie de rotation perdue par seconde 
est : 

^'^-"'^ -laa = 471^1 (P/P^) erg s-\ (1.2) 



"™^-"di^2^^ 
Avec un ralentissement P de 4.2 x 10^^^ 

d'inertie I de ~ IC^^gcm^, le taux d'energie de rotation perdue par le pulsar du Crabe est 



s s ( Manchester et al 



2005 ) et adoptant un moment 



de 4.6 X 10^^ ergs ^rl Les lignes pour les differentes valeurs de E en fonction de (P,P) sont 



dessinees sur la figure 2.2 



1.2.3 Champ Magnetique 

Durant I'effondrement gravitationnel le flux magnetique 4) oc BR^ est egalement conserve 
a I'interieur de I'etoile. Le champ magnetique polaire a la surface de I'etoile a neutrons B^s 
est ainsi determine par le champ magnetique de I'etoile initiale Bq, par : 



B ~ 

L'ns — 



Bo( 



.2 



(1.3) 



Utilisant les valeurs standards precedentes pour le rayon et adoptant un champ magnetique 
de 100 G (10~^ Tesla) pour I'etoile initiale, le champ magnetique resultant est de I'ordre de 
10^^ G. Ostriker and Gunn (1969) montrent qu'une etoile a neutrons de rayon R^s tournant a 



la vitesse angulaire D. dans le vide, avec un champ magnetique dipolaire non parallele a I'axe 
de rotation rayonne de I'energie. La puissance rayonnee est donnee par : 



dW 



q4 

3^RnsBnsSm'cX 



(1.4) 



ou (X est Tangle entre I'axe du dipole et I'axe de rotation (voir la figure 1.6) et c est la vitesse 
de la lumiere. En reliant les equations 1.2 et |1.4 , on pent estimer, avec les seules observables 
(P la periode de rotation et P sa derivee), le champ magnetique a la surface de I'etoile : 



Bns =3.2 X lO^^(PP)^/^ G. 



(1.5) 



2. 1 erg = 6.24150974 x 10" electron-volt 
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Pour le pulsar du Crabe, le champ magnetique de surface obtenu est de I'ordre de 4 x 10^^ G, 



en accord avec le champ moyen estime apres I'effondrement. 
valeurs de B^s en fonction de (P,P) sont dessinees sur la figure 



niveau du cylindre de lumiere (voir la section 1.4.1) est quant a lui decrit par : 



lies lignes pour les differentes 
2.2 Le champ magnetique au 



= (^^i^)V2 ^ 2.943 X 10«(PP-^)^/^ G. (1.6) 

1.2.4 Indices de Freinage et Age du Pulsar 

Rappelons que les resultats precedents sont verifies dans I'approximation dipolaire du 
pulsar. Nous allons voir que la loi qui determine le ralentissement de la rotation est legerement 



differente de celle observee dans le cas d'un dipole magnetique. A partir des equations 1.2 et 



1.4 la variation de la vitesse angulaire du pulsar est donnee par 



blc^ 

Plus generalement, si Ton suppose que le pulsar est forme avec une grande vitesse angulaire, 
on admet que la vitesse angulaire suit une simple loi de puissance d'indice n de la forme : 

a = -kO^, (1.8) 

ou k est une constante et n est defini comme I'indice de freinage. En derivant I'expression 
precedente, on pent decrire I'indice n a partir de la periode de rotation du pulsar et ses 
derivees : 

O.CL ^ vv , , 

n=^=2-^. (1.9) 

Q2 p2 ^ ^ 

Notons qu'il est tres difficile de mesurer I'indice de freinage. Cela demande de connaitre par- 
faitement la seconde derivee, laquelle est dominee par les instabilites du pulsar. Quelques 
valeurs ont ete mesurees pour les pulsars du Crabe (n = 2.525 ± 0.005), de PSR B 1509— 58 
(n = 2.8 ± 0.2) et de PSR B0540-69 (n = 2.01 ± 0.02). On remarquera I'ecart entre les 
indices de freinages mesures et la valeur theorique n = 3 dans le cas unique de radiation 
par un dipole magnetique. Cependant, plusieurs effets n'ont pas ete pris en compte dans 
ce modele tres simplifie : la dissipation de I'energie par remission d'ondes gravitationnelles, 
la deformation du champ magnetique B par la rotation du pulsar, les effets de glitch et de 
timing noise, ainsi que le vent d'electrons relativistes cree par le pulsar qui heurte le milieu 



interstellaire (voir 1.4.1). 



En integrant I'equation 1.8 et en supposant que la periode P est grande devant la periode 



initiale Pq, on pent deduire I'age caracteristique du pulsar par : 



)n-l 



Toujours dans I'approximation du dipole magnetique (n = 3), I'age caracteristique pent etre 
decrit tres simplement par t = P/2P. Dans le cas du pulsar du Crabe cree par la supernova en 
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1054 (955 ans) avec une periode P de 33 ms et une deriveePde4.2x^0-^^l' age caracteristique 
obtenu est de ~ 1260 ans, ce qui est legerement different (mais raisonnable) si on utilise la 
valeur experimentale n = 2.525 (1600 ans). Dans le cas du pulsar PSR J0205+6449, nous 
verrons au chapitre [6] que I'age estime de cette fagon est en desaccord avec la supernova 
associee a cet objet. 



1.3 Premieres Observations de TEmission Pulsee y 



Un an apres la decouverte du premier pulsar en radio par Hewish et al (1968), les astro 



nomes entamerent une recherche dans les autres domaines du spectre electromagnetique. En 



1969, Cocke et al (1969) publient la decouverte de pulsations optiques du Crabe, tandis que 



Pritz et al (1969) et Bradt et al (1969) annoncent sa decouverte en X etablie grace a une 



fusee placee dans la haute atmosphere. Le rayonnement electromagnetique de haute energie 
tels que le rayonnement X ou y, ainsi que certaines longueurs d'ondes en infra-rouge sont 
absorbes par I'atmosphere terrestre et n'atteignent pas le sol. II est done indispensable, pour 
observer ces domaines d'energie, d'utiliser des instruments situes dans la haute atmosphere 
ou au-dela. 



Third EGRET Catalog 

E > 100 MeV 

+90 




^ Active Galactic Nuclei 

O Unidentified EGRET Sources 



■ Pulsars 
A LMC 

• Solar FLare 



Figure 1.3 - Cette carte indique la 'position des sources du troisieme catalogue d'EGRET 



(Hartman et al 1999). Les sources sont placees suivant leurs coordonnees galactiques et sont 
identifiees suivant un code de forme. Cette figure est extraite de Thompson ^2008 ). 



27 



Chapitre 1. Les Pulsars 



L'astronomie gainina|^ debute en 1972 avec le Small Astronomy Satellite (SAS-2), qui 
met en evidence trois sources ponctuelles le long du plan galactique, toutes trois identifiees 
comme des pulsars. Deux d'entre elles sont mises en evidence par leur coincidence spatiale et 



la pulsation de leur emission avec les pulsars radio du Crabe et de Vela (Kniffen et al 1974 



Thompson et al, 1975), tandis que la troisieme source nommee Geminga n'est identifiee qu'en 



1992 en tant que pulsar, par la decouverte d'une periodicite de 237 ms dans les photons X 
(Halpern and Holt 1992) et les photons y ( Bignami and Caraveo[ 1992). Trois ans plus tard, 
I'Agence Spatiale Europeenne met en orbite le Cosmic Ray Satellite (COS-B). L'instrument va 



reveler la presence de nombreuses sources ponctuelles dans le plan galactique (Hermsen et al 



1977), et cartographie pour la premiere fois le fond diffus y qui resulte de la desintegration 
de rfi produits par I'interaction des rayons cosmiques avec le gaz interstellaire. 

Cependant, c'est en 1991 que l'astronomie y connait un essor important. La NASA lance 
le Compton Gamma Ray Observatory (CGRO) qui compte a son bord quatre instruments : 
BATSE (Burst And Transient Source Experiment) dedie a I'observation des sursauts y entre 
20 et 100 keV, OSSE (Oriented Scintillation Spectrometer Experiment) utilise pour I'etude 
des raies de desintegration des noyaux radioactifs dans les vestiges de supernovae, COMPTEL 
(imaging COMPton TELescope) dedie a I'etude des sources entre 1 et 30 MeV, et EGRET 
(Energetic Gamma Ray Experiment Telescope) pour I'etude des sources entre 30 MeV et 
30 GeV. Ce dernier, precurseur du telescope Fermi, presente une meilleure resolution angulaire 
et une plus grande surface de collection que ses predecesseurs, et dresse ainsi le premier 
catalogue des sources y de hautes energies (Hartman et al 1999). La figure 1.3 presente la 
position dans le ciel des 271 sources du troisieme catalogue d'EGRET, dont 94 associations 
avec des noyaux actifs de galaxie, un sursaut solaire, le nuage de Magellan, et 5 (+1) pulsars]^ 
II est important de souligner que 170 sources n'ont pas ete identifiees en raison de I'incertitude 
sur leur position, notamment au niveau du plan galactique. Ce sujet constitue I'un des enjeux 
majeurs de I'experience Fermi. 



1.3.1 Pulsars CGRO 1991 - 1999 



Les telescopes a bord du satellite CGRO ont identifie au total 7 pulsars avec un haut 
niveau de confiance. Le telescope EGRET a detecte remission y pour 6 pulsars (9 possibles), 
comprenant 5 pulsars radio et le pulsar Geminga dont la pulsation radio n'a jamais ete 
observee. Quant au septieme pulsar PSR B1509— 58, il a ete detecte en dessous de 5 MeV 
par le detecteur COMPTEL (Matz et al, 1994). Le tableau [171] reference les principales pro- 
prietes des pulsars detectes par les deux instruments COMPTEL et EGRET, ainsi que 3 
pulsars detectes par EGRET avec une significativite plus faible, PSR J 1048— 5832 (chapitre 

PSR J0659+1414 (chapitre [9]), et le pulsar milliseconde J0218+4232 ([Kuiper et "al[|2000[ 



3. Le domaine y est defini au-dela de 100 keV (des photons jusqu'a la centaine de TeV sont detectes). Etant 
donne le vaste domaine, les methodes de detection sont differentes suivant I'energie du photon y . On distingue 
le domaine de basse energie entre 100 keV et 30 MeV (detection du y par effet photoelectrique), le domaine 
de haute energie entre 30 MeV et 30 GeV (detection du y par conversion e*), et le domaine de tres haute 
energie au-dessus de 30 GeV (detection de la lumiere Tcherenkov emise par la cascade creee par le photon y 
primaire). 

4. Le pulsar PSR B1951+32 n'est pas inclut dans le 3®"^'= catalogue d'EGRET a cause de son flux trop faible, 
mais est represente sur la carte. 
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chapitre IqI). Ly est la luminosite y, donnee par : 



(1.11) 



ou d est la distance du pulsar depuis la Terre (voir la section 2.2), Ft^obs est le flux energetique 



observe depuis la Terre et moyenne sur toute la phase, et fo(a,CE) est le facteur de correction 



qui prend en compte la geometrie du faisceau donne par ([ Watters et al[ |2009[ ) 



2jFY(a;CE,*)(i(t) 



(1.12) 



ou F-y(a; C,4^) est le flux rayonne en fonction de Tangle d'observation C par rapport a I'axe 
de rotation (ligne de visee) et la phase du pulsar (4)), tandis que fa est le rapport entre 
remission y globale sur tout le ciel et le flux prevu moyenne sur toute la phase pour une 
courbe de lumiere observee depuis la Terre. Enfln, I'efficacite y n otee r\ est le rapport entre 
la luminosite (equation |l.ll[) et le ralentissement E (equation 1.2). 



Nom 


Periode 




E 


Fe 




d 


Ly 


11 




(s) 


(ky) 


(ergs"') 


(erg cm"^ 


s-l) 


(kpc) 


(ergs"') 


(E>1 eV) 


Crab, J0534+2200 


0.033 


1.3 


4.5 X lO^s 


1.3 X 10 


-8 


2.0 


5.0 X 10^5 


0.001 


B1509-58 


0.150 


1.5 


1.8 X 10^7 


8.8 X 10- 


10 


4.4 


1.6 X 1035 


0.009 


Vela, J0835-4510 


0.089 


11 


7.0 X 10^6 


9.9 X 10 


-9 


0.3 


8.6 X 10^3 


0.001 


J1709-4429, B1706-44 


0.102 


17 


3.4 X 10^'' 


1.3 X 10 


-9 


2.3 


6.6 X 1034 


0.019 


J1952+3252, B1951+32 


0.040 


110 


3.7 X 10^6 


4.3 X 10- 


10 


2.5 


2.5 X 1034 


0.007 


Geminga, J0633+1746 


0.237 


340 


3.3 X 10^4 


3.9 X 10 


-9 


0.16 


9.6 X 1032 


0.029 


J1057-5226, B1055-52 


0.197 


530 


3.0 X 10^4 


2.9 X 10- 


10 


0.72 


1.4 X 1033 


0.048 


J1048-5832, B1046-58 


0.124 


20 


2.0 X 10^6 


3.7 X 10- 


10 


2.7 


2.6 X 1034 


0.013 


J0659+1414, B0656+14 


0.385 


100 


4.0 X 10^4 


1.6 X 10- 


10 


0.3 


1.3 X 1032 


0.003 


J0218+4232 


0.002 


460000 


2.5 X 10^5 


9.1 X 10- 


1 1 


2.7 


6.4 X 1033 


0.026 



Tableau 1.1 - Parametres intrinseques et mesures des pulsars et des possibles candidats 
detectes par les telescopes COMPTEL et EG RET a bord du satellite CGRO. P est la periode 
de rotation, x est I' age caracteristique (voir 1.2.4). ^ est le taux de perte d'energie. est 



le flux en energie observe a haute energie (E>1 MeV). d est la distance entre le pulsar et 
la Terre . Ly est la luminosite y, tandis qu'r[ est I'efficacite de conversion entre I' energie de 
rotation E et la luminosite. 



Une grande partie de I'information que Ton connait sur les pulsars derive des proprietes 
temporelles. Pourtant, les parametres E, Bns, t, ne repondent pas a la question de I'origine du 
rayonnement electromagnetique. C'est a partir des proflls pulses et des spectres en energie 
resolus en phase que Ton va pouvoir identifier les processus d'emission qui produisent le 
rayonnement pulse (voir 1.4.3). Les sous-sections suivantes presentent les profils pulses et les 
spectres en energies des pulsars detectes par la mission CGRO. 
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1.3.2 Courbes de Lumiere 



La figure L4 montrent les courbes de lumierej^ des 7 pulsars detectes par la mission 



CGRO pour 5 intervalles en energies : de haut en bas, radio, optique, rayonnement X mou 
(<lkeV), X dur/y mou (lOkeV - IMeV), et rayonnement y (au-dessus de 100 MeV). Les 
particularites de ces courbes de lumieres sont, d'une part leurs profils differents en fonction 
de la longueur d'onde observee, ce qui illustre une geometrie et des mecanismes d'emission qui 
dependent de I'energie, emission radio coherente, rayonnement X thermique, et rayonnement 
y non thermique. D 'autre part, pour la plupart des pulsars, les courbes de lumieres affichent 
une structure en double pics couvrant environ 50% de la rotation, ce qui pourrait signifier un 
large cone d'emission y. PSRs J1709-4429 (B1706-44) et J1057-5226 (B1055-52) semblent 
avoir un profile y plus complexe. Ainsi, comprendre la difference entre les profils radio (la 
plupart sont formes d'un pic) et les profils y constitue un defi pour la mission Fermi. Notons 
egalement que PSR B1509— 58 qui a le champ magnetique le plus intense des 7 pulsars, n'a 
pas ete observe par EGRET. 

1.3.3 Spectres Multi-longueurs d'Onde 



La figure \L5\ presente les spectres multi-longueurs d'onde, de la radio au rayonnement y 
pour les 7 pulsars detectes par la mission CGRO. La caracteristique commune est une lumi- 
nosity maximale dans I'intervalle d'energie y, ainsi qu'une reelle difference entre remission 
radio (processus coherent) et remission de haute energie (acceleration de particules chargees), 
notamment pour les pulsars du Crabe et de Vela. Pour les pulsars de Geminga, B1055— 52, 
et Vela, le spectre X indique une emission thermique qui provient de la surface de I'etoile 
a neutrons, tandis que le spectre y pour les 7 pulsars est relativement plat (indice spectral 
autour de 2), et montre au moins pour les pulsars EGRET une coupure autour du GeV. Seul 
PSR B1951-I-32 n'en possede pas. Rappelons egalement qu'une detection pulsee du pulsar du 



Crabe, au-dessus de 25 GeV, a ete observee par I'experience Cherenkov MAGICf] ( [MAGIC 



2008), soulignant I'importance de ces coupures dans les mo deles des pulsars (voir la section 



1.4.3). Une limite superieure de remission pulsee avait ete determinee quelques annees au- 



paravant par la collaboration CELESTE indiquant une coupure dans le spectre y d'EGRET 



(de Naurois et al 2002). 



1.4 Consequences et Modeles d'Emission Gamma 

Les astronomes comprirent rapidement que les pulsars ne pulsaient pas, mais que I'ef- 
fet de pulsation provenait de la rotation d'un faisceau d'onde electromagnetique emis dans 
I'environnement du pulsar. Cependant, comment expliquer une emission pulsee sur tout le 
domaine d'energie et notamment a haute energie, comme pour les pulsars du Crabe et de 
Vela? Quels sont les processus physiques mis en jeu pour produire des photons jusqu'au 
GeV ? Ou et comment les particules chargees sont-elles accelerees dans la magnetosphere du 

5. La courbe de lumiere ou profil pulse est rhistogramme du nombre de coups en fonction du temps, modulo 
la periode du pulsar. 

6. Major Atmospheric Gamma-ray Imaging Cherenkov telescope : http ://magic. mppmu.mpg.de/ 
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Crab B1 509-58 Vela B1 706-44 B 1951 +32 Geminga B1 055-52 
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Figure 1.4 - Courbes de lumieres pour 5 intervalles en energie des 7 pulsars detectes par 
la mission CGRO : radio, optique, rayonnement X mou (<1 keV), X dur / y mou (lOkeV - 
1 MeV), et rayonnement y (> 100 MeV). Les pulsars sont ranges de gauche a droite par age 
croissant. Chaque case montre une rotation de I 'etoile a neutrons. Cette figure est extrait de 
la revue de Thompson ^2008) . 



pulsar ? Est-ce que les processus sont les memes pour tous les pulsars ? Quel est le rapport 
entre les pulsars radio emetteurs y tels que le Crabe ou Vela, et les pulsars y silencieux en 
radio tel que Gemingaj^? 

Seule une petite partie du faisceau electromagnetique balaye notre ligne de visee. La conse- 
quence est qu'il n'est pas simple de construire la geometrie complete de remission a partir 



7. Dans la litterature, le terme radio-loud est employe pour definir les pulsars emetteurs radio et y, tandis 
que le terme radio-faint est employe pour les pulsars y dont remission radio n'a pas ete mis en evidence. 
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Figure 1.5 - Spectres multi-longueurs d'onde des 7 pulsars detectes par la mission CGRO 
Thompson\ ( fgOOgp . 
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des observations (voir la figure 1.4). Theoriquement, I'idee generale est que les emissions 
de la radio jusqu'aux hautes energies sont produites dans differentes regions de la magneto- 
sphere du pulsar. Les ondes radio polarisees sont emises depuis la calotte polaire par emission 
synchrotron, zone au-dessus des poles magnetiques de I'etoile. Quant aux photons de haute 
energie (optique, rayonnement X, et rayonnement y), ils sont produits principalement par le 
rayonnement de courbure et la diffusion inverse Compton des electrons (et positons) accele- 
rees le long des lignes de champ magnetique. Cette propriete est un reel avantage, puisque les 
particules accelerees rayonnent leur energie tres vite et sont liees aux sites d'acceleration. De 
plus, la luminosite des 6 pulsars detectes par EGRET atteint son maximum dans le domaine 



des y (voir la figure 1.5), laquelle est une fraction de 0.1% a 10% de I'energie cinetique de 
rotation (~ 10^"* — 10^^ pour la radio). Dans ce sens, I'etude des mecanismes d'acceleration des 
particules passe principalement par I'observation du rayonnement y qui en resulte. Ainsi, I'un 
des enjeux du telescope Fermi est de contraindre la geometrie et les mecanismes des zones 
d'acceleration, en etudiant en details les spectres resolus en phase d'une large population de 
pulsars y. 

Les trois sections suivantes presentent la theorie d'emission standard des pulsars et les trois 
modeles d'emission y, lesquelles se differencient par I'emplacement de I'origine du faisceau, 
le modele de la calotte polaire ou Polar Cap (PC), le modele Slot Gap (SG) ou Two Pole 
Caustic (TPC) et le modele de la cavite externe ou Outer Gap (OG). 

1.4.1 Modele Standard 

Comme nous I'avons vu precedemment, il est parfaitement raisonnable de penser que les 
pulsars sont des etoiles a neutrons en rotation rapide, fortement magnetisees, formees dans 
I'explosion d'une supernova. Ceci est le principe de base des theories d'emission des pulsars. 
Nous introduisons ici certaines proprietes importantes sur la magnetosphere des pulsars, ou le 
champ magnetique et le moment angulaire vont jouer un role preponderant dans sa structure. 
Notons que ces proprietes sont decrites dans le cas relativement simple ou I'axe de rotation 
du pulsar est aligne avec son axe magnetique (modele du pulsar aligne). 

Si Ton considere que I'etoile a neutrons est un parfait conducteur, le champ magnetique 
dipolaire B tournant avec une vitesse angulaire II induit pres de I'etoile un champ electrique 
E : 

E+-(dxf)xB=0. (1.13) 
c 



Goldreich and Julian (1969) montrent dans le cas ou I'axe de rotation est aligne a I'axe 
magnetique, que les forces electromagnetiques en dehors de I'etoile dominent completement 
les processus physiques, et sont suffisamment intenses pour arracher les particules de la surface 
de I'etoile et creer une magnetosphere composee d'un plasma charge. La densite de charge, 
appelee egalement la densite de charge de Goldreich- Julian pcj, est definie par : 

1 ^ ^ d-B 

tandis que la densite de particules est donnee par : 

n = 7 X 10"^ B^ P"^ particules cm"^ (1.15) 
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Figure 1.6 - Illustration des modeles Polar Cap, Slot Gap et Outer Gap. d represents I'axe 
de rotation, tandis que B represente I'axe magnetique. L'inclinaison magnetique a est I'angle 
entre D. etB. Les zones colorees representent les differents sites d' acceleration et d' emission 
y, avec les surfaces de densite de charge nulle (pc] representees par les lignes bleues. 



avec la composante axiale du champ magnetique en gauss, et P la periode de rotation. 

Le champ electrique intense cree influence la region comprise entre le pulsar et le cylindre 
de lumiere tc = c/Cl, distance ou les lignes de champ magnetique et les particules en co- 



rotation avec I'etoile a neutrons atteignent la vitesse de la lumiere. La flgure |1.6| montre la 
geometrie de la magnetosphere et la zone d'emission des modeles y dans le cas ou I'axe de 
rotation et I'axe magnetique ne sont pas alignes. La zone fermee (partie grisee) contient le 
plasma en corotation, tandis que la zone ouverte permet au plasma de s'ecouler le long des 
lignes de champ, et de s'echapper dans le milieu interstellaire creant un "vent" relativiste|^ 
Les zones au niveau de la surface de I'etoile d'ou proviennent les lignes ouvertes sont deflnies 
comme les calottes polaires, lesquelles sont delimitees par la derniere ligne de champ fermee, 
soit par un cone de demi-angle d'ouverture 9cp : 

sinOep) = (^)V2, (1.16) 

8. Ce vent en contact avec le milieu interstellaire pent creer une source : une nebuleuse a vent de pulsar 
(PWN : Pulsar Wind Nebula). 
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avec Rns le rayon de I'etoile a neutrons. En supposant que le rayon du cylindre de lumiere 
est tres superieur a celui de I'etoile a neutrons, le rayon de la calotte polaire tcp a la surface 
de I'etoile est donnee par : 



I cp 



i^ns»-'cp — 



(1.17) 



Le flux de particules chargees accelerees le long des lignes de champs ouvertes au niveau de 
la surface de I'etoile est donne par : 



' ^cp — 



(7tr^p)(^: 



2ec ' 



(1.18) 



avec Bs le champ magnetique a la surface de I'etoile. La difference de potentiel entre le pole 
de I'etoile et la derniere ligne de champ fermee est proportionnelle a la perte de I'energie de 
rotation |1.2[ et pent etre exprimee a partir des equations |1.13| et |1.16| par : 



AO = 



E- ds 



(1.19) 



Les grandeurs donnees aux equations 1.16 , 1.17[ et 1.19 deflnissent les principales ca- 
racteristiques du modele standard des pulsars. En appliquant ces equations au pulsar du 
Crabe, le demi-angle d'ouverture de la calotte polaire obtenu est de Q^p = 4.8 ° avec un rayon 
0.84 km, tandis que la difference de potentiel atteinte entre le centre et le bord de la 



' cp 



calotte est de AO ~ 2 x 10^^ V. 



1.4.2 Les Modeles d'Emission Gamma 



D'apres Goldreich and Julian (1969), le champ electrique induit par le champ magnetique 



dipolaire de I'etoile a neutrons est si intense que la magnetosphere se remplit de particules 



chargees, arrachees a la surface de I'etoile. L'equation 1.14 implique a I'interieur de la ma- 
gnetosphere une densite de charge en corotation avec le champ magnetique et qui ecrante 
la composante du champ electrique parallele au champ magnetique (E ■ B = 0). Cependant, 
il existe des zones deflnies par O ■ B = ou cette densite est nuUe, et ou le champ elec- 
trique n'est pas ecrante (E ■ B 7^ 0) (voir la flgure 1.6). Cette propriete joue un role essentiel 



puisque a I'interieur de ces cavites les particules chargees sont accelerees a des energies ultra- 
relativistes. Les processus physiques mis en jeu pour creer des photons de haute energie sont 
le rayonnement de courbure, remission synchrotron, et la diffusion inverse Compton. Pour 



plus details sur ces processus, on se reportera a I'article de Rudak et al (2002). 

Les tentatives pour modeliser les accelerateurs de particules dans la magnetosphere du 
pulsar sont reparties en deux scenarios. Les modeles Polar Cap et Slot Gap placent la source 
d'emission au-dessus des poles magnetiques jusqu'a quelques rayons de I'etoile a neutrons, 
tandis que le modele Outer gap propose une emission pres du cylindre de lumiere. 



Polar Cap et Slot Gap 



Le modele Polar Cap fut propose par Sturrock ( 1971 ) pour expliquer remission radio. II 
suppose que les particules arrachees a la surface de I'etoile sont accelerees dans une petite 
zone au-dessus du pole magnetique par le champ electrique (voir la flgure 1.7). A partir 
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de rapproximation p = (l/47t)V^O et de Tequation 1.14 
de particules chargees est donne par 



CD 



le potentiel responsable du flux 
(2nBs/c)h77~^vec h ~ r^p, la hauteur de cette 
region. Sturrock ( 1971 ) demontre que ce potentiel est sufRsant pour accelerer les particules 
a des vitesses relativistes (y < 10'') le long des lignes de champ magnetique. Ces particules 



primaires (e.g. electrons) rayonnent des photons y par rayonnement de courbure (Ruderman 
and Sutherland! 1975), dont I'energie caracteristique est : 



3 y^ 
= -cTi — eV, 

2 Qrr 



(1.20) 



avec y le facteur de Lorentz et Pcr le rayon de courbure de la ligne de champ. Lorsque ces 
photons se deplacent a travers le champ magnetique intense, ils creent des paires electron- 
positon (y + B ^ e+e^)|^ Ces nouvelles particules chargees sont elles-memes accelerees et 
produisent a leur tour des y qui se materialisent en paires e^ et ainsi de suite. Ces cascades 
electromagnetiques s'echappant le long des lignes de champ ouvertes, alimentent le plasma 
de la magnetosphere et ecrantent le champ accelerateur Ej|, limitant la zone acceleratrice. Les 
paires secondaires (voire tertiaires) produites le long de I'axe magnetique, a basse altitude, 
sont a I'origine de remission radio observee. 
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Figure 1.7 - Schema de la magnetosphere d'un pulsar dans le cas du modele de Goldreich- 
Julian. La figure de droite illustre les particules chargees arrachees a la surface de I'etoile. 
Les particules accelerees par un champ electrique intense emettent des photons y par rayon- 
nement de courbure. 



L'emission observee a haute energie est due, quant a elle, a la diffusion inverse Compton 
des paires e^*" sur les photons infra-rouges (Daugherty and Harding 1986). Ce mecanisme est 



tres efiicace dans le perte d'energie des particules chargees, produisant des photons y plus 



9. Un photon d'energie Ey se propageant selon un angle \\> par rapport a la direction du champ magnetique, 
pent se materialiser en paire e* si EySin(t|)) > 2meC^. 
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energetiques que dans le processus de rayonnement de courbure. Seulement, etant donne 
que remission est situee au-dessus des poles, le champ magnetique est intense. Le modele 
predit ainsi une coupure de type super-exponentielle dans le spectre pulse du pulsar en rai- 
son de I'absorption des photons y par le processus de creation de paires. Notons finalement 
que ce modele reproduit correctement le spectre resolu en phase observe pour le pulsar de 
Vela, mais sous certaines conditions, telles qu'une zone d'acceleration a trois rayons de I'etoile. 



Pour attenuer les effets du champ magnetique sur les photons y, Arons (1983) a examine 



un autre site d'acceleration compris entre la calotte polaire et le cylindre de lumiere. Le 
site est relativement fin et confine le long de la derniere ligne de champ ouverte (zone rose 



sur la figure 1.6). Ce type d'accelerateur est appele "Slot Gap" (SG) ou "Two-Pole Caustic" 



(TPC) (Dyks and Rudak 2003; Harding 2006). II reproduit correctement les proprietes de 
polarisation, les profils pulses et les spectres resolus en phase pour le pulsar du Crabe entre 
0.3 keV et 10 GeV, ainsi que le spectre a haute energie du pulsar milliseconde PSR J02 18+4232 



(Kuiper et al 2000). 



Modele de la Cavite Externe : Outer Gap 



Dans le modele Outer Gap (Cheng et al 1986 Romani and Yadigaroglu 1995), remission 



y se produit pres du cylindre de lumiere, entre la ligne de densite de charge nuUe (O ■ B =0), 



la derniere ligne de champ fermee, et le cylindre de lumiere (section 1.4.1). La figure 1.6 



illustre la zone d'emission (zone bleue) dans le cas du modele OG. HoUoway (1973) appuie 
I'idee que les charges qui s'echappent de cette region ne sont pas remplacees, et laissent une 
cavite vide separee par des densites de charges opposees. La consequence est I'apparition d'un 
champ electrique parallele aux lignes de champ magnetique. 

Les electrons et les positons se propageant le long des lignes de champs pres de cette cavite 
emettent des photons par rayonnement de courbure ou par diffusion inverse Compton sur des 
photons de plus basse energie. Ces photons crees vont penetrer dans la cavite, et a leur tour 
se materialiser en paires e^. La difference de potentiel presente dans la cavite va accelerer 
les particules chargees, lesquelles vont emettre des photons de haute energie. Comme dans 
le modele PC, ce fiux de particules chargees va delimiter la cavite acceleratrice et maintenir 
un intense courant charge a I'origine de remission y. Generalement, ces cavites externes ne 
sont pas prevues pour etre associees a remission radio, bien que le Crabe soit peut-etre 
une exception. De plus, le champ magnetique local etant beaucoup plus faible qu'aux poles, 
I'attenuation spectrale a haute energie est moins importante que dans le modele PC. 

Selon ce modele, on distingue deux types de pulsars, associes au Crabe et a Vela. Pour 
les pulsars de type Crabe relativement jeunes, la cavite produit des photons au GeV par 
rayonnement de courbure, puis produit des electrons et des positons, lesquels rayonnent par 
emission synchrotron et par diffusion inverse Compton. Dans le cas des pulsars de type Vela, 
les paires sont creees dans la cavite, et accelerees a des energies relativistes dans le sens 
oppose au vent. Ces particules chargees primaires vont produire des photons y par diffusion 
inverse Compton sur des photons infrarouge (IR). Les photons y suffisamment energetiques 
vont a leur tour se materialiser en paires e^, par interaction avec la composante IR. Suite a 
cela, remission synchrotron de cette deuxieme population de particules chargees va donner 
un faisceau de photons X et y, qui vont interagir pour produire une troisieme population de 
paires moins energetique. Les photons IR seraient crees par remission synchroton de ces 
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paires tertiaires. Notons egalement que ce modele predit une diffusion inverse Compton des 
photons IR sur les particules primaires ; cette diffusion serait par ailleurs a I'origine d'une 
emission au TeV qui n'a pas ete encore observee. 

Finalement, les deux modeles sont assez differents. Dans le cas des pulsars de type Crabe, 
les paires primaires perdent la quasi totalite de leur energie par le rayonnement de cour- 
bure, tandis que dans le cas des pulsars de type Vela, les paires perdent leur energie par 
emission synchrotron. La consequence est que la luminosite y predite dans les deux modeles 
est differente. Pour les pulsars de type Crabe, elle est donnee par : 

= 1.5 X 10^^ P erg s"\ (1.21) 

avec P la periode de rotation. Pour les pulsars de type Vela, elle est donnee par : 

oc P"^ erg s-\ (1.22) 

avec P la periode de rotation et Bs le champ magnetique de surface. 



1.4.3 Predictions des Modeles 

Courbes de Lumiere 



La figure 1.8 presente le profil pulse y du pulsar de Vela observe par EGRET (en haut 
a gauche), ainsi que les courbes de lumieres predites par les modeles PC, SO, et GO. Les 
trois modeles arrivent a produire un profil pulse avec deux pics. Cependant, dans le cas du 
modele PC, I'observation d'un tel profil demande un angle d'inclinaison tres petit (< 10°). 
Les deux autres modeles generent des caustiques p°| ( [Morini 1983) en raison de remission y a 



plus haute altitude dans la magnetosphere. L'emission n'a pas la forme d'un faisceau comme 
en radio, mais pent etre vue en forme d"'eventail". Pour le modele SG, on observe l'emission 
de la caustique des 2 poles, tandis que pour le modele GG, seule l'emission d'un des deux 
poles est visible. 

Spectres en Energie 

Le modele Polar Cap predit une coupure rapide (de type super-exponentielle) dans le 
spectre pulse a haute energie, due a la fois a la limite energetique des processus d'emission 
y et a I'attenuation des photons par la creation de paires avec le champ magnetique. Les 
modeles Slot Gap et Guter Gap annoncent, quant a eux, une coupure moins brutale, car leur 
emission est a plus haute altitude. L'energie de ces coupures est estimee entre 1 et 20 GeV et 
depend principalement du champ magnetique a la surface de I'etoile, de la periode du pulsar, 
et de I'altitude d'emission. Le spectre en energie des pulsars pent etre decrit par une loi de 
puissance multipliee par une coupure exponentielle ou super-exponentielle : 

^ = KE-'^e'^/^^--^'', (1.23) 

avec K un facteur de normalisation, V I'indice spectral, Ecutoff la coupure en energie, et |3 
I'indice d'attenuation (:ii 1 pour les modeles GG et SG, et ~ 2 pour le modele PC). Cette 

10. L'effet de caustique est du aux photons emis a des altitudes differentes de I'etoile le long des lignes de 
champ magnetique. Ceci est d'autant plus vrai que I'inclinaison magnetique a est grande. 
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Figure 1.8 - En haut a gauche : Profit pulse du pulsar de Vela, observe par le telescope 
EGRET au-dessus de 30 MeV. Autres figures : Profils pulses prevus par les modeles Polar 
Cap, Slot Gap, et Outer Gap. L'encadre du haut represente la phase en fonction de V angle 
d'observation (Q a partir de I'axe de rotation du pulsar, tandis que l'encadre du has repre- 
sente le nombre de coups en fonction de la phase cf) pour une inclinaison donnee oc. Les 
chiffres 1 et 2 definissent les 2 poles de I'etoile. La figure est extraite de ( Harding 2007b). 



interpretation proposee par Nel and De Jager (1995) est la consequence des differents proces 
sus d'emission cites precedemment. La figure 
dans le cas des modeles PC et OG. 



1.9 presente la simulation du spectre de Vela 



Finalement, I'un des enjeux de la mission Fermi et de cette these est de mesurer la 
coupure de ces spectres pour pouvoir discriminer I'un des modeles, voire de les rejeter. La 
partie Resultats de cette these presente, entre autres, I'analyse spectrale des pulsars detectes 
par le Large Area Telescope. Cependant, la determination des parametres spectraux depend 
d'une bonne calibration de I'appareil, et notamment de la surface efiicace (Chapitre |4]). 
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Figure 1.9 - Spectre d haute energie du 'pulsar de Vela. Les croix sont les donnees du telescope 
EGRET. La ligne pointillee represente les predictions du modele Outer Gap, tandis que la 
ligne discontinue represente les predictions du modele Polar Cap. Les barres d'erreurs sont 



estimees pour un an d' observation. La figure est extraite de (McEnery et al 2004)- 
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2.1 Population Galactique 



2.1.1 Distribution Galactique 



Plus de quarante annees apres la decouverte du premier pulsar radio (Hewish et al 1968), 
la population de pulsars connus a ce jour excede 1800 objets. La majorite provient des 



observations du ciel par les radiotelescopes de Jodrell Bank (Angleterre, Hobbs et al 2004), 



Green Bank Telescope (Eta ts-Unis, [K aplan et al 2005), Nangay (Prance, Theureau et al 
2005 D, Arecibo (Porto Rico, |Dowd et al„ 2000l), et Parkes (Australie, |Manchester„2008D . En 



fait, la moitie des decouvertes resulte de la campagne multi-faisceaux du telescope de Parkes 



durant ces 8 dernieres annees (Manchester et al 2001). La figure 2.1 presente la repartition 
en coordonnees galactiques des 1827 pulsars[J decouverts dans notre Galaxie, et les nuages de 

1. Ce chiffre provient de la base de donnees ATNF (Australia Telescope National Facility). Parmi ces detec- 
tions, 8 pulsars ont ete observes dans les nuages de Magellan, et environ une centaine dans les amas globulaires. 
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Magellan. II faut noter que cette distribution est limitee par la sensibilite des instruments, 
et biaisee par la selection de la region observee. Par exemple, la concentration de pulsars 
vers la longitude 50° provient d'une recherche approfondie du telescope d'Arecibo. D'autre 
part, les pulsars observes sont fortement selectionnes par leur flux apparent F oc L/d^. En 
effet, les pulsars dont la luminosite L est faible, ne sont vus qu'a des petites distances d. 



De plus, les effets du milieu interstellaire sur la propagation du signal radio (voir 2.2.3) 



sont tres importants dans le plan galactique, d'autant plus qu'on se rapproche du centre. La 
consequence de ces nombreux effets est que I'echantillon actuel represente juste une petite 
part de la population totale des pulsars actifs dans la galaxie. Celle-ci est estimee entre 10^ 



et 10^ (Lyne and Graham-Smith 1998) 



180. 




Longitude Gala-clique 



Figure 2.1 - Distribution des pulsars connus dans le del en coordonnees galactiques, en 
projection de Hammer- Aitoff (source :ATNF). 



La figure 2A_ indique que les pulsars se repartissent essentiellement dans le plan galactique 
comme les etoiles massives O et B, sur une epaisseur d'environ 1 kiloparsec (kpc) pour un 
rayon de 10 kpc autour du centre galactique. Cette observation conforte I'hypothese que les 
etoiles a neutrons naissent de I'effondrement gravitationnel des etoiles massives. |Kaspi| ([l998) 
associe une vingtaine de pulsars radio a des restes de supernova. II apparait egalement que 
de nombreux pulsars sont en dehors du plan (| b | grand). Ceci pent etre du a la proximite du 
pulsar par rapport a la Terre, ou il est possible que les explosions violentes de supernova ne 



soient pas totalement symetriques, donnant au pulsar une impulsion a leur naissance (Janka 



et al 2005). La mesure des mouvements propres d'un echantillon de 233 pulsars montre une 
Vitesse moyenne de propagation de 246 ± 22kms^^ (Hobbs et al 2006b). Les plus rapides 
peuvent atteindre 1000 km s^^ La consequence est qu'une fraction importante des pulsars 
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s'echappe de leur reste de supernova, et subsequemment du plan galactique. 
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Figure 2.2 - Diagramme P — P, distribution des pulsars en fonction de leur periode de 
rotation et de leur derivee. Les lignes vertes representent I'dge caracteristique . Les lignes 
bleues representent le champ magnetique a la surface de I'etoile, tandis que les lignes rouges 
representent le taux de perte d'energie de rotation E. Les points noirs sont les pulsars listes 
dans le catalogue ATNF. Les etoiles bleues sont les pulsars detectes par la mission EGRET 
avec un bon niveau de confiance, tandis que les triangles rouges sont ceux detectes avec 
un faible niveau de confiance. Les principales caracteristiques des 9 pulsars EGRET sont 
enumerees dans le tableau 
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2.1.2 Diagramme P-P 

Dans le chapitre [ij nous avons vu que seules les mesures de la periode de rotation P 
et du taux de freinage P pouvaient nous renseigner sur revolution de la rotation et les 

appelee le diagramme P — P, montre la 



2.2 



proprietes intrinseques des pulsars. La figure 
distribution des pulsars en fonction de la periode de rotation et de sa derivee. Deux groupes 
sont clairement distin cts. O n y retrouve en haut a droite les pulsars normaux (P ~ 0.5 s et 



P ~ 10^ ss^ , section 2.1.3), tandis que le petit g roup e en bas a gauche represente les pulsars 
milliseconde (P ~ 3 ms et P ~ 10^^°ss^\ section 2.1.4). La difference dans P et P implique 
fondamentalement des differences d'age et de champ magnetique pour ces deux populations. 
Le champ magnetique B oc (PP)°-^, I'age caracteristique t oc P/P sont indiques par differentes 
lignes. Le taux de perte d'energie de rotation E oc P/P^ est egalement indique. Les pulsars 
naissent avec une periode relativement petite et subissent un ralentissement. L'evolution 
naturelle pour un pulsar normal est ainsi de migrer des petites periodes de rotation (region 
en haut a gauche) vers les plus grandes periodes dans un intervalle de temps de 10^ ^ ^ annees. 
Au dela de 10^ annees, le flux du pulsar est certainement trop faible pour etre detectable. 
Remarquons que le pulsar du Crabe et de Vela font partie de la classe de pulsars la plus 
jeune, la plus energetique, avec un champ magnetique intense. 

Les etoiles bleues representent les pulsars y detectes par la mission EGRET, tandis que les 



triangles rouges representent les detections marginales (1.3.1). Ces pulsars sont relativement 



jeunes et sont concentres (a I'exception du pulsar milliseconde candidat) dans la region ou le 
champ magnetique est intense. Tous ont egalement un taux de perte d'energie de rotation E 
eleve. Pour le Crabe, il est de 4.5 x 1 0^^ erg s^^ , tandis que pour Vela, il est de 7 x 1 0^^ erg s^^ . 
On constate ainsi que les pulsars observes a haute energie ont un profil P-P similaire, ce qui 
va permettre de selectionner nos pulsars candidats y. 



2.1.3 Pulsars Normaux 

La premiere classe definit les pulsars "normaux" que Ton observent apres la phase de 
supernova, tels que PSR B1919-I-21 (le premier pulsar detecte), et les pulsars du Crabe et de 
Vela. lis ralentissent a cause du rayonnement dipolaire et s'eteignent au bout d'une dizaine 
de millions d'annees (voir 1.2.3). lis sont concentres le long du plan galactique et represente 
90% des pulsars observes. Signalons que dans cette these, nous analyserons essentiellement 
des pulsars normaux (chapitres 6 a 9). 



2.1.4 Pulsars Milliseconde 

La seconde classe identifie les pulsars "milliseconde". 80% d'entre eux sont presents dans 
des systemes binaires (moins de 1% pour les pulsars normaux). Initialement, le systeme est 
compose de deux etoiles de la sequence principale. L'etoile la plus massive va evoluer plus vite 
et s'effondrer pour former une etoile a neutrons. Dans moins de 10% des cas, malgre la violence 
de I'explosion de la supernova, l'etoile reste liee gravitationnellement a son compagnon. Si les 
parametres orbitaux sont favorables, le champ gravitationnel intense de l'etoile a neutrons 
attire de la matiere de l'etoile compagnon et forme un disque d'accretion. Le resultat de ce 
transfert de masse et de moment angulaire, est un rayonnement X intense du systeme, ainsi 
qu'une acceleration de la rotation de l'etoile a neutrons jusqu'a des periodes de I'ordre de la 
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milliseconde (P < 30ms). Les particules chargees dans renvironnement du nouveau pulsar 
vont etre accelerees et produire un faisceau d'onde electromagnetique. 

L'age caracteristique des pulsars milliseconde est compris entre 10^ et 10^° annees. Leur 
taux de freinage est relativement petit (P < lO^^^ss^^) compare a celui des pulsars normaux. 
Quant au champ magnetique de surface, il est possible qu'il se dissipe lors de la phase d'ac- 
cretion (Bisnovatyi-Kogan and Romberg 1974). II est 10"^ plus petit que pour les pulsars 



normaux. Finalement, leur vitesse de propagation est de 87 ± 13kms en moyenne. 

Le premier pulsar milliseconde a ete decouvert en 1982 a Arecibo avec une periode de 
1.6ms (Backer et al, 1982). A ce jour, 168 pulsars sont listes dans le catalogue ATNF avec 
une periode inferieure a 30ms et un freinage P < 10^^''. Parmi eux, 96 se trouvent dans les 
amas globulaires ( |Camilo and Rasio 2005), et 8 dans les nuages de Magellan. Pour plus de 
details, on se reportera a la these intitulee "Detections de pulsars milliseconde avec le FERMI 
Large Area Telescope" (|Guillemot| |2009D. 



2.2 Distance des pulsars 



Dans cette section, nous presentons les differentes methodes utilisees pour estimer la dis- 



tance des pulsars. Cette quantite difficile a mesurer de fagon directe (2.2.1) est essentielle a 
la fois pour les etudes statistiques des regions de formation des etoiles massives, mais ega- 
lement pour comprendre les proprietes physiques des pulsars, notamment la luminosite. A 
noter qu'aucune relation ne permet encore d'evaluer la distance a partir des parametres des 
pulsars tels que leur luminosite y. Cependant, la situation pourrait changer, etant donne le 
nombre grandissant de pulsars y detectes avec la mission Fermi. La distance est egalement 
primordiale pour selectionner nos candidats y. Rappelons que le flux observe est une fonction 
de la luminosite du pulsar, mais aussi de la distance au carre (F oc L/d^). 



II existe plusieurs methodes pour evaluer la distance des pulsars, soit par la mesure de leur 



parallaxe (2.2.1), soit par I'absorption d'hydrogene neutre et I'association d'objet (2.2.2), ou 



encore par la mesure de la dispersion du signal radio couple a un modele de la distribution 



de densite des electrons ; les deux les plus utilises sont le modele de Taylor and Cordes (1993) 



et le NE 2001 (Cordes and Lazio, 2002) (2.2.3). Les 2 premieres methodes pour estimer 



la distance permettent de contraindre les modeles de densite des electrons. Le tableau [271 
presente la distance des pulsars y du premier catalogue de pulsars du Large Area Telescope 
(Chapitrelgl). 



2.2.1 Parallaxe 



Les distances des pulsars les plus proches, situes a moins de 1 kpc du Soleil, peuvent etre 
obtenues a partir de la mesure de leur parallaxe. La methode consiste a mesurer la courbure 
des fronts d'onde radio provenant du pulsar au long d'une annee complete. La courbure est 
evaluee par la mesure tres precise du temps d'arrivee des pulsations a differents endroits de 
I'orbite terrestre. 
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2.2.2 Absorption d'Hydrogene et Association d'Objet 



Le spectre des pulsars situes a des latitudes galactiques faibles, a I'interieur des bras 
spiraux de la Galaxie, peuvent presenter une rale d'absorption de I'hydrogene neutre. L'hy- 
drogene present sous forme de nuage absorbe le rayonnement a la longueur d'onde de 21 cm 
(1420 MHz). Cette raie Hi peut montrer un decalage spectral du a I'effet Doppler cause par le 
deplacement du gaz dans les bras spiraux. A partir des modeles dynamiques de la rotation de 



la Galaxie, le decalage spectral de I'liydrogene peut etre relie a une distance (Fich et al 1989). 
Cette technique peut etre egalement utilisee pour estimer la distance d'un objet associe au 
pulsar (ex : restes de supernova, nuages de Magellan, amas globulaires), par I'absorption ou 
remission de I'hydrogene neutre. Dans le cas du pulsar PSR J0205+6449, la distance est 
determinee a partir du PWN decouvert avant le pulsar ( [Roberts et al 1993). 

Les mesures d'absorption de I'hydrogene neutre (Hi) a 21 cm sont une importante estima- 
tion des diverses proprietes du milieu interstellaire, et permettent de contraindre les modeles 
de densite des electrons. 



2.2.3 Mesure de Dispersion 

Avant d'arriver sur Terre, la large bande d'emission radio emise par le pulsar traverse 
un plasma compose principalement d'electrons. Celui-ci provient du milieu interstellaire. La 
consequence est que les ondes radio de basse frequence y-i arrivent avec un retard par rapport 
aux ondes radio de haute frequence yz- Celle-ci peut en principe etre corrigee avec exactitude. 
Le retard en temps relatif entre deux signaux pulses pour deux frequences differentes (vi et 
V2) peut etre enonce par : 

lrl-l)DM, (2.1) 



At 



ou K 
et 



2.410331 X 10 ^MHz ^cm ^s ' pc est la constante de dispersion (Backer et al 1993) 



DM 



TLedl, 



(2.2) 



est la mesure de dispersion (pc cm^^) avec L la distance entre la Terre et le pulsar en parsec 
et Tie la densite d'electrons en cm^^. L'acronyme DM provient du terme anglais "Dispersion 
Measure". Generalement, la dispersion est determinee durant la decouverte du pulsar, a par- 



tir du temps d'arrivee des photons et de deux frequences differentes (voir I'equation 2.1). La 



figure 2.3 montre le signal regu pour chaque canal de frequence en fonction de la phase du 



pulsar pour PSR J0437— 4714. Signalons que pour les photons X ou y, ce plasma n'a aucun 
effet. 



La mesure de dispersion est un excellent outil pour non seulement identifier remission 
radio des pulsars, mais aussi pour modeliser la structure du gaz ionise dans le milieu in- 



terstellaire. Taylor and Cordes (1993) et Cordes and Lazio (2002) ont construit des modeles 
quantitatifs sur la distribution des electrons libres dans la galaxie. Ces modeles sont contraints 
a partir des mesures de dispersion des pulsars et de leurs distances obtenues par des mesures 



independantes (voir 2.2.1 et 2.2.2). lis sont modelises suivant la forme des bras spiraux, 
qui derivent d'observations radio et optiques des regions Hn . Generalement, on suppose des 
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Figure 2.3 - Gauche : Signal radio regu pour chaque canal de frequence en fonction de la 
phase du pulsar pour PSR 304.37—4.714.. En V absence de correction de dispersion du signal, 
les basses frequences radio arrivent apres les hautes frequences. Droite : Profil radio apres 
V application de la mesure de dispersion. Les donnees proviennent du telescope radio de 
Parkes en Australie. La figure est extraite de Han (2009) 



densites uniformes d'electrons dans les bras spiraux et des transitions entre les bras. La conse- 
quence est qu'a partir d'un modele de densite d'electrons et de la mesure de la dispersion, on 
peut estimer la distance du pulsarj^ Cependant, I'incertitude de ces distances derivees pent 
exceder 50% pour certains pulsars. N'oublions pas que le milieu interstellaire est irregulier. 
Le signal radio peut done etre deteriore en fonction du temps et de la frequence observee. 



2.3 Les Candidats Pour Fermi 



2.3.1 La Campagne pour les Ephemerides 



Parmi les 1827 pulsars connus ( Manchester et all 2005), seulement sept pulsars ont ete de 



tectes comme emetteurs y par la mission CGRO. Avant le lancement du Large Area Telescope 
(Chapitre [3]), les previsions sur les decouvertes de pulsars y par le telescope tablaient entre 
une dizaine et plusieurs centaines de detections selon les modeles CHarding et al 2007 Jiang 



and Zhang, |2006) . Les deux methodes pour la detection de pulsations y sont : d'une part 



la recherche de pulsations en aveugle en utilisant uniquement les donnees du LAT (Atwood 



et al 2006 Abdo et al 2009-BS), et d'autre part la recherche de pulsations pour des pulsars 



deja connus dans une autre longueur d'onde. 

Nous avons fait le choix a Bordeaux d'etudier les pulsars prealablement observes a une 
autre frequence, le plus souvent en radio. La collaboration Fermi a ainsi demarre, a partir de 
2007, une campagne de suivi des pulsars, avec la cooperation des radiotelescopes enumeres a 



2. http ://rsd-www.nrl.navy.iiiil/7213/lazio/ne_niodel/ 
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Nom du Pulsar 


Distance (kpc) 


Methode 


J0007+7303 


1.4it0.3 


HI 


JOOciO+0451 


0.300±0.09U 


P 


<J VJ ^ U >J U 1 ± 


2.6—3.2 


HI 


J0218+4232 


2.5-4 


A 


J0248+6021 


2-9 


A 


J0437-4715 


0.1563±0.0013 


P 


J0534+2200 


2.0±0.5 


A 


J0613-0200 


0.48lg]^ 


P 


J0631+1036 


0.75-3.62 


A 


J0633+1746 


0.250+O-J^o 


P 


J0659+1414 


0.288+00^3 


P 


J0742— 2822 


2 07+] -J? 
-1.07 


DM 


Tn7m -1-1 R07 


^■°-0.2 


p 


J0835— 4510 
J1028— 5819 


n oQ'7+0-019 

0-287_o.oi7 
2.33±0.70 


P 
DM 


J1048— 5832 


2.71it0.81 


DM 


J1057— 5226 


0.72 ±0.2 


DM 


J1124— 5916 


. n + O 7 

4.8_,2 


A 


J1418-6058 


2-5' 


A 


J1420— 6048 


5.6±1.7 


DM 


J1509— 5850 


2.6±0.8 


DM 


J1614— 2230 


1.27it0.39 


DM 


J1709— 4429 


1.4—3.6 


A 


T1'71D QQOC 


O QO_L 1 1 C 

o.cszihl.lb 


JJM 




n 9s-t-n 1 1 
u.oomu. 1 1 


UiVi 


T1 '7 A A 1 1 Qyl 

J1747— 2958 


n ciC7+0.043 
2.0it0.6 


r 

DM 


J loUy — Z66z 


1.7ibl.O 


rli 


J1833-1034 


4.7±0.4 


HI 


J1836-I-5925 


<0.8 


A 


J1952-I-3252 


2.0±0.5 


HI 


J2021-I-3651 


2 1+21 


A 


J2021-I-4026 


1.5±0.45 


HI 


J2032-I-4127 


3.6±1.08 


DM 


J2043-I-2740 


1.80±0.54 


DM 


J2124-3358 
J2229-I-6114 


2B+°-^^ 
""^^-0.08 

0.8-6.5 


P 
A 



Tableau 2.1 - Estimation de la distance des pulsars du catalogue (chapitre^. 
HI : distance evaluee par I'observation d'HI ; P : mesure de la parallaxe ; DM 

\20M); A 



dispersion ( Cordes and Lazio 



mesure de 

autres mesures. La reference des distances pour 



chaque pulsar est indiquee dans le papier Abdo et al (2009-psrcat) 



la section 2.1.1 et des telescopes XMM et RXTE ( Smith et al[ 20081). Les communautes radio 
et X se sont engagees a fournir les ephemerides des pulsars (voir la section 5.2.2), c'est-a-dire, 
les parametres de rotation qui permettent en controlant les variations de la rotation du pulsar, 
d'integrer a la bonne phase rotationnelle les photons y qui nous parviennent. Seulement, le 
LAT regoit du pulsar de Vela, qui est la source la plus emissive du ciel y, en moyenne un 
photon tous les 4 minutes. Cela est tres peu compare aux autres longueurs d'onde. Pour les 
pulsars les plus faibles, le temps entre deux evenements pent etre de quelques jours. II est done 
indispensable pour une detection efficace des pulsars y et une analyse rigoureuse, d'avoir des 
ephemerides de qualite. Les pulsars montrent egalement des irregularites dans leur rotation 
{timing noise et glitch). La consequence est qu'il est tres difficile d'extrapoler la periode de 
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rotation et ses derivees sur une tres grande periode. Un exemple concret est la decouverte 
des pulsars avec un fort t, PSR J2229+6114 et PSR J1410-6132, apres la fin de la mission 
CGRO. Une recherche de pulsations dans les archives des donnees EGRET a ete effectuee, 
mais a echoue pour trouver un signal pulse car les photons etaient trop peu nombreux et trop 



vieux compares aux ephemerides (Thompson et al 2002 O'Brien et al 2008). L' observation 
des pulsars par les telescopes radio et X doit etre reguliere, mais represente un cout, ce qui 
justifie de determiner une liste bien precise de pulsars candidats. La section suivante presente 
la selection de nos pulsars candidats y. 



2.3.2 La Selection 



Arons (1996) montre qu'au-dessus d'une certaine valeur du potentiel au niveau des lignes 
de champs ouvertes V, les cascades de paires qui generent les photons y dans la magne- 
tosphere vont se produire. Ce potentiel proportionnel a E pent etre decrit en fonction de la 
periode d e rota tion du pulsar et de sa derivee par V = 4 x 10^° p-3/2 pi/2 ~ 3,18 x 10^^\/e 
(equation 1.19 en fonction de P et P). De plus, I'efiicacite r\y = L-^/E augmente quand ce 

potentiel au niveau des lignes de champs ouvertes decroit, soit riy ~ 1 /V ~ 1 / La conse- 
quence est qu'il existe une limite de production du rayonnement y. Basee sur les 7 pulsars 
detectes a haute energie, cette limite semble etre proche de E 3 x 1 0^'* erg s^^ . Suivant 
ces observations empiriques, une liste de 224 pulsars candidats y a ete selectionnee avec 

E > 1 X lO^'^ergs^^ Les sources ont ete classees selon la quantite \/E/d^, qui correspond au 
flux du pulsar incident sur Terre, et souligne I'importance sur I'estimation de la distance d. 
Ces criteres de selection ainsi que la liste complete des sources sont references dans le papier 



Smith et al (2008) 



Le tableau 2.2 presente parmi la liste des 224 objets, une sous-selection des pulsars nor- 
maux associes a une source EGRET non identifiee (^Hartm an et al[ 1999), ou/et a un reste 
de supernova / nebuleuse a vent de pulsars (Kaspi ,1998) , ou/et a une source detectee au 
TeV par I'experience Cherenkov HESSj^ Le but de ce tableau est de pre-definir des objets 
potentiellement interessants pour la suite de la mission (et de la these), et de signaler que I'ob- 
servation y n'est pas une finalite, mais que la comprehension de ces sources passe egalement 
par les autres longueurs d'onde et les sources associees aux pulsars. Signalons au lecteur que 
ce tableau est base sur des observations ou des references de I'annee 2007 (ou anterieures). 

Les pulsars sont ranges selon \/E/d^, ou d est la distance obtenue dans le catalogue ATNF 
("Dist"). Ce tri ne prend pas en compte la geometrie du faisceau par rapport a la ligne de 
visee terrestre. De plus, la colonne "EGRET events" presente les photons y detectes par 
EGRET au-dessus de 10 GeV dans un rayon de 1° autour de la position du pulsar, sur toute 
la duree de la mission. Sur les 1506 photons detectes au-dessus de 10 GeV sur tout le ciel, 35 
sont localises a moins de 1 degre de la position des 6 pulsars EGRET. Cette etude est basee 
sur les travaux de (jThompson et al 2005). Le tableau resume egalement remission pulsee 
detectee dans les differentes longueurs d'onde, de la radio au rayonnement y (O : optique, 
Xs : X mou, X^ : X dur, Gj : gamma mou, Gh : gamma dur). Les deux dernieres colonnes 
presentent respectivement la densite du flux moyen radio pour une observation a 1400 MHz, 
et I'initiale du radiotelescope observant le pulsar (P : Parkes, N : Nangay, J : Jodrell). 

3. High Energy Stereoscopic System : http ://www.nipi-hd.inpg.de/hfni/HESS/ 
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Precisons finalement que la collaboration avec les astronomes des radiotelescopes est un 
veritable succes. Les groupes bases a Manchester (Jodrell Bank) et a Nangay nous ont fourni 
en 2008, suite a nos nombreuses decouvertes, plus de 500 ephemerides des pulsars qu'ils 
observent regulierement. La plupart d'entre eux presente un ralentissement E inferieur a 
10^^ergs^\ Cela nous permet de couvrir un plus grand espace de phase P-P et de ne pas 
imposer une limite theorique sur E pour la production de photons y. Les resultats sont 
resumes au chapitre [9| (catalogue de pulsars). 
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3.1 Le Telescope Spatial Fermi 



Le telescope spatial Fermi Gamma-ray Space Telescope (figure 3.1) est une mission 
scientifique Internationale dediee a I'observation du ciel des hautes energies. II est compose 
de deux instruments realises par la collaboration de divers instituts, agences spatiales et uni- 
versites d'AUemagne, des Etats-Unis, de France, d'ltalie, du Japon et de Suede. Cinq equipes 
fran^aises de riN2P3/CNRS, de I'lNSU/CNRS et de I'lRFU/CEA contribuent a ce projet. 
Le satellite de 4,2 tonnes a ete lance le 11 juin 2008 par la NASA a Cap Canaveral en Floride 
par une fusee Delta II Boeing. II suit actuellement une orbite quasi circulaire autour de la 
Terre, de periode 95 minutes, a une altitude d'environ 550 km. La duree de la mission est 
de 5 ans, prolongeable jusqu'a 10 ans. Avant le decoUage et la phase de verification, le te- 
lescope portait le nom de GLAST pour Gamma-ray Large Area Space Telescope. La NASA 
I'a rebaptise Fermi, en I'honneur du celebre scientifique italien Enrico Fermi, pionnier des 
mecanismes d'acceleration des rayons cosmiques. 

L'objectif de la mission Fermi est d'etudier I'univers a travers le rayonnement y de haute 
energie. Cela inclut : (1) la comprehension des mecanismes d'acceleration des particules dans 
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les noyaux actifs de galaxie, les pulsars, les nebuleuses a vent de pulsar, et les restes de super- 
nova, (2) I'identification des sources non-identifiees et I'origine de remission diffuse relevee 
par EGRET, (3) de comprendre le comportement des sursauts y , des sources transitoires, et 
des sursauts solaires, et finalement (4) d'utiliser I'observation du rayonnement y comme une 
sonde de la matiere noire. 

Sa construction est justifiee par le fait qu'aucun detecteur au sol n'observe de photons 
dans sa gamme d'energie. Les detecteurs Cherenkov comme HESS et MAGIC, ont un seuil 
en energie superieur a environ 30 GeV malgre leur grande surface collective. lis observent les 
photons crees par effet Cherenkov, issus des gerbes electromagnetiques produites dans I'at- 
mosphere. Cependant, en dessous de cette limite, les gerbes prennent place trop haut dans 
I'atmosphere pour que les telescopes actuels les observent, et le bruit de fond des particules 
cosmiques chargees noit le signal. L'instrument permet egalement de realiser des campagnes 
multi-longueur d'onde avec d'autres telescopes. 




Figure 3.1 - Illustration du satellite Fermi Gamma-ray Space Telescope. 



Le telescope est compose de deux instruments (figure [372 ) : 



O le Large Area Telescope (LAT), l'instrument principal du telescope (Atwood et al 2009). 
II mesure la direction, I'energie entre 0.02 et 300 GeV[^ et le temps d'arrivee des photons 
ainsi que des particules chargees. II a la particularite d'explorer I'ensemble du ciel en 
trois heures grace a son tres grand champ de vue (20% du ciel a tout moment). De 
nombreuses sources de rayons gamma etant variables, cette surveillance continuelle du 
ciel permet d'alerter la communaute scientifique en cas d'eruptions. Signalons que pour 
I'etude des pulsars qui est le sujet principal de cette these, seules les donnees et les 
performances provenant du LAT seront utilisees. 

@ le Gamma-ray Burst Monitor (GBM) qui vient completer le LAT, est dedie a I'etude 
des sursauts gamma entre lOkeV et 30 MeV. II est le successeur de l'instrument BATSE 
sur CGRO. II est compose de 12 detecteurs scintillateurs Nal et de 2 detecteurs scin- 
tillateurs BGO, ceci avec un angle solide de 8 sr pour une precision de 10 arcmin. II sert 
notamment de declencheur pour le repointage du LAT en cas de detection d'un sursaut 



y. Pour plus d'informations, on pourra se reporter a I'article de Meegan et al (2009). 



1. Au-dela de 300 GeV, la resolution en energie est superieure a 20%, et la reconstruction des evenements 
de vient tres difficile. 
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Notons que la contrainte due a renvironnement spatial impose un cout et une conception 
bien particuliere, lesquels jouent sur les performances de I'instrument telles que la surface du 
detecteur et done la sensibilite de I'appareil. La puissance electrique est egalement limitee 
a 700 W, ce qui contraint le traitement des donnees a bord du satellite et les technologies 
(par exemple, la limitation des photomultiplicateurs). Finalement, le satellite doit pouvoir 
resister aux conditions spatiales. Une serie de tests tels que la resistance aux vibrations, et 
la resistance aux basses et hautes temperatures, a ete effectuee avant le decoUage. 




Figure 3.2 - Le Large Area Telescope (LAT) et le Gamma-ray Burst Monitor (GBM) 



3.2 Principes de Detection des y par le LAT 

Interessons-nous a present aux principes d'interactions du LAT. Nous nous plagons dans 
le cadre de I'electromagnetisme. Les processus d'interactions y-matiere varient en fonction 
de I'energie de I'onde electromagnetique/photon. Aux energies du LAT (au-dessus du MeV), 
les principales interactions sont I'ionisation, le bremsstrahlung ou rayonnement de freinage, 
et la creation de paires e^. 



Au-dessus du MeV, le photon y qui entre dans le trajectographe (section |3.3.1 ), a une 
chance d'interagir avec le champ electromagnetique des atomes du milieu, et de creer une paire 
electron-positon e^. Ce mecanisme est illustre par la figure 3.2 La probabilite d'interaction 



(ou section efRcace) d'un photon depend de la densite du materiau qu'il traverse. La paire e 



± 
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creee va ensuite interagir avec le milieu. Selon les equations de Maxwell, toute charge dont la 
Vitesse varie en valeur absolue ou en direction rayonne d'autant plus facilement que sa masse 
est petite. Ce rayonnement s'appelle le rayonnement de freinage, decrit par + A — > e^ + y, 
avec A le champ electromagnetique du milieu. Les electrons et positons vont ainsi emettre 
des y qui vont a leur tour creer des paires e^, et ainsi de suite. Les premieres populations 
de photons et de paires vont se derouler dans le trajectographe, tandis que le reste de la 
cascade electromagnetique va se former dans le calorimetre (section 3.3.2), a cause de sa plus 
grande densite. A chaque iteration I'energie de depart est divisee, et les photons y sortent peu 
a peu du regime de creation de paires. Ensuite, chaque electron de la gerbe depose de I'energie 
par ionisation et rayonnement de freinage dans le Csl du calorimetre. Le depot d'energie par 
ionisation est decrit par la formule de Bethe-Bloch. 

Notons que les sections efRcaces de creation de paires et de rayonnement de freinage 
peuvent etre decrites a I'aide de la longueur de radiation Xq. Celle-ci depend du materiau 
traverse, soit Xq ~ 716gcm^^A/Z^. La consequence est que la resolution en energie qui derive 
de la reconstruction de la gerbe electromagnetique, depend de la profondeur du calorimetre 
du LAT. Au dela d'une energie initiale de 30 GeV, le calorimetre mesure moins de la moitie de 
la cascade. II est done tres difficile ensuite de reconstruire I'energie de la particule incidente. 
Le deuxieme facteur limitant dans la reconstruction des particules est la diffusion multiple. 
Quand un electron traverse un materiau, il est devie par le champ des atomes, d'autant plus 
que sa vitesse est petite. Si I'electron a une trop grande deviation, le trajectographe ne pourra 
pas identifier correctement la direction initiale de la particule. 

Ces processus physiques vont limiter les performances de I'instrument, et ainsi contraindre 
la surface efficace, la resolution en energie, et la resolution angulaire. Le chapitre suivant decrit 
ces caracteristiques et presente la validation de la surface efficace. 



3.3 Fonctionnement de I'instrument 

Le Large Area Telescope (LAT) est compose de 16 tours formant un carre 4x4 enveloppe 
par un bouclier d'anti-coi'ncidence pour 1' identification des particules chargees provenant du 
rayonnement cosmique et de I'albedo terrestre. Chaque tour est elle-meme composee d'un 
trajectographe, d'un calorimetre, et d'un systeme d'acquisition. Ces differentes parties seront 



decrites par la suite. Le tableau 3.1 enumere les principaux parametres et capacites du LAT. 



3.3.1 Trajectographe 

Le trajectographe est I'un des sous-detecteurs utilise dans la determination de la direction 
des particules incidentes. Cette direction est reconstruite a partir d'algorithmes utilisant les 
coordonnees {xj^^, yj^^, zj^^) des pistes silicium i activees par les paires creees par le 
photon incident. Le trajectographe est egalement la voie principale de declenchement de 
I'acquisition du LAT. 

Un module du systeme est compose de 18 plans x - y, sensibles au passage d'une particule 



chargee (figure 3.3). Chaque plan x et y est compose de 16 detecteurs a piste de silicium 
(Silicon Strip Detector, SSD), qui permettent de visualiser les traces laissees par les particules 
chargees. La dimension z est obtenue en empilant les plans. La taille d'un detecteur est 8.95 
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Parametres 


Valeur ou intervalle 


Intervalle en energie 


0.02 - 300 GeV 


Surface efiicace a incidence normale 




0.1 - lOGeV 


7000 cm^ 


10 - 300 GeV 


8000 cm^ 


Resolution en energie (la gaussien) : 




0.1 - 10 GeV sur I'axe 


< 10% 


10 - 300 GeV sur I'axe 


< 20% 


> 10 GeV (a >60 ° d'incidence) 


< 6% 


Resolution angulaire (rayon de confinement de 68%) : 




> 10 GeV 


< 0.15° 


IGeV 


0.5° 


100 MeV 


3.5° 


Champ de vue 


2.4 sr 


Precision temporelle 


< 1 (xsec 


Temps mort 


26.5 usee 



Tableau 3.1 - Parametres et capacites du LAT pour les fonctions de reponse de I'instrument 
(IRFs) appelees Pass6_v3. 



X 8.95 cm pour une epaisseur de 400 |j.m, et compte 384 pistes de silicium (jonction PN). 
Leur efiicacite de declenchement est superieure a 99%, et le bruit electronique est de 45±8 
pistes declenchees lors du passage d'un evenement, ceci pour I'ensemble des tours. Notons 
qu'un evenement allume quelques dizaines de pistes concentrees sur un cylindre de quelques 
centimetres de rayon. Cette technologie apporte egalement un temps morts beaucoup plus 
court, par rapport a la chambre a ionisation utilisee par le detecteur EGRET (~ 5 ms). 

Entre ces plans sont intercales 16 plans convertisseurs en tungstene permettant de conver- 
tir le photon incident en une paire e^. Etant donne que la probabilite d'interaction du photon 
depend de la densite du materiau traverse, le trajectographe a ete congu pour amorcer la 
cascade electromagnetique qui se developpera pour le reste dans le calorimetre. Ces plans 
convertisseurs sont maintenus par une structure en carbone composite, et ont des epaisseurs 
differentes. Les 12 premiers plans (section mince) ont une epaisseur de 0.03 longueur de radia- 
tion (Xo), pour eviter que les electrons et les positons n'interagissent trop par rayonnement 
de freinage ou par diffusion multiple. Cela permet d'augmenter la qualite des traces, mais 
diminue la probabilite de conversion en paire pour le photon y entrant. Les 4 plans sui- 
vants (section epaisse) ont une epaisseur de 0.18 longueur de radiation, permettant un taux 
de conversion des photons de 65%. Cette caracteristique donne une plus grande sensibilite 
aux sources faibles, mais degrade consider ablement la reconstruction des photons de basse 
energie (<100MeV). Pour ces photons, la gerbe electromagnetique ne se developpe pas jus- 
qu'au calorimetre. Ainsi, les 2 derniers plans ne possedent pas de plans convertisseurs pour 
eviter d'accentuer la conversion. 

La reconstruction du vertex cree par la conversion en paire du photon incident se fait 
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GLAST complet 

Reseau de 4x4 modules 



Photon gamma 




Rejet des rayons cosmiques charges 

systeme anti-coi'ncidences 
Tuiles de scintillateur 




Figure 3.3 - Module trajectographique du LAT. 



au sol, et est determinee par des algorithmes de reconnaissance de formes, tels que le track 
finding et le track fitting. Les traces sont reconstruites en tenant compte du nombre de 
pistes touchees i, et de la diffusion multiple qui donne une indication sur I'energie des paires 



e . Notons qu'un autre algorithme utilise le procede appele filtre de Kalman (Jones 1998), 



dont le principe est de propager un vecteur d'etat d'une couche vers la couche inferieure pour 
reconstruire la trace. Finalement, une seconde analyse basee sur des "Classification Trees" 



(voir la section 4.3) selectionne le meilleur algorithme pour reconstruire la direction de la 
particule incidente notee 



Les principales caracteristiques du trajectographe sont detaillees dans le tableau 3.2 



amsi 



que dans |Atwood et al ( 2007| ). Ce sous-syteme joue un role principal dans la resolution 
angulaire du LAT. Cette performance sera ainsi limitee par la trace, et done par I'imprecision 
sur la direction du photon incident. 



3.3.2 Calorimetre 

Le calorimetre du LAT a pour but principal de mesurer I'energie Ey des photons inci- 
dents au-dessus de 20 MeV. II est egalement utilise pour la trajectographie des evenements et 
I'identification des rayons cosmiques. La figure [3741 presente un des modules du calorimetre et 
expose sa structure en couche des faces laterales : plaque de fermeture en aluminium, circuit 

2. est defini comme Tangle d 'incidence, soit Tangle entre Taxe normal au plan du LAT et la direction du 
photon incident. 
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Caracteristiques 


Performance 


Epaisseur d une SSD 


0.4 mm 


Nombre de pistes par SSD 


384 


"!~\ "J. ' J. 1 "J. 

Distance separant les pistes 


228 i^m 


Efiicacite de declenclienient 


>yy /o 


Pistes declenchees par le bruit electronique 


10-^ 


Rapport hauteur /largeur 


0.4 


Xo des 12 premiers plans convertisseurs (front) 


12 X 0.03 


Xo des 4 derniers plans convertisseurs (back) 


4 X 0.18 


Tableau 3.2 - Caracteristiques du trajectographe 


du LAT. 



electronique de lecture, puis protection thermique et electromagnetique en aluminium. Sa 
masse totale, qui est limitee par la necessite d'une mise en orbite du satellite, est de 1800 kg. 

Le principe du calorimetre est de creer une cascade electromagnetique grace a sa petite 
longueur de radiation (Xo(CsI) = 1.85cm), et de mesurer I'energie (ef'^'") et la position 
(xf^^jXj'l^^jZf^^) des depots d'energie. Les photons et les ionisent les cristaux, qui emettent 
a la suite de la lumiere visible creee par la desexcitation de leurs atomes. La lumiere est ensuite 
recuperee par des photodiodes et convertie en signal analogique. 



Structure en carbone 



CDE : Cristal de Csl equipe 

Fermelure de la stucture 

Electronique 
de lecture 




Plaque de support 



Protection de I 'electronique 
Boucher electromagnetique 



Figure 3.4 - Module du calorimetre du LAT (Bregeon 2005). 



Chacun des 16 modules du calorimetre contient 8 couches alternativement croisees de 12 
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barreaux de lodure de Cesium dope au thalmm CsI(Tl), lesquelles sont soutenues par une 
structure en carbone composite, developpee et construite au laboratoire Leprince-Ringuet 
de I'Ecole Poly technique. Chaque couche dont I'epaisseur est de 19.9 mm, soit 1.08 Xq, est 
orthogonale a la precedente, formant une structure hodoscopique. Les barreaux qui mesurent 
326 mm x 26.7 mm x 19.9 mm, ont a leurs extremites une double photodiode, laquelle est 
composee d'une grande diode (LE : Low energy), et d'une diode six fois plus petite (HE : 
High Energy). De plus, la sortie de chaque diode est connectee a deux circuits electroniques 
permettant d'amplifier le signal. Ce systeme permet ainsi une lecture du depot d'energie sur 
quatre gammes, soit un intervalle s'etendant de 2 MeV a 102.4 GeV, avec un recouvrement 



des gammes deux a deux (voir le schema 3.5) 



Le systeme permet egalement d'estimer la position moyenne du depot d'energie des pho- 
tons et des e^. La mesure se fait grace a I'attenuation de la lumiere lors de sa propagation a 
I'interieur d'un barreau, dont la surface est depolie. Le signal est lu aux deux extremites du 
barreau. Son amplitude est asymetrique en fonction de la position longitudinale du passage 
des particules. En supposant que I'attenuation de la lumiere soit lineaire dans tout le cristal, 
la position du depot d'energie a une distance x de son centre est donnee par : 

^ = (3-1) 
(3L + R' ^ ^ 

avec |3 le coefficient d'attenuation de la lumiere par unite de longueur, et L et R respecti- 
vement la lumiere coUectee a I'extremite gauche et droite du barreau. La resolution sur la 
position de la gerbe est de I'ordre de quelques millimetres. Ces caracteristiques permettent 
ainsi de visualiser la gerbe electromagnetique en trois dimensions, et notamment d'aider a la 
reconstrution de la direction de la particule incidente. 

La reconstruction de I'energie totale de la gerbe est, quant a elle, obtenue par trois me- 
thodes distinctes, ayant des specificites pour chaque population de particules (photons, par- 
ticules chargees, basses et hautes energies, petits et grands angles d'incidence). Signalons que 
seuls les evenements reconstruits par le trajectographe sont traites, et que Ton corrige I'ener- 
gie deposee dans les parties inactives du calorimetre. La premiere methode prend en compte 
la profondeur du maximum de la cascade electromagnetique confinee dans deux barreaux, 
ainsi que I'energie du photon incident. La seconde s'appuie sur le profil longitudinal de la 
cascade mesure par les 12 couches, qui est ensuite ajuste a la forme analytique. Celle-ci ne 
sera utilise qu'a haute energie. Enfin, la troisieme methode utilise une simulation de gerbe 
qui est ensuite comparee au profil observe. Ces methodes sont presentees dans la these de 



D'Avezac (2006). Finalement, une des trois methodes sera selectionnee par un "Classification 



Tree" (voir la section 4.3) pour reconstruire I'energie de la particule incidente E. 

Comme dans le cas du trajectographe, un ensemble de ces variables est present dans les 
fichiers de pre-analyse. EUes nous servent entre autres a valider I'instrument en orbite et a 
identifier les photons. 



3.3.3 Detecteur Anti-CoVncidence 

Le but du detecteur anti-coincidence (ACD) est d'identifier les particules chargees, d'une 
part pour la rejection du fond de rayonnement cosmique (il joue le role de veto), et d'autre part 
pour I'etalonnage en energie du calorimetre avec les ions lourds du rayonnement cosmique. 
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Particule incidente 



Depot 



Csl(TI) ^ 



Luminescence 




Signal Lumineux 



r '\ Signal 

1/ Analoglque 



Detection 




Diodes 



Signal Digitise 

Electronique de Sortie 



Figure 3.5 - Schema de fonctionnement du calorimetre. La particule qui traverse le cristal 
depose de I'energie, par excitation des atomes du cristal. L'energie deposee se transforme en 
lumiere par scintillation. Ensuite, le signal lumineux se propage jusqu'aux diodes ou il est 



converti en signal analogique. Cette figure est extraite de iD'Avezac 2006). 





Caracteristiques 


Performance 


Profondeur, incluant le trajectographe 


10.1 Xo 


Stabilite des voies actives 


> 99% en 5 ans 


Gamme d'energie 


20 MeV - > 300 GeV 


Resolution en energie 


< 20% a 2 MeV < E < 100 MeV 




< 10% a 100 MeV < E < 10 GeV 


... e >60° 


< 6% a 10 GeV < E < 100 GeV 


Tableau 3.3 - Caracteristiques du calorimetre du LAT. 



Le detecteur est compose de 89 detecteurs scintillateurs plastiques, appeles tuiles (figure 
3.6), et couvre 8.3 m^ du detecteur. Les tuiles ont une taille d'environ 30 cm sur 30 cm, 



qui depend de leur emplacement. EUes ne recouvrent pas les faces du calorimetre, pour 
eviter des vetos dus a la retro-diffusion. Quand une particule chargee traverse I'ACD, elle 
depose de l'energie. La lumiere produite dans les tuiles est recuperee par des fibres optiques 
a decalage spectral. Cette caracteristique modifie la longueur d'onde de la lumiere pour une 
meilleure adequation avec les photomultiplicateurs, qui mesurent la lumiere. Pour alimenter 
les 18 photomultiplicateurs, le LAT dispose de chassis electroniques et d'alimentations haute 
tension ajustables, pour prendre en compte le vieillissement des photomultiplicateurs. 

Le taux de rayons cosmiques qui interagit avec le LAT est d'environ 10 kHz, compare aux 
quelques Hz de photons y. Le cahier des charges du LAT impose un bruit de fond residuel 
charge inferieur a 10% de I'intensite du fond gamma diffus. Cette condition impose une 
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Figure 3.6 - Schema d'un detecteur scintillateur plastique. Extrait de Johnson (2002) 



suppression d'un facteur 10^ pour les protons et IC* pour les electrons au niveau de la totalite 
de r instrument. La consequence est que I'efiicacite de detection des particules chargees par 
I'ACD doit etre de 0.9997, sans pour autant reduire le nombre de photons detectes. La 
topologie des evenements dans le trajectographe et la forme des cascades hadroniques dans 
le cas des protons vont aider a rejeter les protons et les electrons non detectes par I'ACD. 
Notons que le seuil de declenchement est de 0.3 MIP ("Minimum Ionizing Particle") sur chaque 
tuile. Finalement, la rejection du fond evite de stocker trop d'evenements compte tenu des 
contraintes telemetriques. 

Une cascade electromagnetique qui se produit dans le calorimetre pent produire des pho- 
tons jusqu'a des angles de 180° par rapport au sens de developpement de la cascade (voir 
la figure 3.7). Ces photons de seconde generation qui ont une energie de quelques keV vont 
interagir par effet Compton avec des electrons, lequels vont generer un signal dans I'ACD. 
Ce phenomene s'appelle la retro-diffusion ou "backsplash" en anglais. Le LAT, contrairement 
a EGRET, est compose d'un detecteur d'anti-coi'ncidence segmente, qui permet de correler 
I'alignement de la trace induite par la particule incidente avec la position les tuiles activees. 
Ce systeme limite les effets de self-veto dus a remission vers I'arriere et preserve la surface 
efiicace, notamment a haute energie. 



3.3.4 Declenchement 

L'instrument est compose de plusieurs niveaux de declenchement pour I'acquisition. L'en- 
semble est gere par le Global trigger Electronics Module (GEM). II utilise pour cela les 
performances du trajectographe, du calorimetre et du bouclier anti-coincidence ; le but etant 
d'avoir d'une part, un rapport signal-sur-bruit maximum en raison de la limitation de la 
bande passante telemetrique, et d'autre part, un temps mort minimum entre deux prises de 
donnees. Le satellite a la capacite d'enregistrer jusqu'a 36 heures de donnees. II transmet ses 
donnees 2 fois par jour a un debit de 3 ■ 10^ bits s^^ 

Le premier niveau (essentiellement declenche par les particules chargees) est compose de 
8 triggers individuels, et determine si la particule a les conditions requises pour etre analysee. 
La fenetre de declenchement est ouverte des qu'un des triggers est active. Elle est ajustee 
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Figure 3.7 - Illustration d'un phenomene de retro- diffusion ("backsplash") d I'interieur du 
LAT. Extrait de llJohnson\ \200^ 



en temps (entre 50 et 1550 ns) pour permettre aux (eventuels) autres triggers individuels 
d'arriver, mais suffisamment etroite pour eviter que d'autres evenements interagissent avec 



le LAT. Nous verrons qu'il peut y avoir des accidents (voir la section 4.5.2). Notons que 
chaque trigger est neutralise si le flux de donnees produit par un evenement est trop grand 
(par exemple, un trop grand nombre de pistes de silicium touchees), ou parce que le taux de 
comptage est trop important. Pour cela, des echelles sont associees a chaque trigger individuel 
pour diviser le nombre d'entrees des evenements a fort taux de comptage. La liste suivante 
presente ces 8 triggers individuels utilises dans le LAT au niveau du hardware. 

O TKR : Ce declenchement, dont la frequence moyenne est de 4 kHz, est donne par le 
trajectographe, si 3 plans x - y consecutifs sont touches. 

CAL-LO : Obtenu par le calorimetre, le systeme declenche la lecture des voies si un 
barreau de Csl mesure un depot de 100 MeV. II reduit le taux de declenchement TKR 
dans le cas ou un photon convertit toute son energie dans le trajectographe. 

© CAL-HI : Identique au declenchement de type CAL-LO, mais avec un seuil de 1 GeV. 

ROI : Certains evenements y peuvent activer les tuiles de I'ACD a cause de la retro- 
diffusion. Pour eviter de rejeter tous les evenements a forte emission vers I'arriere, un 
algorithme analyse le declenchement de type TKR et les tuiles touchees. 

© CNO : Active quand ~ 20MIP sont deposes dans une tuile de I'ACD. II permet la 
calibration du calorimetre avec des ions lourds. CNO : Carbone, Azote, Oxygene. 

© Periodic : Ce declencheur ouvre toutes les voies du detecteur. II permet de verifler les 
piedestaux de I'ACD et du calorimetre, et le taux d'occupation du trajectographe en 
vol. 

Solicited : Utilise pour le software du satellite. 

© External : Rien n'est connecte. Ne permet pas I'ouverture d'une fenetre. Historique- 



ment utilise pour les tests au sol (Bregeon, 2005). 



A la fermeture de cette fenetre, on enregistre I'etat de ces 8 conditions. Le second niveau 
de declenchement, le Global LAT Trigger (GLT), utilise la combinaison de ces triggers 
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individuels que Ton appelle "trigger engine". II contient les instructions de lecture et les 
configurations relatives au type d'evenement que Ton veut enregistrer (photons, ions). Par 
exemple, la selection des photons y demande que la condition TXR soit vraie, tandis que les 
autres doivent etre fausses. Le dernier niveau que Ton nomme le On Board Filter (OBF) 
est un filtre software, integre au systeme embarque d'acquisition des donnees. Son but est 
de rejeter un nombre suffisant d'evenements du fond, sans pour autant eliminer les photons 
y. Le grand avantage de ce filtre est qu'il est reprogrammable compare aux deux premiers 
niveaux. Le taux de particules apres I'OBF est d'environ 300 Hz. 

Finalement, la combinaison de ces trois niveaux triggers reduit le fond par un facteur de 
presque 10^, tandis qu'il preserve 75% des photons qui ont interagit avec le LAT. La suite 
de I'analyse se fait au sol et discrimine les photons et les particules chargees (voir la section 



4.3). 



3.3.5 Etalonnage en Energie du Calorimetre 

L'etalonnage en energie du LAT est un enjeu majeur, d'une part pour nos mesures de 
la coupure en energie des spectres des pulsars (discrimination des modeles), et d'autre part 
pour la nouvelle physique telle que la detection des signaux des particules de la matiere noire. 
Un premier etalonnage en energie a ete effectue au sol en utilisant les particules chargees et 
des injections de charges. La methode consiste a calibrer la gamme d'energie la plus basse 
(LEXSQ avec des muons cosmiques, qui laissent un depot d'energie d'environ 11.2 MeV par 
barreau, proche du minimum d'ionisation. Plusieurs millions d'evenements declenches par le 
calorimetre sont ainsi enregistres, et ensuite analyses. L'etalonnage au sol des autres voies se 
fait par inter-etalonnage grace aux injections de charge. 

L'etalonnage en vol utilise les rayons cosmiques. Un ion depose environ x 1 1 .2cos9 MeV 
au minimum d'ionisation dans chaque barreau traverse, avec 9 Tangle d'incidence. Ainsi, un 
proton permet d'etalonner la gamme d'energie la plus basse en deposant 11.2 MeV, tandis 
qu'un noyau de fer deposant environ 7.6 GeV permet d'etalonner la voie enregistrant les 
depots les plus energetiques (HEXl|^. Notons egalement que pour un meme depot d'energie 
la quantite de lumiere produite depend du type de particule. Cette non-linearite ou effet 
de "quenching" doit etre prise en compte. Une campagne de tests sous faisceau au GSl[^ a 
permis de mesurer ces non-linearites pour I'ensemble des ions et la gamme en energie utilisee 



pour l'etalonnage en vol (Lott et al, 2006). En resume, du proton jusqu'au fer, le spectre des 
rayons cosmiques couvre parfaitement les quatres gammes en energie, ce qui evite d'utiliser les 
injections de charges. Un algorithme ("GCRCalib") permet ensuite de garder les evenements 
bien reconstruits et rejette ceux qui ont interagit par reaction nucleaire. Finalement, comme 
le rayonnement est relativement intense, une journee est sufiisante pour obtenir une bonne 
statistique sur les ions legers, soit un taux de ~ 1.8 Hz pour les evenements CNO, et ~ 
10^^ Hz pour les evenements de fer. De plus, etant donne que les 4 gammes en energie ont un 
recouvrement deux a deux, il est possible d'inter-etalonner en vol les voies touchees par un 



3. LEX8 :Low Energy, grande diode et amplificateur a fort gain. LEXl : Low Energy, grande diode et 
amplificateur a gain faible. 

4. HEX8 :High Energy, petite diode et amplificateur a fort gain. HEXl : High Energy, petite diode et 
amplificateur a gain faible. 

5. http ://www. gsi.de/index. html 
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4.1 Introduction 



Le telescope LAT est congu pour repondre a des objectifs scientifiques tels que la compre- 
hension des mecanismes d'acceleration des particules dans la magnetosphere des pulsars. Ce 
travail passe par I'etude des spectres incluant les parametres enumeres a I'equation [1.23 et 
le flux integre. Pour acceder a ces quantites physiques, on utilise les fonctions de reponse de 
I'Instrument (IRFs, Instrument Response Functions), qui permettent de traduire une quan- 
tite reconstruite (energie E et direction 9 des photons) en une quantite physique vraie. Ces 



fonctions basees sur une simulation Monte Carlo (voir la section 4.2) sont gouvernees par 
la conception du LAT, des algorithmes de reconstruction des evenements (voir les sections 



3.3.1 et 3.3.2), et I'identification des particules (voir la section 4.3). Ce chapitre presente la 
validation de I'une des trois IRFs, la surface efficace. 
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4.2 Simulation de rinstrument 

La simulation detaillee du LAT et des evenements est faite par Gleam (GLAST LAT 
Event Analysis Machine), qui se base sur I'outil Monte Carlo GEANT40 ( [Agostinelli et al 



20031), logiciel de simulation des interactions entre les particules et la matiere. Cela inclut 



la diffusion multiple, la production de paires et la diffusion Compton pour les photons, le 
rayonnement de freinage pour les e^, aussi bien que les interactions hadroniques. L'informa- 
tion sur la geometrie du detecteur, la structure du satellite, et les materiaux utilises (pistes 
silicium du trajectographe, cristaux de Csl et les diodes du calorimetre, tuiles du bouclier 
anti-coi'ncidence) est stockee dans des fichiers XML. Notons que ces fichiers sont egalement 
utilises pour les algorithmes de reconstruction de I'energie et de la direction des evenements. 
La simulation est utilisee pour etudier les algorithmes de reconstruction et la classification 
des evenements, ainsi que les fonctions de reponse de I'instrument telles que la resolution 
angulaire, la resolution en energie, et la surface efficace. 

Avant le lancement du satellite, une caracterisation experimentale de I'instrument a ete 
effectuee, pour verifier que les reponses de I'instrument etaient en accord avec les predictions 
de la simulation Monte Carlo. II est essentiel d'obtenir une bonne reproductibilite entre 
les quantites mesurees directes (energie deposee ef'^'", activation des pistes x\}z^^, etc) et 
I'analyse haut- niveau (reconstruction de I'energie E et de la direction 9). De plus, une mesure 
directe des IRFs est requise pour certaines configurations critiques, telles que les photons de 
basse energie a cause de la diffusion multiple. Pour ces raisons, une campagne de tests sous 
faisceau a ete realisee au pres des lignes du PS^et du SPS^au CERN en ete 2006 ( |Baldini 



2007). Une premiere calibration avait egalement ete effectuee en utilisant des muons produits 
par le rayonnement cosmique dans I'atmosphere. Notons que toutes les configurations n'ont 
pas ete testees, en raison de I'espace de phase important du LAT en terme d'angle, de position, 
et d' energie. 



4.3 Reconstruction et Classification des Evenements 

La figure [47l] presente les differentes etapes de "I'analyse" des evenements du LAT, incluant 
la reconstruction et la classification des evenements. A leur arrivee sur Terre ou en sortie de 
Gleam, en d'autres termes qu'ils soient simules ou vrais, les evenements sont reconstruits et 
classes par la meme chaine d'analyse, ce qui rend les choses comparables ! 

La premiere etape notee (1) est la reconstruction de la direction et de I'energie de la 
particule incidente a partir des donnees brutes des sous-sytemes du LAT, comprenant le 
trajectographe (xyz^*^) et le calorimetre (ef'^'^ et xyzf'^'^). Les sections 



3.3.1 



et 



3.3.2 



res- 



pectivement decrivent brievement les algorithmes de reconstruction de la direction et de 



I'energie d'une particule incidente. On trouvera plus de details dans la these de D'Avezac 



(2006). L'analyse produit un fichier appele "recon" avec une multitude de variables utilisees 



1. http ://geant4.web.cern.ch/geant4/ 

2. PS : Faisceaux de photons "etiquetes", d'energie comprise entre 50 MeV et 3 GeV, et d'electrons et de 
protons entre 300 MeV/c et 15 GeV/c. 

3. SPS : Faisceaux d'electrons et de protons dans un domaine entre 10 GeV/c et 300 GeV/c. 
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GLEAM (GEANT4) 
Simulation des 
evenements et du LAT 
->sortie : digi 



VRAIES DONNEES 
Data Server 
-> sortie : digi 



Application des algorithmes pour la ' 
reconstruction de la direction (Kalman, etc) 
et de I'energie (cascade em, etc) de la particule. 
^ ^sortie : recon i 



(^Application des arbres de classification 
pour evaluer la meilleure direction, energie, et 
la nature de la particule a partir de I'etape 1. 
^sortie : merit 



^ Sous-ensemble des variables de 
I'etape 2 pour I'analyse temporelle et 
spectrale des sources. 

-»gnrrip : fits 



Figure 4.1 - Organigramme representant les differentes etapes de la reconstruction des 
evenements du LAT. fits (Flexible Image Transport System) est un format de fichiers utilise 
en astronomic : http :/ /fits. gsfc.nasa.gov 



pour I'analyse suivante. Notons que le developpement de la reconstruction repose fortement 
sur la simulation Monte Carlo des evenements. 

La deuxieme etape notee (2) associe a I'evenement une probabilite caracterisant le de- 
gre de confiance de la quantite mesuree. Cette analyse est basee sur les outils statistiques 
d'aide a la decision et a I'exploration des donnees, que Ton nomme les arbres de decisions 
(Classification Trees (CT), Breiman et al||1984 ). La construction de ces arbres a ete realisee a 
I'aide de la simulation Monte Carlo Gleam. Les particules qui passent a travers I'analyse sont 
reparties en group es selon un ensemble de variables, tels que la direction de la particule dans 
le trajectographe, I'emplacement de la trace, ou I'energie deposee dans le calorimetre. Cette 
analyse se deroule en trois etapes. Premierement, I'energie finale pour chaque evenement est 
assignee a I'energie determinee par la methode de reconstruction la plus appropriee (methode 
parametrique, correlation avec la derniere couche du calorimetre, estimation du profil de la 
gerbe). Celle-ci est selectionnee par un CT, que Ton nomme "CTBBestEnergy". Ensuite, un 
second CT permet d'estimer la probabilite que I'energie soit bien reconstruite ("CTBBes- 
tEnergyProb"). La deuxieme etape est une analyse similaire qui attribue a la particule la 
meilleure direction disponible ("CTBest[XYZ]Dir"), ainsi qu'une probabilite correspondante 
d'etre dans le coeur de la distribution de la PSF ("CTBCORE"). Finalement, quand I'energie 
et la direction de la particule ont ete evaluees, la probabilite que I'evenement soit un photon 
est examinee. La rejection du fond est une analyse difficile mais essentielle dans la recons- 
truction des particules, en raison de I'espace de phase important du LAT et du faible rapport 
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signal-sur-bruit (~l/300) des donnees entrantes. Cette derniere analyse est divisee en trois ni- 
veaux, chacun associe aux trois sous-sytemes du LAT (ACD, trajectographe, et calorimetre) . 
La probabilite que Tevenement soit un photon est evaluee pour chaque partie. Ces variables 
sont nommees "CTBCPFGamProb", "CTBCALGamProb", "CTBTKRGamProb"|5 

Finalement, trois classes d'evenements sont definies a partir de I'analyse decrite ci-dessus. 
Elles sont basees egalement sur le fond de particules chargees prevu en orbite et la connais- 
sance du ciel y. Ces classes se differencient par leur niveau de discrimination des photons, 
I'efiicacite, et la resolution angulaire et en energie. L'utilisateur pent egalement definir ses 
propres categories d'evenements selon le but de I'analyse, comme par exemple maximiser 
la resolution angulaire pour I'identification des sources, mais au prix de reduire la surface 
efficace. Les trois classes standards sont : 

O la classe transient : elle est definie avec les coupures les plus laches et presente un 
taux de fond residuel d'environ 2 Hz. Ceci equivaut en moyenne a un photon du fond 
toutes les 5 secondes a I'interieur d'une 10° autour de la source. La surface efficace est 
maximisee notamment a basse energie et done est appropriee pour I'etude des sources 
transitoires tels que les sursauts y. 

@ la classe source : elle a ete elaboree pour que la contamination du fond residuel soit si- 
milaire au flux extragalactique y estime par EGRET, soit un taux de photons d'environ 
0.4 Hz. Cette classe est appropriee pour I'etude des sources ponctuelles. 

© la classe diffuse : elle a la meilleure rejection de fond, soit un taux de photons d'environ 
0.1 Hz. Elle a ete elaboree pour etudier les sources diffuses les plus faibles. Notons que 
pour I'instant cette classe est utilisee pour I'analyse actuelle des pulsars y. 

Ces trois classes d'analyse definies pour la classification appelee "Pass6" sont documentees 



dans Atwood et al (2009). Pour I'analyse des sources de cette these, nous utilisons la version 
3 de cette classification, Pass6_v3. Contrairement a la version 1 utilisee pendant les premiers 
mois de la mission, la version 3 prend en compte les effets de traces fantomes observes pendant 



la phase de calibration de I'instrument (voir la section 4.5.2 et la figure 4.7). Un nouveau 
schema appele "Pass7" prevu pour la fin de I'annee 2009 inclut une analyse specifique des 
particules du fond, telle que les electrons, les protons, et les ions lourds. 



4.4 Fonctions de Reponse de I'instrument 

Les fonctions de reponse de I'instrument (IRFs) sont le lien entre le flux de photon vrai 
et les evenements detectes. Comme nous I'avons vu precedemment, la detection depend des 
caracteristiques de I'instrument, mais egalement de I'analyse qui assigne une probabilite que 
I'evenement soit un photon. Les IRFs peuvent etre formulees comme la surface multipliee par 
la probabilite, qu'un photon avec une energie vraie E' et une direction vraie 9' soit detecte 
aux quantites reconstruites E et 9. En negligeant les correlations entre I'energie et la direction 
de la particule et en supposant que la PSF est circulaire autour de la position vraie de la 

4. Cette page liste Tensemble des variables CTB present dans I'analyse : 

http ://www. slac.stanford.edu/exp/glast/ground/software/status/documentation/GlastRelease/latest/GlastClassify/latest/merittuple. html 
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source, les IRFs peuvent etre decrites par : 

R(E,e;E',e') = Aeff(E',e') PpsF(e;E',e') PE(E;E',e') (4.1) 

Le terme Agff (E',9') represente la surface efRcace differentielle, soit I'efiicacite de detection 
multipliee par la surface geometrique en fonction de I'energie et la direction vraie du photon 
incident. La dispersion de I'energie est la densite de probabilite qu'un photon avec I'energie 
E' et la direction 9' soit detecte avec I'energie reconstruite E. Tandis que Ppsp designe la 
fonction d'etalement du point (le terme en anglais est "Point Spread Function"), c'est-a-dire, 
la densite de probabilite qu'un photon d'energie E' et de direction 9' soit detecte avec la 
direction reconstruite 9. 

Les IRFs sont basees sur la simulation Monte Carlo Gleam decrite dans la section 14. 2[ Un 
ensemble de photons et particules chargees est genere isotropiquement sous tous les angles 
d'inclinaisons possibles et sur toute la gamme d'energie du LAT (de 20MeV a 300 GeV), et 
sont analyses separement par la procedure decrite dans la section |4.3[ La comparaison entre 
les parametres reconstruits de I'evenement et le photon entrant, dont on connait exactement 
la direction et I'energie, donne les IRFs. La simulation des particules chargees permet, quant 
a elle, d'estimer la contamination des differentes classes de photons. Les trois fonctions de 
reponse de I'instrument (resolution angulaire, resolution en energie, et surface efficace) sont 
presentees dans les sections suivantes. 

4.4.1 Resolution Angulaire 

La resolution angulaire du LAT est determinee par sa PSF, traduisant le pouvoir de lo- 
calisation d'une source. II s'agit de la precision avec laquelle une source d'emission y va etre 
reconstruite par le detecteur. La PSF est dominee par la diffusion multiple des premieres 
paires dans le trajectographe, principalement a basse energie (E<100 MeV). A haute ener- 
gie, la PSF est dominee par la resolution spatiale des plans de positionnement. Ceci est defini 
par le terme 5 dans I'equation suivante. A ces incertitudes s'ajoute celle sur le pointe du 
satellite (<40arcmin). De fagon generale, la PSF pent etre decrite par : 

PSF(E^) = © -p= © 5 , (4.2) 

avec Ey I'energie du photon incident. Le terme <x definit la dependance en energie de la 
diffusion multiple, tandis que |3 resulte de la distance entre I'entree du photon dans le LAT 



et la premiere conversion e^. La figure 4/2 gauche illustre pour la classe diffuse, la PSF du 
LAT a 68% et 95% de contenance en fonction de I'energie pour les differentes sections de 
conversions dans le trajectographe et pour une incidence normale, tandis que la figure droite 
montre la PSF en fonction de Tangle d'incidence pour des photons de 10 GeV. La PSF peut 
etre decrite, en se basant sur la correction de la formule de Moliere pour la diffusion multiple, 
par : 

PSF(Ey) = a X E7^ (4.3) 

avec E-y en GeV. Dans le cas de I'analyse "Pass6_v3" que Ton utilisera pour I'analyse des 
pulsars des chapitres 6 a 9, les termes a et b sont respectivement egaux a 0.8 et 0.75. 
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Figure 4.2 - PSF du L AT pour la classe d'evenements diffuse, Pass6_v3. Gauche : Angles 
d 68% (ligne continue) et 95% (ligne pointillee) de contenance de la direction reconstruite 
du photon entrant en fonction de I'energie, pour une incidence normale (cos 6 > 0.9), et 
pour des conversions dans la section mince (front), epaisse (back), mince et epaisse (total) 
du trajectographe. Droite : Angles a 68% (ligne continue) et 95% (ligne pointillee) de conte- 
nance de la direction reconstruite du photon entrant en fonction de I' angle d 'incidence pour 
des photons de lOGeV. 



4.4.2 Resolution en Energie 

L'energie reconstruite est le resultat des algorithmes de reconstruction du trajectographe 



et principalement du calorimetre. La figure 4.3 gauche montre pour la classe diffuse, la 



resolution en energie du LAT en fonction de l'energie pour une incidence normale. Une 
part de la cascade se developpe dans le trajectographe limitant la reconstruction par le 
calorimetre a toutes les energies, tandis qu'a haute energie (E > 30 GeV) seule la moitie de 
la gerbe electromagnetique (~ 50%) est presente dans le calorimetre (fuites vers I'arriere). 
Remarquons egalement que suivant I'endroit de la conversion (couche mince ou epaisse) 
et l'energie de la particule incidente, l'energie reconstruite est plus ou moins bien connue. 
Le calorimetre contient seulement la fin de la cascade electromagnetique pour un photon de 



basse energie converti dans les couches minces. La figure 4.3 droite represente la resolution 



en energie en fonction de Tangle d'incidence pour des photons de 10 GeV. 

4.4.3 Surface Efficace 

La surface efficace est la surface vraie acceptant les evenements bien reconstruits. Integree 
en fonction de l'energie E et de Tangle d'incidence 0, elle est le rapport entre le nombre de 
photons y detectes et acceptes par Tanalyse, et le flux incident y de la source pendant un 
temps d'observation. Ainsi de fagon simplifiee, le flux integre d'une source peut etre calcule 
par : 

= ^ ph cm-'s-i (4.4) 

" At X jAeff(E,e)dEde ^ ^ ' 

La surface efficace depend de la surface geometrique du detecteur, de la reconstruction, 



et de la selection des evenements suivants leur classification. La figure 4.4 gauche represente 
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Figure 4.3 - Gauche : Resolution en energie du LAT pour les evenements de la classe 
diffuse, classification Pass6_v3, en fonction de I' energie pour une incidence normale (cos > 
0.9), et pour des conversions dans la section mince (front), epaisse (back), mince et epaisse 
(total) du trajectographe. Droite : Resolution en energie en fonction de I'angle d'incidence 
pour des photons de lOGeV. 



pour les photons de la classe diffuse, la surface efficace en fonction de I'energie pour une 
incidence normale. Au-dessus de 100 MeV, en utilisant les photons convertis a la fois dans les 
sections minces et epaisses (total), la surface efiicace est compris entre 7000 et 8000 cm^. La 
figure droite represente la surface efiicace en fonction de Tangle d'incidence pour des photons 
de 10 GeV. 




Figure 4.4- Gauche : Surface efficace du LAT pour la classe diffuse en fonction de I'energie 
pour une incidence normale, et pour une conversion dans la section mince (front), epaisse 
(hack), mince et epaisse (total) du trajectographe. Droite : Surface efficace en fonction de 
I'angle d'incidence pour des photons de 10 GeV. 
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4.5 Validation de la Surface Efficace du LAT 

4.5.1 Methode d'Analyse 

La methode de validation de la surface efficace a ete developpee au SLAC[^ pendant les 
simulations de la reception des premieres donnees du LAT, entre octobre 2007 et mars 2008. 
Elle a ete ensuite appliquee aux vraies donnees a partir du 30 juin 2008 pendant la phase 
de calibration et de verification de I'instrument. Notons qu'elle est toujours utilisee par la 
collaboration pour valider les nouvelles versions de I'analyse de reconstruction des particules. 

Conceptuellement, la surface efficace des photons pent etre decrite comme : 

M 

Aj„(E,e)=Ageox];];e[(E,e), (4.5) 
avec Ageo la surface geometrique du detecteur (~22500 cm^), et M le nombre de niveaux dans 



I'analyse des evenements decrite en section 4.3 e^(E,9) est I'efficacite de la coupure (i) en 



fonction de I'energie E de la particule et de 9 Tangle entre I'axe normal au plan du LAT (z) 
et la direction du photon incident. On definit I'efficacite par : 

avec Ny[i) le nombre d'evenements restant apres la coupure (i). 

Pour valider la surface efficace Aj^^ des photons, nous proposons de comparer pour les 
donnees simulees et les donnees observees par le LAT, I'efficacite des coupures a chaque ni- 
veau (i) de la rejection du fond de particules chargees. Si les estimations ne sont pas en 
accord, nous determinons de la meme maniere quelles sont les variables responsables, telles 
que les CT associes au trajectographe, au calorimetre, et au detecteur anti-coincidence (sec- 



tion 4.5.3). La problematique avant de comparer I'efficacite des coupures est la selection des 
photons. Comment obtenir un echantillon de photons purs provenant d'une source observee 
par I'instrument ? Comment distinguer un photon d'une particule chargee observee par le 
LAT a chaque niveau de I'analyse? Le choix des sources est ainsi limite par la capacite a 
selectionner des photons. II est egalement necessaire d'avoir une source assez brillante pour 
I'observer avant la rejection du fond hadronique. Les pulsars permettent de selectionner un 
echantillon pur de photons y en utilisant leur information temporelle. De plus, le pulsar de 
Vela est la source la plus brillante du ciel y. 



La figure 4.5 presente le phasogramme (ou courbe de lumiere) du pulsar de Vela pour deux 
coupures (rejection du fond) differentes. Rappelons que la courbe de lumiere d'un pulsar est 
le nombre d'evenements en fonction du temps, modulo la periode. Nous reviendrons sur cette 
notion dans le chapitre |5] La figure gauche montre la courbe de lumiere observee apres le 
filtre OBF et I'application d'une coupure minimale pour une region d'integration (ROI) de 
10° autour de la position du pulsar, tandis que la figure droite montre pour la meme ROI 
la courbe de lumiere de ce meme pulsar pour la classe d'evenements diffuse. Les coupures 



appliquees sont enumerees dans la table 4.1 Ainsi, le nombre de photons pent etre estime 



5. SLAC National Accelerator Laboratory, Stanford, Californie. 
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par la soustraction du nombre d'evenements de la region OFF normalisee (0.64 < ^ < 1.0), 
moins le nombre d'evenements de la region d'emission des pics (region ON). Le premier pic 
est defini entre 0.08 < cj) < 0.18, tandis que le second pic est defini entre 0.5 < (j) < 0.62. 
Cette methode s'appelle la technique ON - OFF. Le nombre de photons Ny est done egal a 
Npi/0.1 + Np2/0.12 — Noff/0.36, avec Npi et Np2 le nombre d'evenements dans les pics et Ngff 
le nombre d'evenements dans la region OFF. 

Notons que I'information temporelle du pulsar (sa phase) n'est pas inclus dans les fichiers 
de donnees. Pour ce travail, nous avons ainsi cree un outil nomme "addPhase" qui permet 
d'ajouter dans les fichiers bas-niveau (fichiers merit) et haut-niveau (fichiers fits) la phase 
du pulsar en se basant sur ses parametres de rotation (voir le chapitre |5]). 
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Figure 4.5 - Gauche : Courbe de lumiere pour le pulsar de Vela apres le filtre OBF et 



V application d'une coupure minimale (voir la table 4-.1). La partie grisee correspond aux deux 



pics du pulsar (region ON ). La ligne pointillee indique le debut de la region OFF (0.64 < ^ < 
].0). Droite : Courbe de lumiere du pulsar de Vela pour les evenements "diffuse" (coupure 
finale). Notons que le fond residuel passe d' environ 5000 coups par intervalle de phase a 
une vingtaine de coups en appliquant les coupures standards. Dans le mime temps, le signal 
passe d'une significativite ( signal/ V fond ) d' environ 5a a 20<j. 



4.5.2 Validation a Haut Niveau 

La validation a haut niveau correspond a I'analyse des efficacites pour les classes d'evene- 



ments transient, source, et diffuse, definies dans la section 4.3 Pour illustrer ce travail, nous 
presenterons dans cette these uniquement la validation de la classe "Diffuse" utilisee pour 
I'analyse des pulsars. L'efficacite des coupures e^(E,6) pour la classe diffuse est : le nombre 
de photons estime par la techinique ON - OFF et passant I'analyse diffuse (rejection du fond 
et reconstruction des evenements), sur le nombre de photons apres la coupure minimale (voir 



I'equation 4.6). Etant donne I'important espace de phase du LAT, l'efficacite de coupures sur 
les donnees LAT et Monte Carlo est evaluee suivant un nombre maximum d'intervalles en 
energie et en angle, lequels sont limites par la statistique des donnees LAT. 



77 



Chapitre 4. Validation de la Surface Efficace du Large Area Telescope 



Nom de la coupure Definition 



Coupures minimes TkrNumTracks>0 && CalCsIRLn>4 &:& CalEnergyRaw>5 



Transient 

Source 

Diffuse 
ROI 



CTBCORE>0. && CTBBestEnergyProb>0. && CTBBestEnergy>10 
&& CTBBestEnergyRatio<5 && CTBClassLevel>0 

CTBCORE>0.1 && CTBBestEnergyProb>0.1 && CTBBestEnergy>10 
&& CTBBestEnergyRatio<5 && CTBClassLevel>l 

CTBCORE>0.1 && CTBBestEnergyProb>0.1 && CTBBestEnergy>10 
&& CTBBestEnergyRatio<5 && CTBClassLevel>2 
(FTlDec-dec)2+(cos(FTlDec*7T/180)*(FTlRa-ra))^<100 



Tableau 4.1 - Definition des coupures standards pour I'analyse des evenements du 
LAT (Transient, Source, Diffuse) et de la coupure minimale permettant d' avoir une 
reconstruction d 'evenement. 



Notes : 

- TkrNumTracks > demande au moins une trace dans le trajectographe. 

- CalCsIRLn > 4 impose que la trace intercepte le calorimetre (>4Xo). 

- CalEnerguRaw > 5 impose au moins 5 MeV d'energie deposee dans le calorimetre. 

- CTBBestEnergTjRatio est le rapport entre la meilleure energie estimee ("CTBBestEnergy") 
et I' energie deposee dans le calorimetre. 

- CTBCLassLevel est une variable qui classe les evenements selon le degre de confiance qu'il 
soit un photon (O=fond hadronique, l=transient, 2=source, 3=diffuse). 

- FT! Dec et FTIRa sont les coordonnees de I'evenement. Les coordonnees du pulsar Vela 
utilisees sont (ra,dec) = (-45.176354° ,128.835881° ). 



La figure 4^ montre pour 8 intervalles en angle defini par "CTBBestZDir"[j, I'efficacite 
y pour la classe d'evenement diffuse en fonction de I'energie definie par la variable "CTB- 
BestLogEnergy" pour la classe diffuse. Les points noirs correspondent aux donnees du LAT, 
tandis que les courbes rouges et bleues correspondent a deux echantillons de photons simu- 
les. La difference entre ces deux simulations est la prise en compte pour la simulation vl7r7 
(courbe rouge) des traces fantomes (figure 4.7). On notera que pour le deuxieme intervalle en 
angle (-0.9< CTBBestZDir <-0.8), I'efficacite a partir de quelques GeV se degradent. L'effet 
est d'environ 2 a, et nous n'en connaissons pas la raison exacte. Dans la selection diffuse, 
il se peut qu'il y ait une dependance en angle que Ton ne prenne pas en compte dans I'analyse. 



Dans les donnees Monte Carlo, chaque evenement est genere independamment, et les 
interactions entre ces evenements ne sont pas pris en compte. Plusieurs cas de ce qu'on 
appelle les traces fantomes ont ete observees dans les donnees reelles pendant la periode de 
mise en route. La plupart des evenements du fond hadronique sont facilement identifies en 
tant que tels, et les triggers individuels du LAT ne sont pas actives pour enregistrer ce type 
d'evenement. Cependant, il arrive qu'un photon interagissant avec le LAT ouvre la fenetre 
d'acquisition des donnees, et au meme moment une particule chargee entre dans le detecteur. 
II est egalement possible que le signal d'un ion soit encore present lors du declenchement du 
detecteur par un photon. Les traces laissees par la particule chargee sont ainsi enregistrees. 

6. Cette variable correspond a —cos 9, avec 6 Tangle entre I'axe normale au plan du LAT et la direction 
incidente du photon. 
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Figure 4.6 - Efficacite des coupures pour la classe diffuse en fonction de I'energie ("CTB- 
BestLog Energy") et de I'angle d'incidence du photon. CTBBestZDir correspond d —cos 6, 
avec 9 I'angle entre I' axe normale au plan du LAT et la direction incidente du photon. 
Les points noirs correspondent aux donnees du LAT (incertitude statistique), tandis que les 
courbes rouge et bleues sont des donnees simulees. (Credit : P.Bruel). 
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La figure 4/7 montre un exemple de trace fantome. La consequence est qu'un certain nombre 
de photons bien reconstruits sont identifies comme des mauvais evenements lors de I'analyse 
de la reconstruction au sol. Ces evenements sont rejetes et conduisent a une degradation des 
performances du LAT et notamment de la surface efficace. Ainsi, pour etre en accord avec les 
donnees reelles, la simulation vl7r7 inclut aleatoirement des faux evenements pour prendre 
en compte ces problemes de coincidence. 




Figure 4.7 - Exemple d'un evenement fantome dans le LAT. Sur la droite, un photon 
interagit dans le trajectographe et dans le calorimetre. Sur la gauche, il apparait une trace 
fantome dans le trajectographe resultant d'une coincidence temporelle accidentelle d'une 
particule chargee. Dans le trajectographe et le calorimetre, seules les parties actives sont 
montrees. Les traces candidates sont egalement dessinees. 



La figure |4.9| presente le rapport entre les vraies donnees et les donnees simulees pour 



les memes intervalles en energie et en angle que la figure 4.6 Les points bleus resultent de 
la comparaison entre les vraies donnees et la simulation vl5r51pl qui n'inclut pas les effets 
d'evenements fantomes, tandis que les points noirs utilisent les donnees simulees prenant en 
compte les traces fantomes. On remarque clairement que la derniere simulation reproduit 
mieux les donnees observees par le LAT. A partir de ces resultats, une enveloppe correspon- 
dant aux incertitudes systematiques induite par la surface efficace a ete estimee. Les resultats 
sont resumes dans la section I4.5.4[ 

4.5.3 Validation des Variables Individuelles 

Pour verifier la procedure de validation a haut-niveau, les principales variables de la clas- 
sification sont egalement analysees. EUes permettent d'affiner la validation et notamment de 
verifier si il n'y a pas un defaut sous-jacent aux trois classes de photons. Nous ne presenterons 
ici que la variable CTBCALGamProb qui represente la probabilite estimee par le calorimetre 
que la particule soit un photon. 



La figure 4.8 (gauche) montre I'efficacite de la variable CTBCALGamProb par intervalle 
de 0.1 entre les donnees du LAT et les donnees simulees vl5r21. Les points bleus representent 
I'efficacite des donnees Monte Carlo, tandis que les points rouges representent les donnees 
LAT avec la prise en compte de I'incertitude statistique. A droite, la deviation de la variable 
est representee. Malgre I'absence de prise en compte des effets fantomes, elle est inferieure a 
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CTBCALGamProb ■ MC(blue) ■ LAT(redt 




CTBCALGamProb 



Deviation : CTBCALGamProb 



■a 
5 



%' ; ndf 

po 



-0.01644 + 0.03264 



CTBCALGamProb 



Figure 4.8 - Gauche : Efficacite de la variable CTBCALGamProb par intervalle de 0.1. Les 
•points bleus representent I' efficacite des donnees Monte Carlo vl5r21. Cette version de la 
simulation ne prend pas en compte les effets fantdmes. Les points rouges representent les 
donnees LAT. Droite : Deviation de la variable CTBCALGamProb. 



20%. La majorite des points (excepte le point 0-0.1) sont en accord aux incertitudes pres 
avec une deviation egale a 0. Ces variables sont evaluees a chaque mise a jour de I'analyse 
de la classification des evenements. Rappelons neanmoins que le but est d'avoir une incerti- 
tude systematique sufiisamment petite sur les quantites mesurees, tels que le flux et I'indice 
spectral des sources. 



4.5.4 Estimation des Incertitudes Systematiques 

La figure [4^ montr e a partir du rapport des efRcacites, I'estimation des incertitudes syste- 
matiques induite par la surface efficace. L'enveloppe superieure representee en noire associee 
a la classification Pass6_v1 montre une erreur de 20% en dessous de 100 MeV, environ 10% 
en dessous de 1 GeV, et de 30% au-dessus de 10 GeV. L'enveloppe rouge (inferieure) associee 
a la classification Pass6_v3, qui prend en compte les effets de traces fantomes, montre une 
erreur de 10% en dessous de 100 MeV, environ 5% en dessous de 1 GeV, et de 20% au-dessus 
de 10 GeV. A partir de ces resultats, deux echantillons d'IRFs ont ete crees pour estimer 
les effets maximaux de la variation de la surface efficace. lis ne prennent pas en compte les 
erreurs sur la PSF et sur la resolution en energie. Ces deux echantillons sont ensuite utilises 
pour evaluer les erreurs systematiques sur le fiux y integre, I'indice spectral F, la coupure en 



energie Ec, et I'indice d'attenuation (3 des pulsars (voir I'equation 1.23). 



Nous verrons aux chapitres 6, 7, et 8 (resultats sur les pulsars individuels) , comment 
est evaluee I'incertitude systematique sur les parametres spectraux et notamment I'indice 
d'attenuation |3, parametre qui permet de discrimer les differents modeles d'emission y des 
pulsars. Brievement, la methode consiste a ajuster les parametres aux donnees en utilisant 
premierement I'lRF standard, et ensuite les deux IRFs modifiees pour estimer un intervalle 
de confiance. 
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Figure 4.9 - Rapport des efficacites entre les vraies donnees et les donnees Monte Carlo. Les 
points bleus resultant de la comparaison entre les vraies donnees et la simulation vlSrSlpl 
n'incluant pas les effets d 'evenements fantdmes. Les points noirs resultent de la comparaison 
entre les vraies donnees et la simulation vl 7r7 prenant en compte les effets de trace fan- 
tome. Les courbes rouges et noirs representent I'estimation des systematiques pour la version 
initiale de la procedure initiale d'analyse des evenements (P6_vl) et la version amelioree 
(P6_v3) prenant en compte les effets de trace fantdme. (Credit : Ph.Bruel). 
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Chapitre 5 



Analyse des Donnees 



Sommaire 



5.1 Introduction 

5.2 Analyse Temporelle 

5.2.1 Barycentrisation des Evenements 

5.2.2 Les Ephemerides 

5.2.3 Construction d'une Courbe de Lumiere 

5.2.4 Significativite 

5.3 Analyse Spectrale Resolue en Phase . 

5.3.1 Introduction 

5.3.2 Selection des Donnees 

5.3.3 Mo dele des Sources 

5.3.4 Maximum de Vraisemblance . . . . 



5.1 Introduction 

Ce chapitre decrit les analyses temporelles et spectrales des pulsars y a partir des don- 
nees du LAT[^ des outils standards developpes par la collaboration Fermi (les ScienceTools), 
ainsi que du programme d'analyse des pulsars radio TEMP02 ( ]Hobbs et al[ 2006). Une pre- 



miere partie presente la construction d'une courbe de lumiere (ou phasogramme, histogramme 
a n intervalles en phase) a I'aide d'une ephemeride de pulsar, et une seconde partie presente 
I'analyse spectrale resolue en phase d'un pulsar en prenant en compte les sources avoisinantes 
et le fond d'emission diffus. 

Notes sur les ScienceTools (ST) : 
Les Sciences Tools sont un ensemble d'outils utilises pour I'analyse des donnees haut niveau 

1. Les donnees sont publiques depuis le 19 aotit 2009 et peuvent etre telechargees ici : 
http ://fernii. gsfc.nasa.gov/cgi-bin/ssc/LAT/LATDataQuery.cgi 
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de Fermi, donnees stockees dans des fichiers an format fits (voir le schema 4.1). Les outils 
sont publiques, developpes en langage C++ et Python, et se repartissent en sous-groupes 
chacun associe a differents types d'analyse : I'extraction et la selection des donnees, I'analyse 
spectrale des sources ponctuelles, les sursauts y, I'analyse des pulsars, la simulation, etc. Le 
lecteur desireux de debuter sa premiere analyse trouvera une documentation complete et mise 
a jour des ST sur : http ://glast-ground. slac.stanford.edu/workbook/sciTools_Home. htm et 
http ://fermi. gsfc.nasa.gov/ssc/ . 



5.2 Analyse Temporelle 

Signalons au prealable, que le temps d'arrivee des evenements du LAT doit etre correct 
et precis afin que cette etude soit menee a bien. Des tests sur les horloges du LAT ont ete 
effectues avant et apres la mise en orbite du satellite pour s'assurer de leur bon fonctionne- 



ment (Abdo et al 2009-calib). Au sol, les tests ont consiste a mesurer la difference des temps 
d'arrivee des muons entre les systemes GPS du LAT et un systeme GPS de reference, tandis 
qu'en orbite nous avons utilise la position des pics y des pulsars les plus brillants. 



Pour construire la courbe de lumiere d'un pulsar (ou phasogramme, voir la figure 1.4), 
il faut determiner a quel stade de la rotation de I'etoile a neutrons le photon detecte par le 
telescope a ete emis. En d'autres termes, il faut attribuer a chaque photon detecte une phase 
rotationnelle. Signalons que I'abcisse d'un phasogramme est tres souvent definie en tant que 
phase du pulsar. II s'agit en fait de la phase normalisee a 2n, c'est-a-dire la fraction de tour 
de I'etoile a neutrons que Ton represente entre et 1. 

La premiere etape est de remettre en ordre le temps des evenements, en transferant les 
temps d'arrivee du photon observe par le LAT vers le barycentre du systeme solaire (section 



5.2.1), ceci pour s'affranchir des mouvements propres de I'observatoire et du pulsar. Ensuite, 



a I'aide de ces temps barycentres et des parametres de rotation du pulsar decrits par une 



ephemeride (section 5.2.2), une phase est assignee a chaque photon. Dans cette these et dans 
une grande partie des publications scientifiques de la collaboration Fermi, la construction 
des phases pour les pulsars est obtenue avec le programme d'analyse TEMP02 et le ST 
"gtpphase". Notons cependant que I'outil gtpphase ne pent prendre en compte que les deux 



premieres derivees de la frequence pour calculer la phase (voir I'equation 5.2). Etant donne 
que les pulsars que Ton etudie sont relativement jeunes et presentent une frequence de rotation 
instable {timing noise et glitch), il est souvent necessaire de recourir a la troisieme voire la 
quatrieme derivee de la frequence. Cette situation m'a amene a utiliser uniquement I'outil 
TEMPO! pour eviter toute confusion dans I'analyse. 

5.2.1 Barycentrisation des Evenements 

Pour analyser le temps des evenements emis par le pulsar, ou en d'autres termes suivre sa 
rotation de fagon precise, il est necessaire idealement de se situer dans un referentiel au repos 
par rapport a I'objet. Cependant, notre referentiel d'observation n'est pas inertiel vis-a-vis 
du pulsar, que ce soit pour des telescopes situes sur Terre en mouvement autour du soleil, ou 
pour le satellite Fermi en orbite autour de la Terre. Le meilleur referentiel est le barycentre du 
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systeme solaire (SSB, Solar System Barycenter), point proche du soleil. La "barycentrisation" 
consiste ainsi a transferer des temps d'arrivee topocentriques (enregistres par un observatoire) 
vers le centre de masse du systeme solaire qui varie suivant la position des astres dans le 
systeme solaire. Ceci revient a dater un evenement suppose observe au barycentre du syteme 
solaire, et permet d'avoir, pour tous les evenements, une echelle de temps homogene notee 
TDB, pour "Temps Dynamique Barycentrique". Ce temps barycentre tssB est calcule par 
TEMP02 en fonction du temps tMET tt de I'evenement enregistre par le LAT. Les temps 
"tjviET TT sont exprimes en Mission Elapsed Time (MET) TT (Temps Terrestre)[^ Le temps 
tssB est calcule de la fagon suivante : 

tsSB — "tjVlET TT + ^Einstein + ^Roemer + ^Shapiro • (5-1) 

O AEinstein '■ cette eclielle regroupe la dilatation du temps introduite par le mouvement 
de la Terre dans le potentiel gravitationnel du systeme solaire, ainsi que le decalage 
vers le rouge gravitationnel induit par les corps du systeme solaire, et permet de passer 
du referentiel geocentrique (temps en TT) au referentiel barycentrique (temps TDB). 
Les positions du soleil, de la Lune, et des 8 planetes de notre systeme solaire sont 
definies dans des ephemerides du systeme solaire DE200 ou DE405 publiees par le 
JPL0 ( |Standish[ |1998D . 



ARoemer : cet effet geometrique correspond au decalage separant I'arrivee des photons 
du LAT par rapport au barycentre du systeme solaire. 

© Ashapiro '■ cet effet est une correction temporelle induit par le passage proche des photons 
dans I'espace-temps courbe d'un objet massif, tel que le Soleil. 

Pour plus de precisions sur la "barycentrisation" ainsi que sur les echelles de temps, on se 



reportera a I'excellent ouvrage de |Lorimer and Kramer (2005). 



5.2.2 Les Ephemerides 

Les ephemerides des pulsars decrivent les parametres de la dynamique de rotation des 
pulsars en fonction du temps (i.e. la frequence de rotation et ses derivees) dans le referentiel 
du barycentre du systeme solaire (TDB), ainsi que les parametres astrometriques necessaires 
pour la barycentrisation, tels que la position du pulsar, le mouvement du pulsar dans le plan 



du ciel, et les parametres orbitaux pour les systemes multiples. Rappelons que la section [2.3.1 
decrit brievement la campagne de chronometrie des pulsars qui fournit les ephemerides. En 
pratique, une ephemeride est construite en minimisant I'ecart entre la phase observee, calculee 
a partir des temps d'arrivee des photons (TOAs, Times Of Arrival) d'une observation donnee, 
et la phase calculee a partir d'un modele de rotation du pulsar, a savoir les parametres decrits 
precedemment (la position, la frequence de rotation, etc). Le calcul de la phase d'un pulsar a 



partir du temps d'arrivee des photons est donne par 1 'equation ] 5. 2 [ La plupart des parametres 
de I'ephemeride sont etablis pendant la decouverte du pulsar, tandis que la frequence et ses 
derivees sont reajustees a mesure que les observations progressent. L'observation d'un pulsar 



2. Le temps TT est le temps que mesurerait une horloge parfaite, situee au niveau de la geoide : surface 
equipotentielle coincidant avec le niveau de la mer en premiere approximation. 

3. Jet Propulsion Laboratory 



85 



Chapitre 5. Analyse des Donnees 



pour un radiotelescope peut durer de quelques minutes a quelques heures, cela depend de 
I'intensite du signal radio. 



La figure 5.1 presente les residus des phases (phase predite moins phase observee) pour 



trois pulsars observes par GBT, dont PSR J2229+6114 presente au chapitre [7} Remarquons 
que ces residus sont disperses autour de 0. Cette dispersion peut s'expliquer principalement 
par du Timing Noise, par I'incertitude des observations radio, ou encore par des effets non 
pris en compte dans Tajustement de la dynamique de la rotation du pulsar. II est done 
essentiel pour certains pulsars tres bruites (i.e. un timing noise important), notamment 
pour les pulsars jeunes tels que le pulsar du Crabe, de les observer regulierement. 
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Figure 5.1 - Residus des phases •pour PSR J1747-2958, PSR J1930+1852, et 
PSR 2229+6114, observes par le Green Bank Telescope a 2 GHz (Credit : F.Camilo). 



La Table presente les principaux parametres de I'ephemeride pour le pulsar J2229+6114 
basee sur les observations des telescopes GBT et Jodrell Bank. Cela inclut : la position du 
pulsar fixee aux coordonnees des observations X (Halpern et al 2001b), la frequence et ses 
derivees, une epoque Tq de reference (point central de I'ajustement du developpement de 
Taylor, voir I'equation 5.2), une date de reference ("TZRMJD") pour laquelle la phase est 
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egale a 0, Tintervalle de validite de rajustement, le nombre de TOAs utilises pour construire 
I'ephemeride, la mesure de dispersion, et le modele d'ephemeride du systeme solaire utilise 
pour la "barycentrisation". 



Parametres Valeurs 

Ascension Droite (J2000) 22^29"^05.28000084'^ 

Declinaison (J2000) +6ri4'09.3000042" 

Frequence (s^^) 19.363953152(1) 

Derivee premiere de la frequence (s^^) -2.92006(5) xlO^^^ 

Derivee seconde de la frequence (s^^) 1.50(8) x 10^^^ 

Epoque, To (MJD[TDB]) 54774.0 

Date de reference pour le 4)o (MJD[TDB])'' 54780.986659550693183 

Intervalle de validite en MJD 54634.286 - 54913.622 

Nombre de TOAs 69 

Mesure de dispersion, DM(cm"^pc) 204.97 ± 0.02 

Modele d'ephemeride du S.S DE405 



Tableau 5.1 - Ephemeride radio du pulsar PSR J2229+61 14 basee sur les observations 
des telescopes Green Bank Telescope (GBT) et Jodrell Bank. Les chiffres entre parentheses 
representent les incertitudes donnees par TEMP02. 
Date a laquelle la phase c}) = (TZRMJD). 



Les ephemerides radio et X utilisees pour la detection de pulsations y avec le LAT sont 
stockees sur les serveurs du CENBG Bordeaux-GradignanQ et de la NASAQ 

5.2.3 Construction d'une Courbe de Lumiere 

Soit 4)(t) la phase du pulsar a I'instant t. En utilisant les parametres d'une ephemeride, 
la courbe de lumiere est I'histogramme de la phase calculee a partir des temps d'arrivee des 
photons, prealablement transformes dans le referentiel du centre masse du systeme solaire 
suivant le developpement de Taylor : 

j=0 



ou To est la date de reference de I'ephemeride du pulsar (voir la table 5.1 ), fj est la derivee de 
la frequence de rotation a I'ordre j, et (j)o est la phase de reference. Cette derniere est donnee 
par la date de reference ("TZRMJD") pour laquelle la phase est egale a : 

*o = -Z (5-3) 

4. http ://www. ceiibg.iii2p3.fr/ephein/ 

5. http ://fermi.gsfc.nasa.gov/ssc/data/access/lat/epheins/ 
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5.2.4 Significativite 



A partir des evenements barycentres et empiles suivant Tequation 5.2, la distribution de 



phases est ensuite analysee pour rechercher un eventuel signal pulse a la frequence radio 



attendue. Dans le cas du pulsar de Vela (voir la figure droite 4.5), la mise en evidence de la 
pulsation est claire. On pent simplement, pour estimer la significativite du signal, compter 
le nombre de coups pulses au-dessus du fond que Ton compare au nombre de coups du fond 
(signal/ y/ fond) . Cependant dans certains cas, la fraction pulsee pent etre beaucoup plus 
faible et le fond d'emission y plus important. Les tests statistiques suivants sont codes et 
utilises pour rechercher un signal periodique a partir d'une courbe de lumiere. 

Test Statistique 

De fagon generale, le test mesure la signification de I'ecart entre des distributions 
observees et des distributions attendues, soit : 

X^ = i^^^. (5.4) 

1=1 ^ 

ou N est le nombre d'intervalles, est le nombre d'evenements dans un intervalle i, est la 
distribution attendue de coups, et a? est la variance associee. Dans le cas d'un phasogramme, 
le test est utilise pour tester I'hypothese selon laquelle la distribution de phases est plate, 
c'est-a-dire, que les phases sont uniformement distribuees, soit \Li = (x), avec (x) la valeur 
moyenne des x^. Pour un nombre sufEsamment grand d'evenements, la distribution des Xt suit 
la loi de Poisson, ainsi 0? est equivalent a (x). Finalement, le test x^ pour un phasogramme 
est defini par : 

= i: ^^^^^^ > (5-5) 
La consequence est que pour un fort signal pulse, la deviation par rapport a une distribution 



plate sera importante, ainsi que le resultat du x^- On trouvera dans Leahy et al (1983) la 
probabilite associe a un x^ d'observer une distribution plate. Notons que ce test depend du 
nombre d'intervalles utilise pour construire la courbe de lumiere. 

Test Statistique 

De la meme fagon que le test x^, (Buccheri et al 1983) s'appuie sur la mesure de 



la difference entre une distribution plate et mesuree. II se base sur le calcul du spectre de 
puissance de Fourier de la distribution de phases, en utilisant les m premieres harmoniques. 
De plus, il est independant du choix du nombre d'intervalles du phasogramme construit. 
Ainsi, la deviation de la distribution de phase mesuree par rapport a une distribution plate 



est calculee a partir des Ni valeurs de phases obtenues par I'equation 5.2 (provenant des Ni 
temps d'arrivee des photons), soit : 

' k=l j=l 5=1 

La probabilite de Z^ dans le cas d'une absence de signal est en fait une probabilite de x^ 
avec 2m degres de liberte. Notons que suivant le choix de I'harmonique m choisie, le test est 
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plus ou moins performant selon la forme du pic. Pour les pics etroits une harmonique elevee 
est plus appropriee, tandis que pour des pics larges une harmonique petite est preferable. II 
faut done avoir une idee a priori de la forme de la courbe de lumiere, ce qui n'est pas tres 
satisfaisant. 

Test Statistique H-test 

Base sur le Z^, De Jager et al ([1989") proposent le test statistique H-test qui calcule la 



valeur optimale m de 1' harmonique. L'avantage de ce test est qu'il n'est pas biaise par la 
forme du pic. L'harmonique optimale correspond a I'harmonique m comprise entre 1 et 20, 
pour laquelle la valeur H est maximisee, soit : 

H = Z^-4m + 4. (5.7) 

L'expression est applicable pour un nombre d'evenements superieur a 100. Finalement, la 
probabilite associee a H d'obtenir une distribution plate est donnee par : 

0.9999755e-°-^'^^°^^ si < H < 23 

P(H) = { i.210597e(-°-4590iH+o.oo29H2) si 23 < H < 50 (5.8) 

<4xlO-^ siH>50 



5.3 Analyse Spectrale Resolue en Phase 

5.3.1 Introduction 

L'analyse spectrale des pulsars (sources ponctuelles) se base sur une liste de photons 
detectes par le LAT a I'interieur d'une region d'integration (ROI) qui est la region d'extraction 
des donnees en fonction de leur energie, de leur direction, et de leur temps d'arrivee. Cette 
liste de photons provient de la reconstrution et de la classification des evenements a partir 



des donnees mesurees du LAT (voir la section 4.3). Cependant, le faible nombre de photons 



observes pour une source ponctuelle apres I'application des coupures d'une classe de photons 
donnee signifie que les methodes statistiques telles que les tests ne sont pas applicables 
pour evaluer les parametres spectraux du pulsar. II est necessaire d'estimer ces parametres par 
un maximum de vraisemblance en supposant une distribution discrete des evenements (loi de 
Poisson). L'estimation a pour but de trouver les meilleurs valeurs possibles des parametres 
pour que la fonction de distribution modelisee des evenements s'ajuste le mieux aux donnees 
disponibles. Dans notre cas, les parametres des sources sont le flux integre et les parametres 



spectraux du pulsar (voir I'equation 1.23). Finalement, les termes "resolue en phase" signiflent 
que les photons sont selectionnes dans un intervalle de phase deflni. A condition que la 
statistique soit sufflsante, ceci permet d'etudier, par l'estimation des parametres spectraux, 
les differents composants d'une courbe de lumiere, tels que les pics ou la region entre les 
pics. Avant de presenter la methode, il est essentiel de rappeler quelques points importants 
speciflques au telescope qui doivent etre pris en compte dans l'analyse : 

O Le LAT ne flxe pas un point flxe du ciel comme le telescope EGRET et les telescopes 
X, il opere selon un mode de balayage lui permettant d'observer le ciel entier en deux 
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orbitesj^ La consequence est que les sources ne sont pas observees continument en temps 
et avec le meme angle d' incidence. 

@ Le champ de vue du LAT est important (2.4 sr). De nombreuses sources sont vues en 
meme temps. 



@ La PSF du LAT est relativement large (voir la figure 4.2). De nombreuses sources, 
notamment au niveau du plan galactique, se recouvrent spatialement. De plus, le fond 
d'emission diffus du plan galactique est fortement structure. II est done impossible, 
sauf eventuellement a haute latitude galactique, d'analyser une source ponctuelle isolee. 
L'analyse implique done I'ajustement simultane de plusieurs sources ponctuelles et du 
fond d'emission diffus. 

Finalement, il n'est pas simple d'analyser une source ponctuelle dont les photons arrivent 
generalement discontinus dans le temps et avec un angle d'incidence different. Ainsi, l'analyse 
suppose de connaitre sufiisamment bien les reponses de I'instrument (necessite de les valider, 
voir le chapitre |4|), ainsi que la liste des sources avoisinantes et remission diffuse du fond 



(voir 5.3.3). L'analyse se compose ainsi en trois etapes : (1) la selection des donnees en tenant 
compte de la PSF, (2) la creation d'un fichier contenant la liste des sources ponctuelles et 
diffuses de la region selectionnee, (3) I'ajustement des parametres spectraux par un maximum 
de vraisemblance. 

5.3.2 Selection des Donnees 

La premiere etape pour l'analyse spectrale d'un pulsar est la selection des donnees. La 
PSF etant large a basse energie (e.g. 68% des photons d'une source ponctuelle sont contenus 



dans une region de 4.9 degres a 100 MeV centree sur la source, voir la figure 4.2), il est 
indispensable de prendre une region (ROI) sufiisamment grande pour analyser I'ensemble des 
photons provenant de la source et ajuster correctement le fond diffus (environ 15° a I'interieur 
du plan galactique et 10° a I'exterieur). La contrepartie est qu'a I'interieur de cette region, 
les sources voisines contribuent a la liste de photons. II faut done inclure ces contributions 
dans l'analyse. La position et le spectre des sources ainsi que le fond diffus proviennent d'un 



catalogue (voir 5.3.3). Notons que la selection des donnees stockees dans un fichier fits (FTl) 
selon un intervalle de temps, un intervalle en energie, ou une region spatiale, est effectuee 
par le ST gtselect. De plus, specifiquement a I'etude des sources ponctuelles galactiques, 
seuls sont selectionnes les evenements de la classe diffuse qui ont la plus grande probabilite 
d'etre un photon. Enfin, les evenements possedant un angle zenithal superieur a 105° sont 
exclus de l'analyse. Au-dela de cette limite, une grande partie des photons y observes par le 
LAT resulte de I'interaction des rayons cosmiques avec I'atmosphere de la Terre produisant 
un fond brillant intense. Signalons que ces deux dernieres coupures s'appliquent egalement a 
l'analyse temporelle, mais nous tacherons de le rappeler dans les sections concernees. On pent 
egalement appliquer a cette etape une selection en phase des donnees. Cela permet d'estimer 
les parametres spectraux du pulsar en fonction de la phase, pour eventuellement identifier 
les processus physiques (i.e. rayonnement de courbure, synchrotron, et inverse Compton) qui 
correspondent aux differentes parties d'une courbe de lumiere. 

6. Rappelons qu'une orbite dure environ 95 minutes. 
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5.3.3 Modele des Sources 



La deuxieme etape dans I'analyse spectrale d'une source est la creation d'un fichier conte- 
nant le modele des sources contribuant a la ROI et au-dela. Ce modele code en format XMLQ 
pour etre utilise par les Science Tools inclut la position spatiale de la source a analyser, la 
position des sources avoisinantes, la forme fonctionnelle des spectres des sources, et la valeur 
des parametres spectraux. Toutes ces informations sont extraites des catalogues de sources 



generes par la collaboration Fermi (par exemple Abdo et al 2009-bslj). Rappelons brieve- 
ment que la figure [LS] presente la position dans le ciel des 271 sources du troisieme catalogue 
d'EGRET ( Hartman et al[ 1999). L'analyse spectrale des sources etudiees dans cette these 
utilise le catalogue du LAT (6 mois de donnees), lequel comptabilise au total 778 sources 
ponctuelles. Le catalogue base sur 11 mois contient plus de 1000 sources! 

Le modele inclut egalement les modeles d'emission diffuse galactique et extragalactique qui 
represente environ 90% des evenements du ciel y. L'emission diffuse galactique est modelisee a 
partir de GALPROPQqui est un modele numerique de propagation des rayons cosmiques dans 
la Galaxie, et est elabore pour calculer remission y ( [Strong et al] |2004a||bl ). Ce modele utilise 
les distributions de gaz de la Galaxie basees sur des observations spectrales de HI (21 cm) et 
CO (115 GHz), ainsi qu'une modelisation detaillee de la distribution stellaire galactique et 
de la distribution de la poussiere dans le milieu interstellaire. En resume, remission diffuse y 
provient de la desintegration des pions produits par I'interaction des rayons cosmiques avec 
le gaz interstellaire, de I'interaction des rayons cosmiques avec une population de photons 



de basse energie produits par les etoiles (processus inverse Compton, Porter et al 2008), 



ainsi que du rayonnement de freinage et la diffusion Compton inverse des electrons. Notons 
que ce modele est mis a jour regulierement. Ainsi, la version de GALPROP utilisee pour 
chaque analyse spectrale des pulsars etudies dans cette these sera indiquee dans des sections 
appropriees. L'emission diffuse extragalactique isotrope est, quant a elle, modelisee par une 
simple loi de puissance d'indice spectrale fixe a -2.25. Cela inclut egalement le fond residuel 
instrumental. II est possible que ce fond trouve son origine dans diverses classes d'objets 
extragalactiques non resolus comme les noyaux actifs de galaxie. 

En pratique, seuls les parametres spectraux de la source etudiee, des sources avoisinantes 
a moins de 10° de la source, ainsi que de l'emission diffuse sont laisses libres pendant I'ajuste- 
ment des parametres de la source. Cela permet d'avoir une limitation des parametres libres. 



Pour l'analyse des pulsars, tous les parametres de I'equation 1.23 peuvent etre ajustes. 



5.3.4 Maximum de Vraisemblance 



L'analyse des donnees du LAT par une methode de maximum de vraisemblance (Mattox 



et al 1996) est utilisee pour determiner les parametres spectraux de notre source. Cette 



methode est implantee dans les ST sous le nom de gtlike. Pour la decrire, prenons le cas 
simplifie d'une source ponctuelle, sans fond diffus, dont les photons sont observes avec un 
angle d'incidence 9 par le LAT. La source a un spectre en energie S(E',9') qui ne varie pas 
au cours du temps. La premiere etape est de creer une carte d'exposition. Les fonctions de 



reponse de I'instrument (voir 4.4) sont une fonction de Tangle d'incidence 6 et de I'energie E 



7. XML : extensible Markup Language. 

8. http ://galprop. stanford.edu/web_galprop/galprop_honie. html 
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mesures. Le nombre de photons observes d'une source depend done de la quantite de temps 
passe a observer la source pour un angle d'incidence donne. Ce temps depend seulement 
de I'orientation du LAT pendant la periode d'observation et non du modele de la source. 
Habituellement, la carte d'exposition est essentiellement I'integration de la surface efiicace 
qui depend de E et 9 sur le temps. Dans le cas du LAT, la carte d'exposition e(E',9') est 



obtenue en integrant les reponses de I'instrument (voir I'equation 4.1) dans I'espace de phase 
des quantites mesurees : 



efE',e' 



dE de dt R(E,e;E',e',t) , (5.9) 



SR 



ou E et 9 sont respectivement I'energie et Tangle d'incidence mesures. La "Region des sources" 
SR est une partie du ciel contenant toutes les sources qui contribue de maniere significative 
a la ROI. Ceci est necessaire pour s'assurer que les photons des sources a I'exterieur de la 
ROI sont pris en compte, en raison de la large PSF du LAT notamment a basse energie. 
A partir de cette fonction d'exposition, le nombre predit d'evenements observes pour une 
source donnee est : 



Npred 



dE' d9' S(E',9') e(E',9') (5.10) 

ROI 



En pratique, etant donne que le meilleur modele a la plus grande probabilite de s'ajuster 
aux donnees, les parametres de S(E',9') sont ajustes pour maximiser le logarithme de la 
fonction de vraisemblance considerant une densite de probabilite discrete (loi de Poisson) des 
evenements 1 : 

ln(L) =^ln(MO-Np,ed , (5.11) 

i 

OU Mi^ est la fonction de distribution des evenements a la coordonnee de I'evenement (i), soit : 

M,(Ei,9i,tO = [ dE' d9' R(E^,9i,tr,E(,9()S(E(,9() (5.12) 

J ROI 

La quantite IrL(L) est, en d'autres termes, une fonction non-lineaire des parametres spectraux 
de la source. On suppose que la vraisemblance L est derivable et que L admet un maximum 
global pour une valeur precise des parametres spectraux. Le maximum de ln(L) s'obtient en 
recherchant pour quelle valeur de ces parametres la derivee de L s'annule. On obtient en sortie 
de I'analyse la valeur avec I'incertitude des parametres spectraux laisses libres dans le modele 
des sources, le nombre de photons predits provenant de la source, ainsi que sa significativite. 
Notons que le ST gtlike propose plusieurs algorithmes pour maximiser le logarithme de la 
fonction de vraisemblance. lis se differencient par leur vitesse d'execution et leur precision. 



92 



P ARTIE III 



RESULT ATS 



93 



94 



Chapitre 6 



PSR J0205+6449 
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6.1 Introduction 



Ce chapitre rapporte la decouverte de pulsations provenant du pulsar PSR J0205+6449 
dans la nebuleuse a vent de pulsar 3C 58 avec le telescope Fermi Gamma-Ray Space Teles- 



cope. Dans la section 6.3 nous decrivons les observations radio et y utilisees pour I'analyse, 



tandis que les sections 6.4. 1| et 6.4.2 incluent respectivement I'etude de la courbe de lumiere 
y fournissant I'alignement avec remission radio et X, et I'etude spectrale de la source. Dans 
la section |6.5[ nous discutons de la courbe de lumiere, de la luminosite, de I'efiicacite y, de 



la geometrie de remission, et de I'age du pulsar. Ce chapitre est tire de I'article Abdoetal 



(2009-J0205) paru dans ApJL en 2009. 
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6.2 Pulsations dans la Nebuleuse a Vent de Pulsar 3C 58 



La source radio 3C 58 a ete reconnue en 1970 comme etant un reste de supernova 
(SNR G130.7+3.1; |Caswell||l970 '), et plus t ard classifiee comme une nebuleuse a vent de 
pulsar (PWN, egalement appelee plerion) par Weiler and Panagia (1978). Becker et al p982) 
identifierent une source ponctuelle dans le coeur de 3C 58 s'apparentant a un pulsar. La 



distance de la PWN est estimee a 3.2 kpc ( [Roberts et alj |1993[ ). Le pulsar PSR J0205+6449 a 
finalement ete decouvert dans les donnees du Chandra X-ray Observatory avec une periode 
de 65.7ms, tandis que les archives du telescope Rossi X-ray Timing Exp lorer ont permis 
une mesure de la derivee de la periode P = 1.93 x 10^^^ (Murray et al, 2002). Cela a ete suivi 



par la decou verte d'un tres faib le signal pulse en radio, dont la densite de flux moyen est 
de 



45 [iJj ( Camilo et al 



2002). Le pulsar a un taux d'energie de rotation perdue de E de 
2.7 X 1 0^^ erg s^^ qui en fait le troisieme pulsar le plus energetique de la galaxie, un champ 
magnetique de 3.6 x 10^^ G, et un age caracteristique t de 5400 ans. On se reportera a la 
section |1.2| pour le calcul de ces proprietes. II presente egalement un haut niveau d'instabi- 
lites temporelles (Ransom et al 2004). Deux glitchs ont ete observes depuis sa decouverte 



(Livingstone et al, 2008). On trouvera dans le tableau 9.1 les principales proprietes du pulsar. 

L'association possible du systeme 30 58/J0205+6449 et de la supernova historique SN 
1181 (828 ans) est un veritable debat. Si la PWN 30 58 coincide spatialement et est energe- 



tiquement compatible avec la supernova historique SN 1181 ^Stephenson 1971 Stephenson 



and Green 2002), de recents travaux sur les modeles de PWN (Ohevalier 2005), sur la vitesse 



d'expansion de la nebuleuse radio, et la vitesse des noeuds optiques impliquent un age pour 
30 58 de plusieurs milliers d'annees, plus proche de I'age caracteristique du pulsar que de SN 
1181. La nebuleuse 30 58 est egalement tres similaire a la nebuleuse du Orabe agee de 955 
ans (reste de supernova SN 1054), composee d'une nebuleuse radio avec un spectre plat, une 
emission X etendue non thermique, et une source ponctuelle X due au pulsar central. Malgre 
cela, les deux objets different significativement en luminosite et en taille. La nebuleuse radio 
30 58, bien que deux fois plus grande, est moins lumineuse par un ordre de grandeur que 



la nebuleuse du Orabe (Ivanov et al, 2004), tandis que sa luminosite en X est 2000 fois plus 



petite (Torii et al 2000). Oes disparites pourraient etre expliquees par un age different. 



6.3 Observations 

6.3.1 Observations Radio 

Le comportement rotationnel de PSR J0205+6449 durant les 8 premiers mois de la mission 



Fermi est decrit par la solution temporelle (ephemeride) presentee dans le tableau 6.1 Le 
pulsar a ete observe a la fois par le NRAO Green Bank Telescope (GBT) et le Lovell telescope 
a Jodrell Bank. GBT a founi des donnees plus precises, tandis que des observations ont ete 
realisees plus frequemment a Jodrell Bank. L'ephemeride utilisee pour empiler les photons y 



est basee sur les TOAs (voir la section 5.2.2) obtenus par les deux telescopes entre le 17 juin 
2008 et le 9 mars 2009. Dix-sept TOAs proviennent du radiotelescope GBT a la frequence de 
2.0 GHz avec une incertitude moyenne de 0.27ms, chacun est base sur 1 heure d'integration. 
Quant aux 51 TOAs de Jodrell Bank dont I'incertitude moyenne est de 0.37 ms, ils derivent 
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d'observations a 1.4 GHz basees sur 2 heures d'integration. Notons que la chronometrie de 
PSR J0205+6449 est tres bruitee. Dans le but de decrire correctement sa rotation pendant 
les 9 mois d'observation y, Tephemeride est composee des sept premieres derivees de la 
frequence. Ce travail a ete effectue avec le programme d'analyse TEMPOj^ Un ecart-type 
(RMS) de 0.4 ms a ete obtenue. La meilleure determination de la mesure de dispersion (DM 



140.7 ± O.Spccm ^) est citee dans Camilo et al (2002), laquelle est utilisee pour corriger 



les TOAs a une frequence infinie, avec une incertitude de 0.3 ms. 



Parametres 


Valeurs 


Ascension Droite ( J2000) 


02^05^37.92^ 


Declinaison (J2000) 


+64°49'42.8000014" 


Frequence FO (s^^) 


15.2144837749(6) 


^1 (s j 


-4.4ozUU^^o J X iU 


F2 (s-3) 


-1.43(2) xlO-^° 


F3 (s-4) 


-2.7(5) X 10-28 


F4 (s-5) 


6.22(7) X 10-33 


F5 (s-^) 


2.05(6) X 10-33 


F6 (s-^ 


-1.43(1) xlO-4^ 


F7 (s-«) 


-7.3(1) X 10-^2 


Epoque, To (MJD[TDB]) 


54767.0 


Date de reference pour le cf)o (MJD[TDB])'' 


54768.07703065005798 


Intervalle de validite en MJD 


54634.253 - 54899.693 


Nombre de TOAs 


68 


Mesure de dispersion, DM(cm^^ pc) 


140.7 ± 0.3 


Modele d'ephemeride du S.S 


DE405 



Tableau 6.1 - Parametres rotationnels du pulsar PSR J0205+6449 bases sur les observations 
des telescopes Green Bank Telescope (GBT) et Jodrell Bank. Les chiffres entre parentheses 
representent les incertitudes donnees par TEMPO. 
" Date a laquelle la phase cf) = (TZRMJD). 



6.3.2 Observations Gamma 



Les donnees coUectees pour ce travail ont ete obtenues par le Fermi LAT dans deux 
modes d'observations differents, du 30 juin au 29 juillet 2008 durant la phase de verification 
des instruments, et du 3 aout 2008 au 9 mars 2009 quand le LAT operait en mode de balayage. 
Nous avons utilise les evenements diffuse ayant la probabilite la plus elevee d'etre un photon 
et exclu ceux possedant un angle zenithal superieur a 105° a cause de I'albedo terrestre. 
La source qui a une significativite statistique de 9.5 o (< 10 a) de 0.2 a 100 GeV pour les 
trois premiers mois de la mission (mode de balayage) n'apparait pas dans le catalogue des 



sources brillantes ( |Abdo et al 2009-bsl). La figure 6.1 montre la distribution des photons de 



1. TEMPO est I'ancienne version de TEMP02 : http ://www.atnf.csiro.au/research/pulsar/tenipo 
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la region entourant le pulsar. Notons la presence du micro-quasar LSI-|-61°303 (Abdo et al 



2009-LSI61), et de deux autres pulsars emetteurs y, PSR J0248+6021 (Cognard et al 2009) 



et PSR J0007+7303 le pulsar du reste de supernova CTAl ( |Abdo et al[ [2008j). 
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Figure 6.1 - Distribution des photons de la region qui entoure le pulsar PSR J0205+6449- 
Les croix representent les sources du catalogue pour 6 mois de donnees. Le nombre de degre 
par pixel est de 0.25. 
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6.4 Analyse et Resultats 



6.4.1 Courbe de Lumiere 

Le pulsar est localise dans le plan galactique ou remission diffuse y est intense, et a 5.3° de 



la source brillante binaire LSI+6r303 (Abdoetal, 2009-LSI61). Pour I'analyse temporelle, 



un ensemble de photons avec une energie superieure a 100 MeV a ete selectionne dans un cone 
dependant en energie de rayon Ogg ^ 0.8 x Eq^v^ degres, mais avec un rayon maximum de 
1.5° centre sur la position X du pulsar (I = 130.719°, b = 3.085°). Ce choix prend en compte 
les capacites de I'instrument et maximise le rapport signal-sur-bruit sur un large intervalle 
en energie. Cela tronque la PSF a basse energie et reduit le nombre d'evenements du fond. 
Un total de 2922 photons restent apres ces coupures. A partir du programme d'analyse 
TEMP02, les temps d'arrivee des evenements ont ete transferes au barycentre du systeme 
solaire en utilisant I'ephemeride du systeme solaire JPL DE405, et les evenements ont ete 
empiles suivant I'equation |5.2| en utilisant les parametres de I'ephemeride radio decrit dans 
le tableau 16. 1[ 



La figure 6.2 (cadre du haut) montre le phasogramme y au-dessus de 100 MeV, sur lequel 
nous avons estime la position et la largeur des pics. Le premier pic (PI) est decale par rapport 
au pic radio (cadre du bas) par 0.08 ± 0.01 ± 0.01 selon I'ajustement d'une lorentzienne 
de largeur a mi-hauteur (FWHM) de 0.15 ± 0.01. Pour la position, la premiere incertitude 
provient de I'ajustement du pic, et la seconde provient de I'incertitude de la mesure de 



dispersion en extrapolant les TOAs radio vers une frequence infinie (voir I'equation 2.1). Le 



second pic est asymetrique. II a done ete evalue avec deux demi-lorentziennes pour prendre en 
compte les pentes differentes a gauche et a droite. L'ajustement place le pic a 0.57±0.01 ±0.01 , 
avec une largeur a mi-hauteur de 0.13 ± 0.04. Les deux pics sont separes par 0.49 ± 0.01 en 
phase. Nous avons defini la region non pulsee (OFF) comme etant le minimum pulse entre 0.65 
et 1.0 en phase. La ligne pointillee represente le fond estime a partir d'un anneau entre 2° - 3° 
autour du pulsar. Un exces entre PI et P2 (cf) = 0.14 — 0.46) apparait avec une significativite 
( signal/ V^fond) de 5 a. Cependant, les donnees ne peuvent contraindre la dependance en 
phase de cet exces. 

Pour examiner le profil y en fonction de I'energie, des phasogrammes avec 50 intervalles 
en phase ont ete traces pour trois rangs en energie, 0.1 - 0.3 GeV, 0.3 - 1 GeV, et >1 GeV 
(cadres du milieu). Entre 0.1 et 1 GeV, deux pics et un possible epaulement apres le premier 
pic sont distinctement presents, tandis qu'au-dessus d'lGeV seul le pic P2 est significatif. 
Nous observons egalement que le rapport P1/P2 (somme des coups) diminue quand I'energie 
augmente, avec un rapport de 0.63 ± 0.06 entre 0.1 - 0.3 GeV, 0.55 ± 0.05 entre 0.3 - 1.0 



et 0.24 ± 0.06 au-dessus de 1 GeV. Les pulsars de Vela ( |Abdo et all |2009-Vela[ ), du Crabe 
(Abdo et al 2009-Crab), de Geminga, et de PSR B1951-I-32 montrent le meme comportement 



(Thompson 2001). Cette tendance generale suggere une dependance spectrale en energie de 



la courbe de lumiere y. Finalement, nous notons que le photon de plus haute energie se situe 
dans le pic P2 avec une energie de 8.6 GeV. 



La figure 6.2 (derniers cadres) montre egalement la courbe de lumiere X entre 2 et 60 



keV mesuree par le telescope RXTE (Livingstone et al, 2009), ainsi que le profil radio a 2 GHz 
provenant du telescope GBT. Le pic radio definit la phase de reference cj) = 0. Livingstone 
et al (2009) montrent que le pic radio devance le premier pic X par 4) = 0.10 ±0.01, et notent 
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Figure 6.2 - Cadre du haut : Courbe de lumiere de PSR J0205+6449 au-dessus de 0.1 GeV. 
Deux rotations sont representees avec 50 intervalles en phase par periode (V = 65.7 ms). La 
ligne pointillee montre le niveau de fond estime a partir d'un anneau autour du pulsar (46 
coups/intervalle). Trois cadres suivants : Phasogrammes pour trois intervalles differents 
en energie, chacun est represents avec 50 intervalles en phase par periode. Second cadre a 
partir du bas : Nombre de coups dans la bande en energie 2—60keV a partir des donnees 
RXTE (Livingstone et al 2009). Cadre du bas : Profil pulse radio base sur 3.8 heures 



d' observations du telescope GET a la frequence de 2 GHz avec 64 intervalles en phase. Cette 



figure est extraite de ^Abdo et al (2009-J0205) . 
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6.4. Analyse et Resultats 



une separation entre les deux pics etroits de Acj) = 0.5. Le bon alignement entre les profils 
X et y suggere une origine commune entre les deux composantes. Cette caracteristique est 
egalement observee dans les donnees LAT pour le pulsar de Vela, ou le plus intense pic X est 



aligne avec le premier pic y (Abdo et al 2009-Vela). 



6.4.2 Spectre et Flux Resolus en Phase 

Pour etudier le spectre resolu en phase de PSR J0205+6449, une analyse spectrale par 
un maximum de vraisemblance a ete effectuee en utilisant une carte y (ROI) de 20° centree 
sur la position du pulsar entre 0.1 et 200 GeV. Rappelons que les erreurs systematiques sur 
la surface efficace estimees pour les IRFs "Pass6_v3" sont < 5% autour de 1 GeV, 10% en 
dessous de 0.1 GeV et 20% au-dessus de 10 GeV. L'emission diffuse pour le plan galactique a 
ete modelisee en utilisant une carte basee sur le modele GALPROP : version gll_iem_v01. 
Quant a remission extragalactique et au bruit de fond instrumental, ils ont ete modelises par 
une simple loi de puissance. Les sources ponctuelles a I'interieur de la ROI et au-dela (25°) 



ont ete incluses (voir la figure 6.1). Pour le pulsar, la forme du spectre a ete modelisee par 



une loi de puissance avec une coupure exponentielle (voir I'equation 1.23). 

Premierement, nous avons determine les composantes spectrales du fond diffus en selec- 
tionnant la region OFF du pulsar. Ensuite, nous avons ajuste les donnees de la region ON 
(cj) = 0.00 — 0.65) du pulsar pour ameliorer le rapport signal-sur-bruit. Le meilleur ajustement 
est decrit par : 

= No E"^ g-E/Ecutoff cm"2g-lQgy-l /g 

dE 

avec E en GeV, le terme Nq = (1 .4±0.1 ±0.1 ) x 1 0"*^ cm"^ s^^ GeV^^ , I'indice spectral T = 2.1 ± 
0.1 ± 0.2, et I'energie de coupure Ecutoff = ^-O-o 7 ^ ^-^ GeV. Les erreurs sont respectivement 
les incertitudes statistiques et systematiques. Nous obtenons de cet ajustement entre 0.1 - 200 
GeV un flux de photon integre de [13.7±1 .4±3.0) x 10^^ cm^^ s^^ et un flux energetique integre 
de pE.obs = (6.7 ± 0.5 ± 1.0) x 10^^^ ergscm^^s^\ Les erreurs systematiques ont ete obtenues 



avec les IRFs modifiees (voir la section 4.5.4). La figure 6.3 montre le spectre d'energie multi- 
longueurs d'onde du pulsar. Les points Fermi ont ete obtenus en effectuant une analyse 
spectrale dans chaque intervalle en energie en supposant une simple loi de puissance. 

Pour verifier I'hypothese d'une coupure dans le spectre, nous avons ajuste le meme lot 
de donnees avec une simple loi de puissance de la forme dE/dE = lMo(E/1 GeV)^'^. Le modele 
spectral utilisant la coupure est mieux contraint avec une difference dans les logarithmes du 
maximum de vraisemblance de ~ 4.5 0, defavorisant I'hypothese d'une simple loi de puissance. 
Nous avons egalement ajuste les donnees avec un indice d'attenuation |3 librej^ L'ajustement 
donne (3 = 0.8 ± 0.2 ± 0.3, ce qui exclut le modele avec une coupure super-exponentielle. 

Une recherche d'emission dans la region OFF du pulsar correspondant a une possible 
emission de la PWN 3C 58 a ete effectuee. Nous avons ajuste les donnees de la region OFF 
par une loi de puissance et en supposant la nebuleuse ponctuelle. Aucun signal n'a ete observe 
de la PWN. Finalement apres une renormalisation du signal sur toute la phase, nous avons 
derive une limite superieure a 95% de confiance sur le flux integre au-dessus de 200 MeV de 
1.7 X lO-^cm-^s-i. 

2. Rappelons que le modele Polar Cap predit un spectre avec une coupure en energie de la forme Exp[- 
(E/Ec)^], avec (3-2. 
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Figure 6.3 - Spectre en energie multi-longueurs d'onde pour le pulsar PSR J0205+6449- 



Shearer and Neustroev ( 2008 ) ; Chandra : Murray et al 



Les references pour cette figure sont : radio : Camilo et al ( 2002 ) ; optique ( WH Telescope ) : 



(donnees non publies) ; Fermi : cette these; Veritas : (Aliu 2008^ 



(2002); RXTE : Kuiper et Hermsen 



6.5 Discussion 



6.5.1 Courbe de Lumiere 



Les courbes de lumieres multi-longueurs d'onde sont essentielles pour localiser remission 
pulsee dans les regions des lignes de champs ouvertes et, par consequent, pour comprendre 
les mecanismes d'acceleration de particule. Le profil y de PSR J0205-I-6449 couvre un large 



intervalle de phase et est similaire a la courbe de lumiere de Vela (Abdo et al, 2009-Vela), 
incluant un alignement des pics y avec les pics X. Le decalage 5 radio-y et la separation 
A des deux pics y par ~ 0.5 est en train de devenir un modele preponderant (lAbdo et 



2009-psrcat). Les pulsars de Vela, B1951+32, et J2021+3651 (Halpern et al 



|2008j .Abdol 



\bdo| 



et al 2009-J2021) ont les memes caracteristiques, un premier pic y a 0.13, 0.16, et 0.17 



respectivement par rapport au pic radio, et une separation des pics y compris entre 0.4 et 
0.5. Cela s'ajuste assez bien aux predictions des modeles magnetospheriques externes, que 
ce soit le modele traditionnel Outer Gap (OG) ou le modele plus recent Two-Pole Caustic 
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6.5. Discussion 



(TPC), et exclut done le modele Polar Cap. Anticipons ! La figure 9.4 presente la distribution 
5 — A des pulsars y detectes par le LAT pour 6 mois de donnees. 



6.5.2 Luminosite, Efficacite, et Geometrie de rEmission 

Pour connaitre I'efficacite y, il est essentiel de determiner au prealable la luminosite totale 
Ly (voir I'equation Cependant, une incertitude importante pour evaluer la luminosite 

est la determination de la distance D. Des observations de I'absorption de I'hydrogene neutre 
(H I) par Roberts et al| ( [l993| ) donnent une estimation de la vitesse radiale de la PWN 3C 58. 
Cette vitesse est convertie en une distance cinematique de 3.2 kpc en supposant une courbure 



rotationnelle galactique plate (Fich et al 1989). Les incertitudes sur la distance avec cette 



methode sont de I'ordre de 25%, rendant ce resultat compatible avec la distance de 2.6 kpc 



precedemment reportee par Green and Gull (1982|. La distance derivee du DM est quant a 
elle de 4.5 kpc selon le modele NE 2001, mais avec une incertitude qui pent exceder 50%. 
Pour ce travail de these, nous adopterons une distance de 3.2 kpc etant donne I'incertitude 
de cette valeur. 

A partir du flux energetique integre, nous estimons une luminosite : 



= 8.3 X 10^4 (D/3.2kpc)^ fa ergs 



-1 



(6.2) 



ou fo est le facteur de correction (voir I'equation 1.12), et deduisons pour un moment d'inertie 



de lO^^gcm^ une efficacite ^ = L^/E = 0.003 fo (D/3.2kpc)^ 

A present, il est necessaire d'evaluer le facteur de correction qui prend en compte la 
geometrie de remission. En supposant les modeles magnetospheriques externes suite a la 
discussion sur la forme de la courbe de lumiere et I'exclusion du facteur d'attenuation |3 de 
la forme spectrale, nous pouvons deduire I'inclinaison magnetique a. et Tangle d'observation 
C a partir de la separation des pics y et en utilisant les courbes de lumieres du papier "Atlas" 
de |Watters et al| ( |2009] ). Pour les modeles OG, nous estimons fo ^ 0.9-1 .0 avec a - 60°-90° 

0.95-1.25 avec a la fois le couple 
45°-90°. De plus, I'examen de la 



et C 80°-85°, tandis que pour les modeles TPC, f^ 

a - 50°-90°, C ^ 85°-90° et le couple cc - 85°-90°, C 

geometrie de la PWN 3C 58 utilisant les donnees Chandra ( |Ng and Romani[ |2004[ |2008D 
mene a un angle d'observation de C = 91 .6 ± 0.2 ± 2.5 ° (tore interne) base sur Tangle 
d'inclinaison du tore au plan du ciel. Cette valeur est en accord avec les resultats du modele 
OG et plus conforme avec le premier intervalle estime pour le modele TPC. Si on adopte une 
valeur moyenne de fn = 1 (compatible avec les deux modeles) et une distance de 3.2 kpc, 
nous en deduisons une efficacite r\ = 0.3% pour convertir sa perte d'energie rotationnelle en 
photons y. Cette valeur de I'efficacite montre que PSR J0205+6449 suit relativement bien la 
relation approximee par r\ cc E^^^^ a 1 /V (voir la figure 



9.2). 



6.5.3 Age du Pulsar 



Thompson et al ( |1999 ) soulignent que pour les pulsars EGRET, le maximum de lumino- 



site pour les jeunes pulsars y tels que le Crabe ou PSR B1509— 58 est au niveau du rayon- 
nement X, tandis que pour les pulsars moins jeunes tels que Vela ou Geminga le maximum 



de luminosite est a haute energie, dans Tintervalle d'energie du LAT (voir la figure 1.5). Le 
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tableau 6^ presente, pour les jeunes pulsars y connus|^ la luminosite X et y, ainsi que le 
rapport Lx/Ly et I'age caracteristique. Considerant une luminosite X entre 0.5 et 8keV de 
1.51 X 10^^ ergs s^^ ( jKargaltsev and Pavlov 2009) et une luminosite Ly de 8.3 x 10^^ ergs s^^ 
avec fo = 1 , PSR J0205+6449 semble avoir un rapport Lx/Ly plus proche de celui des 
pulsars tels que Vela et PSR J2229+6114, que des tres jeunes pulsars tels que le Crabe et 
PSR B1509— 58. Par consequent, cette etude suggere que I'association de la PWN 3C 58 qui 
est alimentee par le pulsar, avec la supernova historique SN 1181, est probablement incorrecte. 



Pulsar 


Luminosite X'' 


Luminosite y 


Lx/Ly 






Refs. 




Lx (ergs-^) 


Ly (ergs-i) 




(an) 




Crab 


1.55 X 10^^ 


6.2 X 10^^ 


2.5 




1300 


1,2 


B1509-58^ . 


1.09 X 10^5 


1.6 X 10^5 


0.7 




1500 


1,2 


J0205+6449 


1.51 X 10^3 


6.1 X 10^4 


2.5 X 10" 


-2 


5400 


1,6 


J2229+6114 


6.76 X 10^2 


2.6 X 10^5 


2.6 X 10- 


-3 


10500 


1,7 


Vela 


3.02 X 10^^ 


8.7 X 10^4 


3.5 X 10- 


-4 


11000 


1,4 


J2021+3651^ 


9.19 X 10^2 


8.2 X 10^5 


1.1 X 10- 


-4 


17000 


3,5 


B 1706-44 . . . 


1.10 X 10^2 


2.9 X 10^5 


3.8 X 10- 


-4 


17000 


1,2 


J1048-5832 . 


2.6 X 10^^ 


2.1 X 10^^ 


1.2 X 10- 


-4 


20400 


1,7 



Tableau 6.2 - Luminosites Lx et Ly des jeunes pulsars connus. x — P/2P est I'age 
caracteristique des pulsars. La luminosite est calculee en supposant le facteur de correction 
de remission fn = 1, o, I'exception de (h). 

References : 

1. JKargaltsev and Pa'uZot^ |gOQgp ; 2. ^Abdo et al\\2009-psrcaU^ ; 3. ([Van Etten et al\ \2008^ ; 
4. \Abdo et aA \2009-VeM ; 5. (Abdo et al\ \2009-J2021\ ) ; 6. jAbdo et al\ \2009-J0205 ) ; 



7. (Abdo et al 



2009-J2229) 



Notes : 

°- Luminosite entre 0.5 et 8keV. 

^ Luminosite totale (optique et au-dessus) calculee en supposant fo 
Luminosite calculee en supposant une distance de 4 kpc. 



l/4n. 



3. Nombre de pulsars y connus au moment de I'etude. 
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7.1 Introduction 



Ce chapitre rapporte la decouverte de pulsations provenant du pulsar PSR J2229+6114 
avec le telescope Fermi Gamma-Ray Space Telescope. Ce pulsar est sur beaucoup d'aspects 
similaire au pulsar de Vela[^ tant par ses emissions pulsees radio, X, et y, que par la presence 
d'une PWN autour du pulsar et une instabilite rotationnelle importante. Kramer et al (2003) 



fournissent une discussion sur les pulsars de type Vela. Le pulsar PSR J1048— 5832 au chapitre 
suivant montre les memes caracteristiques. Dans la section [773 | nous decrivons les observations 
radio et y utilisees pour I'analyse, tandis que les sections [7. 4 . 1 1 et [7. 4 . 2 1 incluent respectivement 
I'etude de la courbe de lumiere y fournissant I'alignement avec remission radio et X, et I'etude 
spectrale de la source. Dans la section |7.5[ nous discutons de la courbe de lumiere, de la 
luminosite, de I'efiicacite y, et de la geometrie de remission. Ce chapitre est tire de I'article 



Abdo et al (2009-J2229), a paraitre dans ApJ en 2009. 



1. Ce pulsar est classe comme un pulsar "vela- like". 
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7.2 Candidat pour la Source EGRET 3EG J2227+6122 



PSR J2229+6114 se trouve aux coordonnees galactiques (l,b)=(106°,2.9°) a I'mterieur de 
la boite d'erreur de la source EGRET 3EG J2227+6122 ( |Hartman et al[ [1999^ . L'objet a 
premierement ete detecte dans les donnees X de ROSAT et ASCA comme etant une source 
compacte, et ensuite identifie comme etant un pulsar grace a la decouverte de pulsations 



radio et X de periode 51.6ms (Halpern et al, 2001b). Le profil radio montre un simple pic 
etroit, tandis que la courbe de lumiere X entre 0.8 et lOkeV se compose de deux pics, separes 
par Acf) = 0.5. La collaboration A GILE ^ a re cemment annonce la decouverte de pulsations y 
au-dessus de 100 MeV (Pellizzoni et al 2009). Le pulsar est aussi jeune que le pulsar de Vela 
(age caracteristique t = P/2P = 10000 ans), et aussi energetique (E = 2.2 x lO^^ergs^^). II est 
probablement la source energetique de la PWN "Boomerang" G106.65+2.96, qui appartient au 



reste de supernova G106.3+2.7 decouvert par Joncas and Higgs (1990). Notons que recemment 



la PWN a ete detectee au TeV par le telescope MILAGRO ( |Abdo et al^ |2009-mil| ). Des 
etudes sur la vitesse radiale de I'hydrogene neutre et de I'environnement moleculaire placent 
le systeme a ~800pc (Kothes et al 2001), tandis que Halpern et al (2001a) suggerent une 



distance de 3kpc estimee a partir de I'absorption X. La distance determinee par le DM en 
utilisant le modele de la distribution d'electron dans la galaxie NE 2001 donne une distance de 
7.5 kpc, significativement au-dessus des autres estimations. Etant donne que les incertitudes 
sur cette distance sont grandes (jusqu'a 50%), nous adopterons pour ce travail de these une 
distance de 3 kpc. 



7.3 Observations 

7.3.1 Observations Radio 

PSR J2229-I-6114 est observe par le telescope GBT et le telescope Lovell a Jodrell Bank. 
L'ephemeride utilisee ici pour le calcul des phases a partir des temps d'arrivee des photons 
y est basee sur des TOAs obtenus par les deux telescopes entre le 17 juin 2008 et le 23 
mars 2009. Vingt-cinq TOAs proviennent du GBT avec une incertitude moyenne de 0.2 ms, 
chacun represente 5 minutes d'observation a la frequence centrale de 2.0 GHz, tandis que 
quarante-quatre TOAs avec une incertitude moyenne de 0.3 ms proviennent des observations 
de Jodrell Bank obtenus a partir de 30 minutes d'observations a la frequence de 1.4 GHz. 
Pour I'ajustement des parametres de rotation avec TEMPO, la position du pulsar a ete fixee 



a la position X connue (Halpern et al 2001b). PSR J2229-I-6114 montre du timing noise sur 
les 9 mois de donnees, et done demande un ajustement de la frequence et des deux premieres 
derivees. De plus, un petit glitch s'est produit le 12 novembre 2008 (54782.6 ±0.5 MJD) avec 
un pas de frequence fractionnaire AF/F = (4.08 ± 0.06) x 10"^ et AF/F = (2.0 ± 0.4) x 10"'*. 
Celui-ci a ete pris en compte dans le modele rotationnel. La rms de I'ajustement posterieur 
est de 0.24 ms. La DM = (204.97 ± 0.02) pccm^^ a ete mesuree a partir d'observations de 
GBT pour trois frequences differentes. Celle-ci est utilisee pour corriger le temps d'arrivee a 
2 GHz vers les frequences infinies, pour la comparaison avec le profil y. 

2. AGILE : Astro-rivelatore Gamma a Immagini LEggero 



106 



7.4. Analyse et Resultats 



7.3.2 Observations Gamma 

Les donnees utilisees pour I'analyse spectrale ont ete coUectees a partir du 4 aout jusqu'au 
23 mars 2009 pendant que le telescope etait en mode de balayage. Pour I'analyse temporelle, 
nous avons additionne les donnees de la phase de calibration de I'appareil, periode s'etendant 
du 25 juin au 3 aout 2008. Pendant cette phase, plusieurs configurations ont ete testees qui 
ont affectees la resolution en energie et la reconstruction des evenements du LAT. Cependant, 
ces changements n'ont pas d'effets sur la chronometrie du LAT. De plus, seuls les photons 
marques "diffuse" ont ete selectionnees, et en addition nous avons exclu ceux possedant un 



angle zenithal superieur a 105°. La figure 7.1 montre la distribution des photons de la region 



entourant le pulsar. Notons la presence du pulsar y PSR J2238+59 (voir le chapitre^ a ~2° 



de PSR J2229+6114 (Abdo et al 2009-BS) 



7.4 Analyse et Resultats 

7.4.1 Courbe de Lumiere 

Pour I'analyse temporelle de PSR J2229+6114, nous avons selectionne les photons d'ener- 
gie superieure a 100 MeV dans un rayon de 1 degre autour de la position du pulsar radio. 
Ensuite, nous avons selectionne les photons dans un cone dependant en energie de rayon Ogg ^ 
0.8 X Eq^y^ degres, gardant tous les photons inclus dans un rayon de 0.35°. Cette selection 
se base sur la PSF du LAT et permet d'augmenter le rapport signal-sur-bruit sur une large 
bande en energie. A partir du programme d'analyse TEMP02, nous avons corrige le temps 
d'arrivee des photons vers le barycentre du systeme solaire en utilisant I'ephemeride solaire 
JPL DE405. La phase des evenements y a ensuite ete calculee a partir de I'ephemeride radio 



delivree par les telescopes de GBT et de Jodrell Bank (voir le tableau 5.1) suivant I'equation 



La figure 7.2 (cadre du haut) montre I'histogramme resultant du nombre de coups en 
fonction de la phase au-dessus de 0.1 GeV. En comparaison, le cadre du bas montre le profil 
radio pulse aligne en phase provenant du telescope GBT. La courbe de lumiere y montre 
un seul et asymetrique pic = 0.15 — 0.65). Nous avons ajuste le pic avec deux demi- 
lorentziennes pour prendre en compte les deux pentes. L'ajustement place le pic a 0.49 ± 
0.01 ± 0.001 avec une largeur a mi-hauteur de 0.23 ± 0.03. Nous avons estime le niveau du 
fond a partir d'un anneau centre sur la position du pulsar entre 1 et 2° en selectionnant le 
minimum pulse de I'histogramme entre 0.65 et 0.15. Cela est represente par la ligne pointillee 
(48 coups/inter valle) et est homogene avec la region OFF. A partir de ces resultats, le nombre 
total de photons pulses est estime a 1365 ± 79, avec une contribution du fond de 2431 ± 49 
evenements. 



Pour examiner la forme du profil pulse en fonction de I'energie (voir la figure 7.2 cadres 
du milieu), quatre phasogrammes ont ete traces suivant differents intervalles (0.1 — 0.3 GeV, 
0.3 — 1 GeV, 1 — 3 GeV et >3 GeV), montrant un possible deplacement du pic. Entre 0.1 - 
0.3 GeV, le pic y est decale selon un ajustement avec deux demi-lorentziennes de 0.51 ± 0.02 
en phase par rapport au pic radio, tandis que les decalages pour 0.3 - 1 GeV, 1 - 3 GeV et 
> 3 GeV sont respectivement 0.48 ± 0.01 , 0.49 ± 0.01 , et 0.45 ± 0.01 . Notons que la position 
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Figure 7.1 - Distribution des photons de la region qui entoure le pulsar PSR 12229+6114- 
Les croix representent les sources du catalogue pour 6 mois de donnees. Le nombre de degre 
par pixel est de 0.25. 



du pic est relativement stable entre 0.1-3 GeV, mais il apparait un leger desalignement 
au-dessus de 3 GeV. Ainsi, les donnees ne suggere que legerement la dependance en fonction 
de I'energie, meme si le desalignement entre le profil X et y est tres prononce. Finalement, 
notons que le photon de plus haute energie est de 10.8 GeV et se situe a la phase ~ 0.30, ou 
Ton observe un faible exces (au-dessus de 3 GeV), qui pourrait correspondre a un second pic. 



La figure 7.2 presente egalement le profil X de PSR J2229-I-6114 obtenu a partir d'une 
observation de XMM entre 1 et lOkeV le 15 juin 2002 (52440 MJD) pour un temps d'ex- 
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position de 20 ks. Les donnees ont ete empilees en utilisant une ephemeride contemporaine 
basee sur des observations de Jodrell Bank et de GBT. Le pic le plus haut du profil X a un 
retard de cf) = 0.1 7 ±0.02 par rapport an pic radio. L'analyse d'une observation du telescope 
RXTE a partir de 52250 MJD donne un decalage similaire. II n'y a pas de dependance en 
energie de la forme pulsee X dans I'intervalle 1-10 keV, tandis que la finesse des pics aussi 
bien que leur forme spectrale indiquent que remission est majoritairement non-thermique. 



7.4.2 Spectre et Flux Resolus en Phase 

Initialement, une region circulaire de 15 ° autour de la position du pulsar et le fond diffus 
galactique base sur le modele numerique GALPROP appele 54_59Xvarli7Sgont ete utilises 
pour ajuster le pulsar. Cependant au vu des residus, certaines structures autour de I'objet 
n'ont pas ete prises en compte dans le modele galactique et pourraient surestimer ou sous- 
estimer le flux du pulsar. Finalement, nous avons adopte un modele different pour remission 
diffuse galactique basee a la fois sur six cartes "anneaux" galactocentriques de H I0et de CO, et 
sur la distribution spatiale de I'intensite Compton inverse modelisee par le Galprop 54_77. 
L'intensite aussi bien que les emissivites y par anneaux ont ete ajustees pour maximiser 
I'accord avec les observations, ceci en tenant compte des sources ponctuelles detectees en y. 
Cette approche est similaire a ce qui est decrit dans Casandjian and Grenier p008) pour 



modeliser les donnees EGRET. Ce modele est a present utilise par la collaboration LAT dans 
le catalogue 9 mois. 

PSR J2229+6114, reference comme la source ponctuelle du LAT OFGL J2229.0+6114 



dans le catalogue de 3 mois (Abdo et al 2009-bslj), a ete modelise par une loi de puissance 



avec une coupure exponentielle de la forme 

^ = No e-^/^^ cm-^s-^GeV-^ (7.1) 
avec E en GeV. L'ajustement a ete fait dans la region ON entre 0.15 et 0.65 en phase. La 



figure 7.3 (ligne continue) montre la distribution d'energie spectrale pour la region ON a 
partir d'un ajustement par un maximum de vraisemblance. Get ajustement donne un terme 
No = (5.2 ± 0.4 ±0.1) X 10-**cm-2s-i GeV-\ un indice spectral V = 1 .85 ± 0.06 ± 0.10 et 
une energie de coupure Ej. = 3.6lo g ± 0.6 GeV. A partir de ce resultat, nous avons estime un 
flux de photon integre de (3.77 ± 0.22 ± 0.44) x 10^^cm^^s^\ Cette valeur est 10% inferieure 
au flux de (4.13 ± 0.61) x lO^^cm^^s"^ de la source EGRET 3EG J2227+6122 obtenue a 
partir de la loi de puissance avec un indice spectral de 2.24 ±0.14, mais plus grande que le 



flux de (2.6 ± 0.4) x lO^^cm^^s"^ mesure par AGILE ( [Pellizzoni et al[ [2009j). Nous avons 
egalement ajuste les donnees avec une simple loi de puissance, menant a un flux de (4.54 ± 
0.14) X 10^'' cm^^ s^^ et a un indice spectral de 2.25 ±0.02, compatible avec l'analyse EGRET. 
Notons que le modele spectral avec une coupure exponentielle est mieux contraint avec une 
difference des logarithmes des maximums de vraisemblance de ~ 9ct, rejetant I'hypothese 
d'une simple loi de puissance. Un modele avec une coupure super-exponentielle de la forme 
Exp[-(E/Ec)'^] a ete egalement teste. L'ajustement donne (3 = 0.8 ± 0.4 ± 0.3, ce qui est en 
accord avec une simple coupure exponentielle |3 = 1 . 



3. Modele standard utilise pour le catalogue 6 mois. 

4. Les densites de colonnes de H I proviennent des observations du L.A.B, http ://www. astro. uni- 
bonn.de/ -webaiub/english/ tools _ labsurvey. php 
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Figure 7.2 - Cadre du haut : Courbe de lumiere de PSR J2229+6114 au-dessus de 0.1 GeV. 
Deux rotations sont representees avec 50 intervalles en phase par periode (V = 51.6 msj. La 
ligne pointillee montre le niveau de fond estime a partir d'un anneau autour du pulsar (48 
coups/intervalle). Quatre cadres suivants : Phasogrammes pour quatre intervalles diffe- 
rents en energie, chacun est represente avec 50 intervalles en phase par periode. Second 
cadre a partir du bas : Nombre de coups dans la bande en energie 1—lOkeV a partir des 
donnees XMM. Cadre du bas : Profil pulse radio base sur les observations du telescope GBT 
a la frequence de 2 GHz avec 128 intervalles en phase. Cette figure est extraite de \Abdo et aZ| 
1^009-12229) . 
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7.5. Discussion 



Une recherche d'une emission non-pulsee de la nebuleuse Boomerang a ete effectuee a 
partir du meme lot de donnees. Nous avons ajuste une source ponctuelle aux donnees de la 
region OFF a la position du pulsar radio entre 0.2 et 100 GeV. Aucun signal n'a ete observe 
de la PWN. Finalement, apres une renormalisation du signal sur toute la phase, nous avons 
derive une limite superieure a 95% de confiance sur le flux integre au-dessus de 200 MeV de 
4.0 X IQ-^cm-^s-i. 
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Figure 7.3 - Distribution d'energie spectrale pour le pulsar PSR J2229+61 14- La ligne 
continue represente I'ajustement globale des donnees par un maximum de vraisemblance en 
supposant une loi de puissance avec une coupure exponentielle (voir I'equation 7.1). Les 
points spectraux sont le resultat de I'ajustement des donnees a partir d'une simple loi de 
puissance pour six intervalles en energie separes logarithmiquement. 



7.5 Discussion 



7.5.1 Courbe de Lumiere 



La courbe de lumiere y de PSR J2229+6114 couvre (voir la flgure 7.2) un large intervalle 
en phase, suggerant que le faisceau y emet selon un grand angle solide. Cette caracteristique 
semble en faveur des modeles d'emission magnetospheriques externe, en particulier le Outer 
Gap (OG) et le Slot Gap/Two-Pole Caustic (SG/TPC) (voir la section [IA2|) . 



Ill 



Chapitre 7. PSR J 2229+6 114 



7.5.2 Luminosite, Efficacite, et Geometrie de rEmission 



A partir de I'analyse spectrale du pulsar, le flux d'energie integre observe est evalue a 
(23.7±0.7±2.5) xlO^" ergcm^^ La luminosite y resultante donne Ly = 2.6x 10^^fa(d/3 kpc) 
erg tandis que I'efficacite donne ri^ = 0.011 fa (d./3 kpc)^, qui e st un facteur 1.6 plu s 



>0.5 



0.007 (Watters et al 



2009) 



grand que I'efficacite estimee par la relation r\ ~ (10 /E 
Notons qu'a partir de la distance estimee de la SNR (0.8 kpc), I'efficacite diminue a 0.001 fa, 
tandis qu'avec la distance derivee a partir du DM (7.5 kpc) I'efficacite augmente a 0.07 fa, 
soulignant I'importance de determiner correctement la distance. 

A present, il est necessaire d'estimer le facteur de correction fa qui prend en compte la 
geometrie de remission. Nous pouvons comparer les courbes de lumiere a partir des modeles 



geometriques de | Watters et al (2009) pour deriver des contraintes sur la geometrie du faisceau. 
Les modeles magnetospheriques externes tels que TPC et OG avec riy ~ 0.01 peuvent decrire 
des proflls y avec un simple pic. Les solutions pour une telle conflguration apparaissent 
pour des inclinaisons magnetiques oc ~ 20° — 55°, et des angles d'observation C ~ 25° — 50° 
(geometrie TPC/SG) et pour cc ~ 45° - 80°, C ~ 35° - 70° et fa ~ 0.47 - 1.08 (geometrie 
OG). Pour PSR J2229+6114, nous avons une contrainte geometrique supplementaire a partir 
de la modelisation des tores des PWN en X entourant le pulsar, tels que mesures par le 



telescope Chandra. Ng and Romani (2004 2008) ont determine a partir des donnees X un 



angle d'observation C = 46° ± 2° ± 6°. Ceci donne un intervalle de recouvrement avec les 
solutions TPC et OG pour un simple pic. Finalement, bien qu'il y ait une certaine incertitude 
sur la determination de I'inclinaison magnetique a partir du profll radio, le decalage du pic y 
par rapport au pic radio (cj) ^ 0.50) fournit une information supplementaire. En examinant 



I'ecliantillon de courbe de lumiere des flgures 9 & 10 de Watters et al (2009), nous trouvons 
que les proflls a un seul pic pour cette phase sont determines seulement quand on pent 
identifler une composante P2, avec un pont ("bridge"), et un pic PI tres faible ou manquant. 
De telles courbes de lumieres apparaissent seulement pour des efficacites r\ < 0.03 et seulement 
pour un intervalle en angle de a ~ 50° - 70°, C ~ 40° - 50° (TPC/SG, avec PI faible) ou 
a > 45°, C = 45° — 50° (OG, avec PI absent). Ceci est precisement I'intervalle C exige par 
I'ajustement du tore de la PWN. En incluant la contrainte (3e = a - a partir de la detection 
de remission radio, nous voyons que les donnees radio, X, et y, sont toutes conformes avec 
a = 55° ±5°, C = 45 ± 5°, impliquant fa ~ 1 pour les deux modeles. Ainsi, pour observer un 
tel profll multi-longueurs d'onde, cela necessite une conflguration geometrique de remission 
bien precise, une faible efficacite y, et une emission a haute altitude. 
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8.1 Introduction 



Ce chapitre presente la decouverte de pulsations provenant du pulsar PSR J1048— 5832 
(PSR B 1046— 58) avec le telescope Fermi Gamma-Ray Space Telescope. Un signal y en pro- 
venance de ce pulsar avait ete prealablement detecte avec un faible niveau de confiance par 
EGRET (voir le tableauO). De plus, au meme titre que PSR J2229+6114, PSR J1048-5832 



est un pulsar de type Vela (voir la section 7.1). Dans la section 8.3 nous decrivons les ob 



servations radio et y utilisees pour I'analyse, tandis que les sections |8.4.1 et 8.4.2 incluent 
respectivement I'etude de la courbe de lumiere y fournissant Talignement avec remission 



radio, et I'etude spectrale de la source. Dans la section [8751 nous discutons de la courbe de lu- 
miere, de la luminosite, de I'efiicacite y, et de la geometrie de remission. Enfin, nous concluons 
ce chapitre par une comparaison des quantites mesurees entre ce pulsar et PSR J2229-I-6114. 



Ce chapitre est tire de I'article Abdo et al (2009-J2229), a paraitre dans ApJ en 2009. 
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8.2 Candidat pour la Source EGRET 3EG J1048-5840 



PSR J 1048-5832 (B 1046-58) est localise dans la region de la nebuleuse de Carina (NGC 
3372) a basse latitude galactique (I = 287.42°, b = 0.58°). II a ete decouvert durant une 
recherche de pulsations a 1.4 GHz dans le plan galactique par le radiotelescope de Parkes, avec 



une periode de rotation P de 123.7ms (Johnston et al 1992). Le pulsar a un ralentissement 
E de 2 X 10^^ erg^^ pour un moment d'inertie I de lO'^^gcm^, un champ magnetique dipolaire 
de surface de 3.5 x 10^^ G, et un age caracteristique Tj. de 20000 ans. Les observations a 
haute resolution d'une source ASCA[^ coi'ncidant avec le pulsar par les telescopes Chandra 
et XMM- Newton ont revele une nebuleuse a vent de pulsar asymetrique de - 6" x 11". 
Cependant, aucune pulsation a partir des donnees n'a ete trouvee ( ] Gonzalez et al 2006). 



PSR J1048-5832 a egalement ete propose comme etant la contrepartie de la source EGRET 
non-identifiee 3EG J1048-5840 (|Fierro[ [l995t iPivovaroffetall [205oj [Nolan et al[ [2003)). Cette 



conclusion a ete suggeree par la coincidence de la position du pulsar avec la source EGRET 
et des proprietes spectrales et energetiques de la source EGRET elle-meme. Une analyse 
detaillee des donnees EGRET montre une possible pulsation y au-dessus de 400 MeV, avec 
une courbe de lumiere composee de deux pics separes par ~ 0.4 en phase (voir la figure 8.2). 
L'observation d'H I dans I'environnement du pulsar donne une distance comprise entre 2.5 



et 6.6 kpc (Johnston et al 1996). Le modele NE 2001 assigne quant a lui une distance de 



2.7 kpc basee partiellement sur la determination de la distance par l'observation H L Pour ce 
travail de these, nous adopterons une distance de 3 kpc. 



8.3 Observations 



8.3.1 Observations Radio 

Le pulsar PSR J1048— 5832 est observe mensuellement par le radiotelescope de 64 m de 



Parkes en Australie depuis 2007 (Manchester 2008). La duree d'une observation est de 2 
minutes a la frequence de 1.4 GHz avec occasionnellement des observations a 0.7 et 3.1 GHz. 
Des informations detaillees sur la procedure et sur I'analyse des donnees peuvent etre trouvees 



dans Weltevrede ([2009). Le pulsar est connu pour avoir une forte instabilite rotationnelle ; 
un glitch a ete enregistre juste avant le decoUage de Fermi (Weltevrede 2009). De plus, 
comme beaucoup de pulsars presentant un fort taux de ralentissement E, PSR J 1048— 5832 



est hautement polarise sur la large bande radio observee ( Karastergiou et al, 2005a; Johnston 



et al 2006). L'ephemeride utilisee pour empiler les photons y a la bonne phase est basee 



sur vingt TOAs et derive du programme d'analyse TEMP02. La solution est ajustee pour la 
frequence de rotation du pulsar et ses derivees, et lissee en utilisant les algorithmes nommes 
fitwaves de TEMP02. L'ecart-type (RMS) resultant de 287|xs presente une amelioration sur 
la solution pre-lissee. Notons que les parametres fitwaves ne sont pas pris en charge par le 
ST gtpphase. Des mesures du retard du a la dispersion a 1369 MHz donnent une mesure 
de la dispersion (DM) de 128.822±0.008pc cm^^, sans indication que celle-ci varie au cours 
du temps. Cette mesure est utilisee pour corriger les TOAs vers les frequences infinies pour 
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calculer la phase des photons y. 



Parametres Valeurs 

Ascension Droite (J2000) 22'^29"^05. 28000084^ 

Declinaison (J2000) +6ri4'09.3000042" 

Frequence (s^^) 8.08403505(1) 

Derivee premiere de la frequence (s^^) -6.282(1) xlO^^^ 

Epoque, To (MJD[TDB]) 54575.0 

Date de reference pour le 4)o (MJD[TDB])'' 54743.974102133058178 

Intervalle de validite en MJD 54565.3 - 54931.3 

Nombre de TOAs 20 

Mesure de dispersion, DM(cm"^pc) 128.822 ± 0.008 

Modele d'ephemeride du S.S DE405 



Tableau 8.1 - Parametres rotationnels du pulsar PSR J1048— 5832 bases sur les observations 
du telescope Parkes. Les chiffres entre parentheses representent les incertitudes donnees par 
TEMP02. 

Date a laquelle la phase cj) = (TZRMJD). 



8.3.2 Observations Gamma 



Les donnees utilisees pour I'analyse spectrale ont ete coUectees a partir du 4 aout 2008 
jusqu'au 10 avril 2009 pendant que le telescope etait en mode de balayage. Pour I'analyse 
temporelle, nous avons additionne a cet echantillon de photons les donnees de la phase de 
calibration de I'appareil, periode s'etendant du 30 juin au 3 aout 2008. Seuls les photons 
marques diffuse ont ete selectionnes, et nous avons par ailleurs exclu ceux possedant un 



angle zenithal superieur a 105°. La figure 8.1 montre la distribution des photons de la region 



entourant le pulsar. Notons la presence des pulsars y PSR Jl 124— 5916 (Abdo et al 2009 



psrcat 


), PSR J1028-5819 ( 


Abdo et al 


(Abdo et al 


2009-EGRET 


. 



8.4 Analyse et Resultats 

8.4.1 Courbe de Lumiere 

Pour I'analyse temporelle de PSR J1048— 5832, nous avons selectionne les photons avec 
une energie superieure a 100 MeV dans un rayon de 1 degre autour de la position du pulsar 
radio. Ensuite, nous avons selectionne les photons dans un cone dependant en energie de 
rayon Ogg ^ 0.8 x Eq^y^ degres, gardant tous les photons inclus dans un rayon de 0.35°. Cette 
selection se base sur la PSF du LAT et permet d'augmenter le rapport signal-sur-bruit dans 
une large bande en energie. A partir du programme d'analyse TEMP02, nous avons corrige 
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Figure 8.1 - Distribution des photons de la region qui entoure le pulsar PSR J1048— 5832. 
Les croix representent les sources du catalogue 9 mois. Le nombre de degre par pixel est de 
0.25. 



le temps d'arrivee des photons vers le barycentre du systeme solaire en utilisant I'eplienieride 
solaire JPL DE405. La phase des evenements y a ensuite ete calculee a partir de Tephemeride 
radio delivree par le telescope de Parkes (voir le tableau 8.1). Le test de periodicite H-test 
(voir la section 5.2.4) sur cet echantillon de photon donne une probabilite ^ 4 x 10^^ que la 
pulsation soit une fluctuation. Cette valeur est par 4 ordres de grandeur plus petite que les 
precedents resultats d'EGRET ( |Kaspi et al 2000) et etablit done fermement la detection du 
pulsar. 
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La figure 8.2| (cadre du haut) montre rhistogramme du nombre de photons y detectes 
en fonction de la phase au-dessus de 0.1 GeV. Le profil radio a la frequence de 1.4 GHz est 
indique dans le cadre du bas pour la comparaison des phases. La courbe de lumiere y est 
composee de deux pics, avec le pic PI a la phase 0.15 ± 0.01 ± 0.000^ et le pic P2 a la phase 
0.57±0.01 ±0.000^. La separation en phase Acf* des deux pics est ainsi de 0.42±0.01 ±0.000^. 
Les erreurs sont respectivement I'incertitude sur Tajustement et I'incertitude sur le DM. Ces 



resultats sont en accord avec ceux trouves par Kaspi et al (2000) 



Les deux pics apparaissent asymetriques, avec un cote montant doux et un cote descendant 
prononce. Nous avons ainsi ajuste chacun des pics avec deux demi-lorentziennes pour prendre 
en compte ces deux differences de pente. PI ((() = 0.05 — 0.17) est plus fin avec une largeur a 
mi-hauteur de 0.06 ± 0.01 , tandis que P2 {(p = 0.45 — 0.65) est plus large avec une largeur a 
mi-hauteur de 0.10 ± 0.02. Notons qu'une structure qui peut etre associee a un epaulement 
ou a un "bridge" entre les pics est presente entre 0.17 et 0.30 en phase. Ce profil est tres 
similaire a la courbe de lumiere de Vela qui elle-meme est composee de deux pics separes par 



0.43 en phase et d'un bridge (Abdo et al 2009- Vela). Nous avons egalement defini la region 



OFF comme le minimum de I'histogramme (cf) = 0.7— 1.05). Pour verifier cette hypothese, 
nous avons estime le fond represents par une ligne pointillee (73 coups/inter valle) a partir 
d'un anneau centre sur le pulsar entre 1° et 2°. Les sources proches ont ete exclues, et nous 
avons renormalise I'anneau au meme espace de phase que notre selection. Le resultat est en 
bon accord avec la region OFF. Ainsi, le nombre total de photons pulses provenant du pulsar 
est estime a 933 ± 93, avec une contribution du fond de 3654 ± 60 evenements. 



La figure 8.2 (cadres du milieu) montre I'histogramme du nombre de photons y en fonction 
de la phase pour quatre intervalles en energie differents (0.1 - 0.3 GeV, 0.3 - 1 GeV, 1-3 GeV, 
>3 GeV). L'evolution en energie de la forme de la courbe de lumiere est visible, bien que plus 
de donnees soient necessaires pour contraindre la largeur des pics pour chaque intervalle. 
En dessous de 0.3 GeV notamment, le premier pic semble plus large qu'a haute energie. 
Cette caracteristique pourrait etre expliquee par la contamination de la source y (EMS 0693) 
a moins de 1° du pulsar. Nous avons egalement regarde revolution du rapport P1/P2 en 
fonction de I'energie, egalement observe par EGRET pour les pulsars de Vela, du Crabe, de 
Geminga, de PSR B1951+32 (Fierro[ fidib) et par le LAT pour Vela (Abdo et al[ |2009-VelaD 



et PSR J0205-I-6449 (voir le chapitre 6j). Dans chaque bande en energie, nous avons calcule la 
hauteur du pic en respectant le niveau de fond estime a partir de la region OFF. Le rapport 
des hauteurs montre qu'il n'y a aucune variation avec un niveau de confiance de 69%. II 
convient d'avoir plus de donnees pour etre plus contraignant. Enfin, notons qu'entre 1 et 
3 GeV et au-dessus de 3 GeV I'epaulement apres PI (cf) = 0.17 - 0.33) ainsi que les 2 autres 
pics sont toujours significatifs, et que le photon de plus haute energie est detecte dans PI a 
19 GeV. 

8.4.2 Spectre et Flux Resolus en Phase 

Afin d'obtenir le spectre resolu en phase et le flux integre de PSR J1048— 5832, une analyse 
spectrale par un maximum de vraisemblance a ete effectuee en utilisant le ST "gtlike" (voir 



la section 5.3). Nous avons utilise les IRFs "Pass6_v3" lesquelles sont, rappelons-le, une mise 
a jour apres le lancement du satellite qui corrige des effets de traces fantomes dans le LAT. 
Une region circulaire de 15° centree sur la position du pulsar a ete modelisee en incluant 
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1.4 1.6 1.8 2 



Pulsar Phase 



Figure 8.2 - Cadre du haut : Courbe de lumiere de PSR J1048—5832 au-dessus de 0.1 GeV. 
Deux rotations sont representees avec 50 intervalles en phase par periode (P — M3.7 ms). La 
ligne pointillee montre le niveau de fond estime avec un anneau centre sur la position du 
pulsar (73 coups/intervalle) . Quatre cadres suivants : Histogrammes du nombre de photons 
detectes en fonction de la phase pour quatre intervalles d'energie. Chacun est represents avec 
50 intervalles en phase. Second cadre a partir du bas : Courbe de lumiere d'EGRET au- 
dessus de 400 MeV (Kaspi et al 2000^ . Cadre du bas : Profil radio pulse enregistre par le 
telescope de Parkes a la frequence de I.4 GHz (1024 intervalles de phase). 
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les sources proches ainsi que remission diffuse. Les pulsars inclus dans la region ont ete 



ajustes par une loi de puissance avec une coupure exponentielle (equation 8.1), tandis que 
les autres sources ont ete ajustees par une simple loi de puissance. Le fond diffus galactique 
a ete pris en compte en utilisant une carte basee sur le modele numerique GALPROP appele 
54_59Xvarli7S. Le fond isotropique a ete quant a lui modelise avec un spectre tabule qui 
derive d'un ajustement des donnees du LAT a haute latitude galactique. Nous avons ajuste 
le spectre de PSR J1048— 5832 par une loi de puissance avec une coupure exponentielle en 
energie entre 0.1 et 0.7 en phase. Le spectre est decrit par I'equation : 

— = No E-^e-^/^^ cm-^s-^GeV-i (8.1) 
dh 

avec E en GeV, le terme No = (5.9 ± 0.3 ± 0.1 ) x 1 0'^ cm"^ s"^ GeV"^ , un indice spectral V = 
1.38 ±0.06 ±0.12 et une coupure en energie E^ — 2.3^q 4 ± 0.3 GeV. La premiere incertitude 
est statistique, tandis que la seconde est systematique (utilisation des IRFs modifiees). A 
partir de ce resultat, nous avons obtenu un flux de photon integre entre 0.1 - 100 GeV de 
(2.1 9 ± 0.22 ± 0.32) x 1 0^^ cm^^ s^^ , lequel correspond a environ un tiers du flux de la source 



EGRET 3EG J1048-5840. La figure [O] montre a la fois I'ajustement global entre 0.1 et 20 
GeV (ligne continue), et les points spectraux obtenus par une analyse dans six intervalles en 
energie espaces logarithmiquement, en supposant a chaque fois une simple loi de puissance. 
Pour verifier I'hypothese de la coupure en energie dans le spectre, nous avons egalement 
modelise le pulsar avec une simple loi de puissance de la forme dE/dE = NoE^^. Le modele 
spectral utilisant la coupure est mieux contraint avec une difference dans les logarithmes du 
maximum de vraisemblance de ~ 1 ff . 



Le pulsar est repertorie dans le Bright Source List du Fermi LAT (Abdo et al 2009-bsl[ ) 



comme la source OFGL J1047.6-5834, qui est localisee a (RA,Dec) = (161.922°, -58.577°) 
avec une extension spatiale de rayon 0.138° a 95% de niveau de confiance. Une autre source 
LAT est localisee a ~1°, OFGL J1045.6— 5937. Nous avons modelise cette source non-identifiee 
par une simple loi de puissance. Le meilleur ajustement donne un indice spectral de 2.2 ± 
0.1 ± 0.1 et un flux de photon integre de (4.49 ± 0.40 ± 0.80) x lO^^cm^^s^^ La somme 
du flux de cette source et de PSR J1048— 5832 est ~ 6.7 xlO^^ cm^^s^\ identique au flux 
de la source EGRET 3EG J1048-5840 de (6.2 ± 0.7) x ]0-^cm-^s-\ La source EGRET 



3EG J 1048-5840 ainsi que la source COS-B 2CG 288-00 ^Swanenburg et al 1981) etaient 



apparemment constitutes de ces deux sources, qui sont a present resolues par le Fermi LAT. 



8.5 Discussion 



8.5.1 Courbe de Lumiere 



La courbe de lumiere y de PSR J1048— 5832 (voir la figure 8.2) couvre un large intervalle 
en phase, suggerant que le faisceau y couvre un grand angle solide. Le profil est tres similaire 



a celui de Vela ( |Abdo et alj |2009-Vela[ ) ou meme de PSR J0205+6449 (voir le chapitre |6]). 
Cette caracteristique semble en faveur des modeles d'emission magnetospheriques externes, 
en particulier le Outer Gap (OG) et le Slot Gap/Two-Pole Caustic (SG/TPC) (voir la 
section 1.4.2). 
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Figure 8.3 - Distribution d'energie spectrale pour le pulsar PSR J1048— 5832. La ligne conti- 
nue represents I'ajustement global par une loi de puissance avec une coupure exponentielle. 
Les points spectraux ont ete obtenus par une analyse dans six intervalles en energie espaces 
logarithmiquement, en supposant pour chaque intervalle une simple loi de puissance. Les 
erreurs sont uniquement statistiques. 



8.5.2 Luminosite, Efficacite, et Geometrie de rEmission 



sur 

1035 



Pour PSR J 1048— 5832, le flux d'energie integre obtenu en integrant I'equation |8.1 
I'energie est (19.4 ± 1.0 ± 3.1) x 10^" ergcm^^s^\ menant a une luminosite y de 2.1 x 
fa (d/3 kpc)^ et une efficacite r| = L^/E = 0.10 fa (d/3 kpc)^. 

La difference de phase des deux pics y et I'efficacite sont utiles pour placer des contraintes 
sur les modeles d'emission magnet ospheriques externes. En utilisant les courbes de lumiere 



y de "1' Atlas" (jWattersetal , 2009), nous avons estime I'intervalle autorise pour I'inclinaison 
magnetique a ~ 60° - 85° , Tangle d'observation C ~ 70° - 80° , et le facteur de correction 
fa ~ 0.7 - 1.1 pour le modele OG, et a ~ 50° - 75° , C ~ 50° - 75° , et fa ~ 1.1 pour le 
modele SG/TPC. De plus, la variation sur Tangle de polarisation radio contraint Tangle |3e = 
Ce- En ajustant le modele de vecteur rotationnel (Radhakrishnan and Cooke, 1969) pour 



les donnees radio de PSR J1048— 5832, une valeur de |3e plus petite que 10° a ete derivee 



( Karastergiou et aH 2005b Weltevrede and Johnston 2008a). Les deux modeles theoriques 



ont ainsi un bon intervalle de solutions pos s ibles, mais cela suppose une efficacite basee sur 



la relation r\ ~ (10^^/E)°-^ de Watters et al (2009), laquelle donne rj ~ 0.02, soit un facteur 



cinq fois plus petit que I'efficacite derivee a partir de la luminosite mesuree pour fa = 1 • En 
supposant notre quantite mesuree r\ = 0.10, I'intervalle autorise reste le meme pour le modele 
TPC et I'intervalle de <x se deplace vers ~ 70° - 90° pour le modele OG. Le decalage de la 
phase entre le pic radio et le premier pic y est ainsi en bon accord avec les deux modeles 
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8.5. Discussion 



theoriques. 



8.5.3 Comparaisons de PSR J1048-5832 et PSR J2229+6114 

Bien que PSR J1048-5832 et PSR J2229+6114 soient deux pulsars de type Vela, leur profil 



y et la geometrie de leur emission sont assez differents. Le tableau 8.2 resume les principales 
quantites mesurees pour les deux pulsars. La courbe de lumiere de PSR J1048— 5832 est 
tres similaire a celle de Vela, tandis que son efficacite est d'un facteur 10 x plus grand a 
celui-ci (en adoptant f a = 1 et d = 3kpc). Au contraire, I'efEcacite de PSR J2229+6114 est 
tres similaire a celle de Vela, mais le pulsar ne presente qu'un seul et large pic, identique a 



PSR J0357+32 decouvert avec les donnees LAT (Abdo et al 2009-BS). 



Analyse 



Parametres 



PSR J1048-5832 



PSR J2229+6449 



Temporelle 



Nombre de photons pulses 
Position des pics (c})) 

Separation des pics (A) 
FWHM(pic) 



Spectrale 



''F (10-^ cm-^s-') 
^Fe (10"" erg cm" 

r 

'^Ec (GeV) 

(10^5 erg s-') 

JTy 



933 ± 93 
0.15 ± 0.01 ± 0.0001 (PI) 
0.57 ± 0.01 ± 0.0001 (P2) 
0.42 ± 0.01 ± 0.0001 
0.06 ± 0.01 (PI) 
0.10 ± 0.02 (P2) 



1365 ± 97 
0.49 ± 0.01 ± 0.001 



0.23 ± 0.03 



2.19 ± 0.22 ± 0.32 
19.4 ± 1.0 ± 3.1 

1.38 ± 0.06 ± 0.12 

2.3+°-3 ± 0.3 
2.1 fa (d/3kpc)2 

0.10 fa (d/3kpc)2 



3.77 ± 0.22 ± 0.44 
23.7 ± 0.7 ± 2.5 

1.85 ± 0.06 ± 0.10 

3.6+°-^ ± 0.6 
2.6 fa (d/3kpc)2 

0.011 fa (d/3kpc)2 



Tableau 8.2 - Ce tableau resume les resultats des analyses spectrale et temporelle des deux 
pulsars de type Vela, PSR J1048- 5832 et PSR J2229+6114. Les erreurs statistiques et 
systematiques sont donnees. 



121 



Chap it re 8. PSR J 1048-5832 



122 



Chapitre 9 



Catalogue de Pulsars Gamma 



Sommaire 

9.1 Introduction [T23l 

9.2 Observations et Analyses 11251 

9.2.1 Analyse Temporelle [125] 

9.2.2 Analyse Spectrale [T26l 

9.3 Associations 11301 

9.4 Resultats et Discussion 11301 

9.4.1 Luminosite Gamma 11301 

9.4.2 Population [T33l 

9.4.3 Courbes de Lumiere et Energie de Coupure 11331 



9.1 Introduction 

Les capacites du Large Area Telescope avec sa sensibilite nettement amelioree compare 
a ses predecesseurs, nous donnent I'opportunite d'etudier une population nouvelle de pulsars 
emetteurs y. Ce chapitre presente les 46 pulsars y detectes avec un haut niveau de confiance a 
partir des six premiers mois des donnees. Cela inclut les six pulsars observes clairement par le 
telescope EGRET, plus les trois pulsars detectes marginalement par EGRET (voir le tableau 



1.1). Nous presenterons essentiellement les proprietes des pulsars normaux detectes (non 



milliseconde). La section 9.2 decrit brievement les observations et les procedures d'analyse 



temporelle et spectrale, tandis que dans la section |9.4| nous discutons des resultats sur la 
luminosite, la population de pulsars, et les courbes de lumiere. Signalons que ce chapitre 
est extrait du premier catalogue de pulsars avec le telescope Fermi Gamma-Ray Space 



Telescope (Abdo et al 2009-psrcat). 
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PSR 


Type, 


I 


b 


P 


P 


age Tc 


E 


Blc 


Sl400 




Ref. 


(°) 


(°) 


(ms) 




(kyr) 


(1034 gj.g g-l •) 


(kG) 


(mJy) 


J0007+7303 


g 


119.7 


10.5 


316 


361 


14 


45.2 


3.1 


' < 0.1 


J0030+0451 


m '^-'^ 


113.1 


-57.6 


4.9 


10 xlQ-'' 


7.7x106 


0.3 


17.8 


0.6 


J0205+6449 


I " 


130.7 


3.1 


65.7 


194 


5 


2700 


115.9 


0.04 


J0218+4232 


mb 


139.5 


-17.5 


2.3 


77x10-6 


0.5x106 


24 


313.1 


0.9 


J0248+6021 


r f 


137.0 


vO.4 


217 


55.1 


63 


21 


3.1 


9 


J0357+32 


g" 


162.7 


-16.0 


444 


12.0 


590 


0.5 


0.2 




J0437-4715 


mb <^ 


253.4 


-42.0 


5.8 


14x10-'' 


6.6xl06 


0.3 


13.7 


140 


J0534+2200 


r 1^ 


184.6 


-5.8 


33.1 


423 


1 


46100 


950.0 


14 


J0613-0200 


mb 


210.4 


-9.3 


3.1 


9nxlO-<' 


5.3x106 


1.3 


54.3 


1.4 


J0631+1036 


r ^ 


201.2 


0.5 


288 


105 


44 


17.3 


2.1 


0.8 


J0633+0632 


g" 


205.0 


-1.0 


297 


79.5 


59 


11.9 


1.7 


2 < 0.2 


J0633+1746 


g*^ 


195.1 


4.3 


237 


11.0 


340 


3.3 


1.1 


< 1 


J0659+1414 


r ^ 


201.1 


8.3 


385 


55.0 


110 


3.8 


0.7 


3.7 


J0742-2822 


r ^ 


243.8 


-2.4 


167 


16.8 


160 


14.3 


3.3 


15 


J0751+1807 


mb 


202.7 


21.1 


3.5 


6.2x10-'' 


8.0x106 


0.6 


32.3 


3.2 


J0835-4510 


I k 


263.6 


-2.8 


89.3 


124 


11 


688 


43.4 


1100 


J1028-5819 


r I 


285.1 


-0.5 


91.4 


16.1 


90 


83.2 


14.6 


0.36 


J1048-5832 


r 


287.4 


0.6 


124 


96. 


20 


201 


16.8 


6.5 


J1057-5226 


r 1^ 


286.0 


6.6 


197 


5.8 


540 


3.0 


1.3 


11 


J1124-5916 


r 


292.0 


1.8 


135 


747 


3 


1190 


37.3 


0.08 


J1418-6058 


g" 


313.3 


0.1 


111 


170 


10 


495 


29.4 


<0.06 


J1420-6048 


r ^ 


313.5 


0.2 


68.2 


83.2 


13 


1000 


69.1 


0.9 


J1459-60 


g" 


317.9 


-1.8 


103 


25.5 


64 


91.9 


13.6 


2 < 0.2 


J1509-5850 


r ^ 


320.0 


-0.6 


88.9 


9.2 


150 


51.5 


11.8 


0.15 


J1614-2230 


mb 


352.5 


20.3 


3.2 


4x10-'' 


1.2xl06 


0.5 


33.7 




J1709-4429 


r 1^ 


343.1 


-2.7 


102 


93.0 


18 


341 


26.4 


7.3 


J1718-3825 


r ' 


349.0 


-0.4 


74.7 


13.2 


90 


125 


21.9 


1.3 


J1732-31 


g" 


356.2 


0.9 


197 


26.1 


120 


13.6 


2.7 


2 < 0.2 


J1741-2054 


gb,t 


6.4 


4.6 


414 


16.9 


390 


0.9 


0.3 


2 0.16 


J1744-1134 


m d 


14.8 


9.2 


4.1 


7x10-6 


9xl06 


0.4 


24.0 


3 


J1747-2958 


r " 


359.3 


-0.8 


98.8 


61.3 


26 


251 


23.5 


0.25 


J1809-2332 


g" 


7.4 


-2.0 


147 


34.4 


68 


43.0 


6.5 


< 0.06 


J1813-1246 


g" 


17.2 


2.4 


48.1 


17.6 


43 


626 


76.2 


2 < 0.2 


J1826-1256 


g" 


18.5 


-0.4 


110 


121 


14 


358 


25.2 


2.3 < 0.06 


J1833-1034 


r ° 


21.5 


-0.9 


61.9 


202 


5 


3370 


137.3 


0.07 


J1836+5925 


g" 


88.9 


25.0 


173 


1.5 


1800 


1.2 


0.9 


4 < 0.007 


J1907+06 


gb,r 


40.2 


-0.9 


107 


87.3 


19 


284 


23.2 


< 0.02 


J1952+3252 


r 


68.8 


2.8 


39.5 


5.8 


110 


374 


71.6 


1 


J1958+2846 


^ b 
g 


65.9 


—0.2 


290 


222 


21 


35.8 


3.0 




J2021+3651 


r P 


75.2 


0.1 


104 


95.6 


17 


338 


26.0 


0.1 


J2021+4026 


gb,S 


78.2 


2.1 


265 


54.8 


77 


11.6 


1.9 




J2032+4127 


gb.t 


80.2 


1.0 


143 


19.6 


120 


26.3 


5.3 


2 0.24 


J2043+2740 


r 1 


70.6 


-9.2 


96.1 


1.3 


1200 


5.6 


3.6 


7 


J2124-3358 


m d 


10.9 


-45.4 


4.9 


12x10-6 


0.6xl06 


0.4 


18.8 


1.6 


J2229+6114 


r 


106.6 


2.9 


51.6 


78.3 


11 


2250 


134.5 


0.25 


J2238+59 


g" 


106.5 


0.5 


163 


98.6 


26 


90.3 


8.6 





Tableau 9.1 - Parametres mesures et intnnseques des pulsars detectes par le LAT. Type : (r) pulsars normaux 
decouverts en radio, (g) pulsars decouverts en gamma, (m) pulsars milUseconde, (b) mdiquent s'lls appartiennent 
a des systemes binaires. S1400 densite de flux radio a I4OO MHz (source : ATNF) ou limite supeneure avec 1 



References : a /[Abdo et al 


\2008jj; b (Abdo et al, 


{ Abdo et al 


\2009-J020S\) ; f 


Cognard et ai, ^2009^ ; h 


1 Abdo et al 


\2009-gemmga^ ; k lAbdo et al\\2009-Vela\l 


Abdo et al 


2009-EGRET) ; (Camilo et alV2009^; } 



Abdo et al\ \ 2009- Crab') ; i (Welievrede and Wright'. 2009) ; j 
TjAbdo et al. 2009-J1028) ; m (Abdo et al. 2009-J2229) ; n 



[eTai, 2009-J1907b) ; s (Abdo et al. 2009-J1907af;t {Camilo et al| \200^ 



Abdo et al^ ^2009-12021( 1 ; q (Noutsos et al^ ^2009) ; r (Abdo\ 
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9.2 Observations et Analyses 

Dans le but d'avoir une analyse uniforme, les donnees y utilisees pour I'etude de cette 
population de pulsars presentent le meme intervalle de temps d'observation. Pour I'analyse 
spectrale, les donnees ont ete coUectees a partir du debut de I'observation du telescope Fermi 
en mode balayage (4 aout 2008) jusqu'au 1 fevrier 2009, tandis que I'analyse temporelle 
commence a partir du premier evenement enregistre par le LAT (25 juin 2008) jusqu'au 1 
fevrier 2009. Comme dans les chapitres precedents, seuls les evenements classes diffuse ont 
ete selectionnes, et nous avons exclu ceux possedant un angle zenithal superieur a 105°. 



9.2.1 Analyse Temporelle 



Deux methodes pour la recherche de pulsations dans les donnees LAT ont ete mises en 
place des le debut de la mission. 

La premiere methode recherche des pulsations uniquement a partir du temps d'arrivee des 
photons en selectionnant des candidats tels que des SNRs et des PWNs suspectees d'abriter un 



pulsar, et sur des sources du catalogue du LAT (Atwood et al 2006). Au total apres six mois 
d'observation, 16 nouveaux pulsars ont ete decouverts parmi lesquels le pulsar de la jeune 
supernova galactique CTAl ( |Abdo et al[ |2008[ |2009-BS[ ). Soulignons deux points : (1) les 



modeles theoriques d'emission predisent une large emission y en eventail et un fin faisceau 
radio, cette methode favorise done les pulsars avec une grande inclinaison magnetique, si 
Ton suppose que remission radio vient du pole magnetique ; (2) les pulsars de cette nouvelle 
population ne sont pas forcement des pulsars de type Geminga {radio -faint). Des recherches 
de pulsations radio ont ete effectuees sur certains de ces pulsars. Parmi eux, trois pulsars 



ont ete observes avec un faible signal radio, PSRs J1741-2054 et J2032+4127 ( [Camilo et al 
20091), PSR J1907+0602 (|Abdo et al[ |2009-J1907b|). 



La deuxieme methode est la recherche de pulsations y en utilisant les ephemerides de 
pulsars deja connus dans une autre longueur d'onde. Comme nous I'avons vu a la section 
2.3. 1[ il est essentiel pour pouvoir suivre et detecter des pulsars d'avoir des ephemerides 



valides sur la periode d'observation y. Au final, les astronomes de la radio et des X nous 
ont fourni 767 ephemerides contemporaines. Une premiere partie des ephemerides est notre 
selection de 224 pulsars candidats y avec E > lO^'^ergs^^ (voir le chapitre [2]), tandis que 
I'autre partie est un echantillon de 543 pulsars avec E < 10^'^ergs^^ permettant ainsi de 
reduire le biais possible dans notre recherche de pulsations. Plusieurs programmes d'analyse 



bases sur les tests de periodicite (section 5.2.4) tel que le H-test ont ete mis en place dans 



la collaboration Fermi et notamment a Bordeaux pour suivre ces pulsars (Guillemot and 



Parent, 2007). Les pulsars presentant une significativite superieure a 5 a pour les 6 premiers 



mois des donnees LAT ont ete selectionnes pour cette etude. Cette methode a I'avantage 
d'etre sensible a des flux plus faibles, et permet de comparer I'alignement en phase des profils 
aux differentes longueurs d'onde, ce qui apporte une excellente contrainte sur la geometrie 
de remission. 

Le tableau [971] liste les parametres mesures et intrinseques des 46 pulsars y detectes. Les 
pulsars sont repertories selon 3 classes : (r) identifie les pulsars normaux decouverts en radio, 
(g) les pulsars decouverts en gamma, et (m) identifient les pulsars milliseconde, (b) indiquent 
s'ils appartiennent a des systemes binaires. La derivee de la periode P est corrigee de I'effet 
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cinematique Shklovskii (Shklovskii, 1970) : P = Fobs ~ l^^Pobsd/c^, ou \l est le mouvement 



propre du pulsar et d la distance. Cette correction s'addresse essentiellement aux pulsars 



milliseconde ( |Abdo et al[ |2009a[ ). 
Courbes de lumiere 

Un echantillon des courbes de lumiere des 46 pulsars detectes par le LAT sont presentees 
dans I'annexe de cette these. La selection des evenements et la construction des phases pour 
chaque pulsar sont similaires a I'analyse decrite dans les chapitres precedents pour les pulsars 
PSRs J0205+6449, J1048-5832, et J2229+6114. Pour un petit nombre de pulsars localises 
pres d'une source y brillante, la ROI a ete redefinie a 0.5° au lieu de 1°. Soulignons que les 
ephemerides utilisees pour les 16 pulsars decouverts en y, ainsi que pour le pulsar Geminga 
et le pulsar radio PSR J1124— 5916[^ derivent de modeles de rotation construits a partir des 



donnees du LAT (Ray and al. 2009). Ces modeles sont bases sur des observations s'etendant 



du 25 juin 2008 au ler mai 2009. Les ephemerides utilisees pour les pulsars emetteurs radio 



ont ete fournies par les telescopes enumeres a la section 2.1.1 



Pour illustrer ce travail, la figure |9.1| represente le resultat de I'analyse temporelle pour 
le pulsar radio PSR J1028-5819 (lAbdo et all |2009-J1028p. Le profil gris en haut represente 



I'histogramme du nombre de coups en fonction de la phase au-dessus de 0.1 GeV, tandis que 
les autres cadres montrent le profil y dans des intervalles d'energies : 0.1 - 0.3 GeV, 0.3 - 
1 GeV, >lGeV (avec les photons superieurs a 3 GeV representes en noir). Pour les pulsars 
emetteurs radio, le profil radio aligne en phase est montre dans le cadre le plus bas. Le nombre 
d'intervalles en phase pour chaque pulsar a ete defini en fonction de la statistique et de la 
RMS des residus de la solution temporelle, telle que la resolution d'un intervalle en phase fft 
soit inferieure a P/N-\/T2, avec IM = 25 ou 50 intervalles et P la periode de rotation du pulsar. 
La ligne pointillee sur le profil au-dessus de 0.1 GeV represente le niveau de fond estime avec 
un anneau entre 1° < 9 < 2° entourant la source. Pour certains pulsars le niveau de fond ne 
correspond pas au minimum du phasogramme. La resolution spatiale du LAT etant large a 
basse energie, dans le cas ou le pulsar se trouve tres proche d'une autre source, I'analyse est 
defaillante. Get effet est observe pour les pulsars PSR J1420-6048 et PSR J1418-6058 situe 
a moins de 0.2 ° I'un de I'autre. II peut egalement s'agir de la presence des PWNe K3 et du 
Rabbit. 

Si Ton regarde la dependance en energie des courbes de lumiere a la fois pour les pulsars 
de type (r) et de type (g), le rapport P1/P2 decroit quand I'energie augmente, ceci pour une 
grande partie des pulsars avec deux pics. Pour les pulsars avec un A ~ 0.5, tels que PSRs 
J1028-5820, J2021+3651, J0633+0632, J1124-5916, J1813-1246, J1826-1256, J1836+5925, 
et J2238-I-59, cette tendance n'est pas observee, suggerant une symetrie du profil y. 



9.2.2 Analyse Spectrale 

Pour etudier le spectre des pulsars et extraire le fiux integre au-dessus de 100 MeV de la 
region pulsee des 46 pulsars detectes, une analyse spectrale par un maximum de vraisemblance 
a ete effectuee en utilisant le catalogue des 6 mois de donnees. La source y Cygnus X-3 tres 

1. PSR J1124— 5916 etant tres faible en radio (une mesure demande environ 4h d'observation au telescope 
Parkes), il est plus facile de determiner des TOAs a partir des donnees LAT. L'ephemeride contient un TOA 
radio necessaire pour avoir le decalage entre les pulsations radio et y. 
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Figure 9.1 - Cadre du haut : Courbe de lumiere de PSR J1028—5819 au-dessus de 0.1 GeV. 

Deux rotations sont representees avec 50 intervalles en phase par periode (P ~ 9]. 4ms). La 
ligne pointillee montre le niveau de fond estime avec un anneau centre sur la position du 
pulsar (42 coups/intervalle). Trois cadres suivants : Histogrammes du nombre de photons 
detectes en fonction de la phase pour quatre intervalles d'energie. Le nombre de coups au- 
dessus de 3 Ge V est represente en noir. Cadre du bas : Profil radio pulse enregistre par le 
telescope de Parkes a la frequence de I.4 GHz (1024 intervalles de phase). 
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proche du pulsar PSR J2032+4127 a egalement ete modelisee, bien qu'elle n'ait pas ete 
detectee automatiquement dans la creation du catalogue. L'analyse etant similaire a celle des 
pulsars etudies aux chapitres 6 a 8, nous n'en rappellerons que les points essentiels. 

Tous les pulsars ont ete modelises par une loi de puissance avec une coupure exponen- 



tielle (equation 6.1 ). Le terme de normalisation Nq, I'indice spectral V, et I'energie de coupure 
tcutoff sont laisses libres dans I'ajustement. L'emission diffuse galactique a ete modelisee en 
utilisant le modele de diffus nomme 54_77Xvarli7S base sur GALPROP. L'emission extra- 
galactique a ete, quant a elle, modelisee par une simple loi de puissance. Pour chaque pulsar, 
nous avons selectionne les evenements au-dessus de 0.1 GeV, a I'interieur d'une region de 10° 
centree sur la position de la source, et inclus dans notre modele tous les objets jusqu'a 17°. 
Rappellons que les sources a I'exterieur de la region analysee peuvent contribuer a l'emission 
de basse energie. L'ajustement des parametres spectraux a ete effectue sur la region pulsee 
(ON) du pulsar, en tenant compte de la region non pulsee (OFF), modelisee par une simple loi 
de puissance (possible PWN). On trouvera la definition de la phase du OFF dans le tableau 



9.4 pour chaque objet. Dans la plupart des cas, l'emission du OFF n'est pas significative. Un 



futur papier est dedie a I'etude de l'emission non pulsee (possible PWN) des pulsars detectes. 



Le tableau |9.2| presente les resultats de l'ajustement des parametres spectraux des 46 
pulsars detectes par le LAT. Pour le calcul de la luminosite et de I'efficacite y, nous avons 
suppose un facteur de correction fo = 1 . Notons que les parametres spectraux du pulsar du 



Crabe proviennent de Abdo et al (2009-Crab), qui prend en compte la modelisation de la 



PWN. Signalons que ce travail a ete effectue par I'equipe du catalogue general Fermi. 

Pour un certain nombre de pulsars, le modele spectral incluant une coupure exponentielle 
n'est pas plus significatif qu'une simple loi de puissance. Nous avons determine cela en cal- 
culant le TScutoff = 2Alog (likelihood) (comparable a une distribution de pour un degre de 
liberte) entre les modeles avec et sans coupure. Les pulsars ayant un TS < 10 presente une 
mesure de I'energie de coupure non significative. 



II faut signaler que la connaissance de l'emission diffuse galactique n'est pas parfaite, et 
pent avoir un impact sur le fiux et les parametres spectraux des sources etudiees. Ceci fait 
partie des erreurs systematiques qu'il est difficile d'estimer, etant donne que le modele de 
diffus evolue en fonction de la statistique grandissante des donnees du LAT. II est ainsi inte- 
ressant de comparer les resultats spectraux des trois pulsars individuels etudies aux chapitres 
6, 7, et 8, avec les resultats du catalogue, etant donne que ces etudes utilisent des modeles 
de diffus differents et des intervalles de temps differents. Aux incertitudes statistiques et 
systematiques pres, I'accord entre les analyses individuelles des pulsars et le catalogue est 
excellent. Les conclusions theoriques a partir des observations ne sont done pas erronees entre 
ces deux analyses. Cependant, notons que ces trois pulsars ont un fiux relativement eleve. 
Des tests comparatifs ont ete effectues pour des pulsars presentant un faible fiux, tel que 
PSR J0248+6021 (F ~ 4 x lO^^phcm^^s"^). Dans ce cas, le changement de diffus pent affec- 
ter le fiux du pulsar d'un facteur 1.5. Les observations spectrales apportees par les pulsars 
faibles sont done a prendre avec precaution. 
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9.3 Associations 



Le tableau 9.3 fournit rassociation spatiale des 46 pulsars detectes par le LAT avec des 
sources connues au GeV et au TeV ; 22 des 46 pulsars sont associes a des sources EGRET, 
11 sont des pulsars de type (g), 2 sont des pulsars milliseconde, et 9 sont des pulsars de 
type (r). Les pulsars PSR J 1420— 6048 et J 1418— 6058 sont associes tous les deux a la source 
EGRET 3EG J 1420— 6038 (complexe de Kookaburra). Le LAT a permis de resoudre ces deux 
sources separemment. Un certain nombre de sources EGRET non-identifiees avaient ete au 



prealablement associees a des SNRs, PWNs, ou autres objets (e.g. Walker et al 2003; De 



Becker et al 2005j). II n'est pas etonnant de voir que 19 des 38 pulsars normaux ont des 
associations avec une PWN ou/et une SNR (Roberts et al 2005 Green 2009). Finalement, 



12 des 38 pulsars non-milliseconde sont associes avec des sources au TeV, dont la plupart (9 
sur 12) sont associes avec une PWN. 



9.4 Resultats et Discussion 



9.4.1 Luminosite Gamma 



Le potentiel au niveau des lignes de champs ouvertes, qui est proportionnel au courant de 
Goldreich et Julian, a ete longtemps reconnu comme un parametre important pour discrimer 



les modeles d'emission des pulsars y (Arons 1996). Ce parametre implique que la luminosite 



y est proportionnelle a la racine ca rree de E, c'est a dire, que I'efiicacite y augmente quand le 



ralentissement E diminue. La figure 9.2 presente la luminosite y en fonction du ralentissement 



E basee sur les flux en energie mesures pour chaque pulsar (voir le tableau 9.2). Signalons 
deux points importants. Premierement, la luminosite est tres sensible a la distance estimee 
du pulsar (equation |1.11[ ). Le tableau 2.1 presente la meilleure estimation de la distance pour 



37 des 46 pulsars detectes par le LAT a partir de la mesure de dispersion, de la mesure de 
leur parallaxe, ou de I'observation de I'absorption de I'hydrogene. Les 9 pulsars restants sont 
des pulsars decouverts en y (type g). Le second point est I'hypotliese d'un faisceau uniforme 
sur tout le ciel (fa = 1), ce qui n'est pas realiste pour certaines configurations geometriques 
des modeles d'emission y. 



Sur la figure |9.2[ nous avons indique une ligne pointillee qui suppose une efRcacite de 
conversion de 100% pour transferer la perte d'energie rotationnelle en photons y (Ly = E), 
ainsi qu'une ligne discontinue theorique qui suppose = (10^^ erg s^^ E)^''^. Etant donne les 
larges incertitudes sur la luminosite, il est difficile pour I'instant d'etudier les details de son 
evolution. Nous discuterons cependant de quelques tendances apparentes. 

Les pulsars avec un fort E semblent suivre relativement bien la tendance E^^^, malgre la 
large dispersion de la distance de PSR J2021+3651 et PSR J1709-4429. Dans I'intervalle 
10^^ ergs^^ < E < 10^^-^ ergs^\ la luminosite semble plus ou moins constante, mais la plupart 
des pulsars de cette zone ont une distance basee sur la mesure de dispersion dont I'erreur 
pent depasser 50%. II faut done etre prudent sur cette observation. De plus, I'une des deux 
distances (9kpc) du pulsar PSR J0248+6021 demande une efiicacite de 100% pour produire 
la luminosite observee, ce qui est deraisonnable compare a la tendance generale. Finalement, 
le pulsar PSR J2021+4026 qui apparait avec une grande efiicacite est associe au reste de 
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OPGL J2020. 8+3649 




PWN G75. 2+0.1 1 








MGRO J2019+37 ^ 


J2021 + 4026 


OPGL J2021. 5+4026 


3EG J2020+4017 


Y Cygni 








SNR G078.1+01.8 ^ 


J2032 + 4127 


OPGL J2032. 2+4122 


3EG J2033+4118 


MGRO J2031+41 ^ 






EGR J2033+4117 




J2043 + 2740 








J2124-3358 


OPGL J2124. 7-3358 




PWN GlO.9-45.4 ' 


J2229+6114 


OPGL J2229. 0+6114 


3EG J2227+6122 


PWN G106.6+2.9 ' 






EGR J2227+6114 


MGRO J2228+61 ^'^ ' 


J2238+59 








Tableau 9.3 - 


Associations spatiales des 4.S pulsars detectes avec des sources observees au GeV et au TeV. 


Notes : ^Sources du 


catalogue LAT pour 3 mois de donnees (Abdo et al 2009-hsl) 


^Sources des catalogues EGRET 


{Hartman et al. 1999), EGRET revise (Casandjian and Gremer. 2008), et GeV fLamb and Macomb, 1997). 


Heferences ; i. , Roberts 


et al,^2005j, 2. {Green. ^i:uuy^, J. ^Abao et ai. 


2uuy-mil), 4. ^Aharoman et ai, t'^UUbcj, b. ^Atiaroman et ai, ,2006b^, 6. 


1 Aharoman et al 20dda\, 7 


\hoppe et ai||200y|l, 8. \Aharonian et ai||2007|l, 


a. \U]annati-Atai et al\\2au'l\, 


lU. \Aharoman et al\\iiuuyij, 11. \Uoodman\ 


J and Sinms^^£009^ 
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supernova y Cygni. L'age caracteristique du pulsar est 10 x superieur a I'age de la SNR. La 
distance associee pour le pulsar doit done etre traite avec attention. II est ainsi essentiel de 
maitriser la distance des objets etudies pour suivre revolution de la luminosite y. 
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Figure 9.2 - Luminosite y en fonction du taux de perte d'energie rotationnelle E. La lumi- 
nosite y est calculee d partir du flux en energie integre (G^oo) si en utilisant un facteur de 
correction du faisceau fn = 1 pour tous les pulsars. Nous avons egalement indique par une 
etoile la luminosite totale a haute energie pour le pulsar du Crabe, Ljot = 1-x + Ly . Carres 
bleus : (g) pulsars decouverts en y. Triangles rouges : (m) pulsars y milliseconde. Cercles 
verts : (r) autres pulsars emetteurs radio et y. Les points vides indiquent les pulsars pour 
lesquels la distance est uniquement connue par le DM. Les pulsars ayant deux distances sont 
indiques avec deux points connectes par une ligne pointillee. Ligne pointillee ; Ly = E. Ligne 

discontinue : Ly proportionnelle a VE- 



Pour E compris entre 10^"* et 10^^ erg s^^ , les luminosites sont distribuees sur un large 
intervalle pour les pulsars normaux. Pour les modeles magnet ospheriques externes OG et SG, 
cette zone correspond a un intervalle de saturation. Pour le modele SG, I'interruption dans 
la production de photons y se produit a environ 10^^ erg s^^ , tandis que pour le modele OG 
la saturation est predite pour E ~ lO^'^ergs^^ Avec la statistique et les incertitudes actuelles, 
il n'est pas possible de discriminer ou de rejeter ces deux modeles. Signalons egalement le cas 
particulier du pulsar PSR J0656+1414, qui a une luminosite 30 x inferieure a la valeur predite 
pour une luminosite proportionnelle a E^^^. Une interpretation possible est I'observation de 
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remission y de la calotte polaire du pulsar, qui est prevue pour C < 20° (Everett and Weisberg 



2001 ). II se pent que le profil pulse constitue d'un seul pic et le spectre mou de ce pulsar, soit 



le resultat d'une autre zone d'emission non predite par les modeles. Finalement, en dessous 
de lO^'^ergs^^ la majorite des objets sont des pulsars milliseconde dont la luminosite semble 
etre proportionnelle a E. 



9.4.2 Population 



La figure 



9.3 



presente le diagramme P — P mis a jour avec les 46 pulsars detectes avec le 
LAT, parmi lesquels 21 pulsars normaux radio et y (type r, cercle vert), 17 pulsars decouverts 
en y (type g, carre bleu), et 8 pulsars y milliseconde (type m, triangle rouge). Les points noirs 
indiquent les pulsars pour lesquels une recherche de pulsations a ete effectuee en utilisant 



les ephemerides, tandis que les points gris indiquent les pulsars connus (Manchester et al 



2005;), qui n'ont pas fait I'objet d'une recherche de signaux pulses. Les 46 pulsars detectes 



sont concentres en haut a gauche pour les pulsars normaux, avec une limite de detection 
inferieure a 10^"^ erg s"^ pour le ralentissement t (PSR J0357+32 et PSR J1741-3054), et en 
bas a gauche pour cette nouvelle population d'emetteurs y, les pulsars milliseconde (jAbdo 
et al 2009a). Cependant, avec la statistique actuelle, il est encore difficile de delimiter une 
zone bien precise de production de photons y en fonction du ralentissement E et du champ 
magnetique au cylindre de lumiere. Dans nos dernieres analyses, un petit nombre de pulsars 
avec un E ~ 10^^^^^ erg s^^ ont des significativites de 2 a 3 a (notre limite de detection est 
definie a 5 cr) . 

Rappelons que certains pulsars decouverts en y (type g) ont ete detectes en radio comme 
les pulsars radio faint PSR J1741-2054 (Li.4ghz = 0.03 mJy) et PSR J2032+4127 (jCamilo 



et al, 2009). Des recherches de pulsations radio sur plusieurs autres objets comme Geminga, 



PSR J0007+7303 (CTAl), PSR J1836+5925, n'ont pas ete conclues par une detection (| Kassim 
and Lazio 1999[ Halpern et al 2004). Ceci suggere que remission y a une extension plus large 
que le faisceau radio, comme prevu pour les modeles OG et SG/TPC. En d'autres termes, 
les modeles magnetospheriques externes predisent un rapport des pulsars de type (r) sur les 



pulsars de type (g) beaucoup plus petit que les modeles Polar Cap ( [Harding et al[ 2007). 



9.4.3 Courbes de Lumiere et Energie de Coupure 

La structure des courbes de lumiere y ainsi que le decalage entre les pics a differentes 
longueurs d'onde peuvent nous aider a comprendre la geometrie et la physique de la region de 



remission. En se referant aux profils y des 46 pulsars (Abdo et al, 2009-psrcat) , on remarque 
que la majorite des pulsars montre deux pics relativement etroits, suggerant que Ton observe 
la caustique du pulsar a partir des bords. Quand un unique pic est observe, il a tendance 
a etre large, suggerant une coupure tangentielle du cone d'emission. Cette description est 
realisee par le modele Outer Gap et les modeles Slot Gap et Two Pole Caustic a haute 
altitude. 

A partir des courbes de lumiere des pulsars detectes, deux parametres, dont le decalage en 
phase 5 entre la position du pic radio principal et le premier pic y, ainsi que I'ecart A entre 
les pics y, ont ete extraits (voir la figure [974| et le tableau 9.4). Ces deux parametres semblent 
correler dans les modeles magnetospheriques externes (Romani and Yadigaroglu 1995). Cette 
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Figure 9.3 - Diagramme P — P, distribution des pulsars en fonction de leur periode de 
rotation et de leur derivee. Les lignes vertes representent I' age caracteristique, tandis que 
les lignes rouges representent le taux de perte d'energie de rotation E. Les points noirs 
indiquent les pulsars pour lesquels une recherche de pulsations a ete effectuee en utilisant les 
ephemerides, tandis que les points gris indiquent les pulsars connus ( Manchester et al 2005^. 
Carres bleus : (g) pulsars decouverts en y. Triangles rouges : (m) pulsars y milliseconde. 
Cercles verts : (r) autres pulsars emetteurs radio et y. 
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distribution 5 — A et en particulier la presence de valeurs de A ~ 0.2 — 0.3 favorisent le modele 



OG si Ton compare cela avec les predictions des modeles OG et TPC dans (Watters et al 
2009), tandis que le nombre important de pulsars avec A ~ 0.4 — 0.5 favorisent les modeles 



TPC. 



IS 0.2 



Pulsars avec 1 Pic 
-■ • • - 

1 . 1 1 1 1 1 1 1 1 1 1 1 1 1 . 1 1 1 1 1 1 1 1 1 1 . 1 1 1 1 1 1 1 . 1 1 . 1 1 1 1 . 1 
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Decalage Radio (8) 

(a) Distribution 5 — A 




(b) Distribution Ecutoff - Blc 



Figure 9.4 - Gauche : Separation de la position en phase A des pics y en fonction du 
decalage en phase 5 entre le principal pic radio et le premier pic y. Les pulsars non detectes 
en radio sont distribues avec 8 — 0, tandis que les pulsars avec un simple pic sont distribues 
avec A = 0. Carres bleus : (g) pulsars decouverts en y. Triangles rouges : (m) pulsars y 
milliseconde. Cercles verts : (r) autres pulsars emetteurs radio et y. Droite : Distribution 
de I 'energie de coupure en fonction du champ magnetique au cylindre de lumiere. 



A present, essayons de degager une observation a partir des resultats spectraux. La cou- 
pure spectrale Ecutoff ne montre pas de correlations avec le champ magnetique de surface 



(equation 1.5), suggerant que remission y ne vient pas de la calotte polaire pres du pole 
magnetique (modele Polar Cap). Cependant, cette coupure spectrale semble legerement cor- 



relee au champ magnetique au niveau du cylindre de lumiere (equation 1.6), c'est-a-dire, a 



une altitude haute de I'etoile (voir la figure 9.4). Les valeurs de Ecutoff s'etendent de 1 GeV 
a 10 GeV, et tous les types de pulsars semblent avoir la meme correlation, cela inclut les 
pulsars milliseconde. Les modeles magnetospheriques externes predisent cette observation oil 
remission y provient principalement du rayonnement de courbure des electrons. 



Cette etude nous a permis d'etablir une premiere avance sur la population de pulsars. Nous 
reviendrons sur les conclusions de ce chapitre dans la conclusion generale de cette these. 
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PSR 


Type" 


Multiplicite 


Decalage Radio Separation des Pics y 


Definition OFF 






des pics 


6 


A 


4> 


J0007+7303 


g 


2 




0.23 ± 0.01 


0.29 - 0.87 


J0030+0451 


m 


2 


0.18 ± 0.01 


0.44 ± 0.01 


0.68 - 0.12 


J0205+6449 


r 


2 


0.08 ± 0.01 


0.50 ± 0.01 


0.64 - 0.02 


J02 18+4232 


m 


2 


0.32 ± 0.02 


0.36 ± 0.02 


0.84 - 0.16 


J0248+6021 


r 


1 


0.35 ± 0.01 




0.71 - 0.19 


J0357+32 


g 


1 






0.34 - 0.86 


J0437-4715 


m 


1 


0.43 ± 0.02 




0.60 - 0.20 


J0534+2200 


r 


2 


0.09 ± 0.01 


0.40 ± 0.01 


0.62 - 0.98 


J0613-0200 


m 


1 


0.42 ± 0.01 




0.56 - 0.16 


J0631+1036 


r 


1 


0.54 ± 0.02 




0.80 - 0.20 


J0633+0632 


g 


2 




0.48 ± 0.01 


0.09 - 0.45 


J0633+1746 


g 


2 




0.50 ± 0.01 


0.24 - 0.54 


J0659+1414 


r 


1 


0.21 ± 0.01 




0.40 - 1.00 


J0742-2822 


r 


1 


0.61 ± 0.02 




0.84 - 0.44 


J0751+1807 


m 


1 


0.43 ± 0.02 




0.63 - 0.99 


J0835-4510 


r 


2 


0.13 ± 0.01 


0.43 ± 0.01 


0.66 - 0.06 


J1028-5819 


r 


2 


0.19 ± 0.01 


0.47 ± 0.01 


0.76 - 0.12 


J1048-5832 


r 


2 


0.15 ± 0.01 


0.42 ± 0.02 


0.64 - 0.04 


J1057-5226 


r 


2 


0.35 ± 0.05 


0.20 ± 0.07 


0.72 - 0.20 


J1124-5916 


r 


2 


0.23 ± 0.01 


0.49 ± 0.01 


0.60 - 0.96 


J1418-6058 


g 


2 




0.47 ± 0.01 


0.54 - 0.90 


J1420-6048 


r 


2" 


0.26 ± 0.02 


0.18 ± 0.02 


0.60 - 0.10 


J1459-60 


g 


2 




0.15 ± 0.03 


0.34 - 0.78 


J1509-5850 


r 


2" 


0.18 ± 0.03 


0.20 ± 0.03 


0.52 - 1.00 


J1614-2230 


m 


2 


0.19 ± 0.01 


0.51 ± 0.01 


0.92 - 0.14 


J1709-4429 


r 


2 


0.24 ± 0.01 


0.25 ± 0.01 


0.66 - 0.14 


J1718-3825 


r 


1 


0.42 ± 0.02 




0.68 - 0.20 


J1732-31 


g 


2 




0.42 ± 0.02 


0.49 - 0.93 


J1741-2054 


g 


2 


0.30 ± 0.01 


0.18 ± 0.02 


0.67 - 0.19 


J1744-1134 


m 


1 


0.83 ± 0.02 




0.08 - 0.44 


J1747-2958 


r 


2 


0.18 ± 0.01 


0.42 ± 0.04 


0.64 - 0.10 


J1809-2332 


g 


2 




0.35 ± 0.01 


0.41 - 0.89 


J1813-1246 


g 


2 




0.47 ± 0.02 


0.56 - 0.90 


J1826-1256 


g 


2 




0.47 ± 0.01 


0.54 - 0.94 


J1833-1034 


r 


2 


0.15 ± 0.01 


0.44 ± 0.01 


0.68 - 0.10 


J1836+5925 


g 


2 




0.48 ± 0.01 




J1907+06 


g 


2 




0.40 ± 0.01 


0.46 - 0.94 


J1952+3252 


r 


2 


0.15 ± 0.01 


0.49 ± 0.01 


U.bo — U.Uo 


J1958+2846 


g 


2 




0.45 ± 0.01 


0.55 - 0.95 


J2021+3651 


r 


2 


0.17 ± 0.01 


0.47 ± 0.01 


0.70 - 0.04 


J2021+4026 


g 


2 




0.48 ± 0.01 




J2032+4127 


g 


2 


0.15 ± 0.01 


0.50 ± 0.01 


0.60 - 0.92 


J2043+2740 


r 


2 


0.20 ± 0.01 


0.36 ± 0.01 


0.64 - 0.08 


J2124-3358 


m 


1 


0.86 ± 0.02 




0.92 - 0.58 


J2229+6114 


r 


1 


0.49 ± 0.01 




0.64 - 0.14 


J2238+59 


g 


2 




0.50 ± 0.01 


0.60 - 0.92 


Tableau 9.4 - Parametres (S,A) des 


courbes de lumiere des 46 pulsars 


detectes par le LAT, 


estimes a partir du profil au-dessus de 100 MeV. La demiere colonne indique la definition 


de la region non pulsee 


phasogramme. 






°- (r) pulsars decouverts 


en radio, (g) pulsars decouverts en gamma, (m) pulsars milliseconde 


^ Pour certains profils, 


les donnees actuelles ne permettent pas de discriminer entre un large 


pic et deux pics 


non-resolus (voir une discussion dans 


Weltevrede et al 


(200^). 
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Le Large Area Telescope a bord du satellite Fermi, en orbite depuis le 11 juin 2008, 
est un instrument destine a I'exploration du del des hautes energies (> 20MeV). Grace a sa 
grande surface efiicace, son pouvoir de localisation, et sa precision temporelle, le telescope 
Fermi ouvre une nouvelle fenetre pour I'etude passionnante de I'acceleration de particules 
chargees a travers les cavites acceleratrices de la magnetosphere des pulsars. Ces objets sont 
une etape dans revolution stellaire et sont d'excellents candidats pour localiser les nebuleuses 
a vent de pulsar et les restes de supernova, sources y possibles pour I'acceleration des rayons 
cosmiques dans notre univers. 

Cette these constitue les premieres observations des pulsars y avec le Large Area Te- 
lescope. Nous avons mis en place une liste de candidats connus en radio ou en X, basee sur 
I'energie de ralentissement et la distance de ces objets. Une analyse sur la recherche de signal 
pulse a partir des donnees du LAT a ete developpee pour analyser les 767 ephemerides de 
pulsars fournies par les radiotelescopes et les telescopes X. Au total 22 pulsars radio relative- 
ment jeunes ont ete detectes dans le domaine des hautes energies pour les six premiers mois 
de la mission Fermi. On ajoutera a cela la detection de nouvelles populations de pulsars, dont 
8 pulsars milliseconde, et 16 peu ou pas emetteurs radio. Le nombre de pulsars y connus a 
done evolue d'un ordre de grandeur depuis la fin de la mission EGRET, fournissant une nou- 
velle voie pour la demographie, la physique des pulsars, et la formation des etoiles a neutrons. 

La forme des courbes de lumieres y des 46 pulsars detectes montre une interessante uni- 
formite. 75% des profils y pulses presentent deux pics separes par plus de 0.2 en phase, 
tandis que les 25% restants semblent montrer un seul et unique pic. Les spectres en energie 
de remission pulsee peuvent etre decrits par des lois de puissance avec une coupure exponen- 
tielle. Cette coupure en energie s'etend de ~ 1 a 6 GeV, limitant actuellement I'observation 
de ces objets par les telescopes y au sol. De plus, leur luminosite a haute energie semble 
occuper une grande fraction de I'energie de rotation, et augmente a mesure que celle-ci di- 
minue jusqu'a une limite de saturation E ~ 1 0^^"^'* erg s~^ Ces observations sur les profils et 
les spectres en energie pulses suggerent que la zone d'emission y se trouve a haute altitude 
par rapport a I'etoile a neutrons et pres du cylindre de lumiere, meme si pour certains cas 
particuliers, remission pourrait venir de la zone au-dessus des poles de I'etoile. 

II s'avere qu'une large fraction de ces decouvertes est correlee a des sources telles que 
des PWN et/ou des SNRs observees en X et/ou au TeV. II sera ainsi essentiel de suivre les 
systemes avec les instruments du TeV et des X pour contraindre les differents modeles d'emis- 
sion. De plus, la maitrise du ciel y ouvre une nouvelle perspective pour explorer les chocs 
entre les vents stellaires des associations OB et des etoiles massives Wolf-Rayet, ou encore les 
jets des systemes binaires X. Ces objets sont des candidats potentiels pour abriter un pulsar, 
lis contribuent egalement au fond galactique pour la recherche de signaux plus subtils prove- 
nant d'autres types de sources telle que la matiere noire. Enfin, cette nouvelle population de 
pulsars devoile peu a peu I'origine d'une partie des sources EGRET non-identifiees a basse 



139 



CONCLUSION 



latitude galactique. II semble en effet que les jeunes pulsars formes au sein des supernovae ou 
par accretion de matiere dominent la population des sources y de notre Galaxie. 

Les resultats de cette these soutiennent I'importance de continuer I'etude de cette popu- 
lation de pulsars grandissante. Le rapport du nombre de pulsars normaux vus en radio et/ou 
en y permet de mettre des contraintes sur la geometrie des faisceaux, sur les sites d'acce- 
leration des particules, et eventuellement sur revolution de leurs proprietes d'emission au 
cours du temps. II sera egalement essentiel de rechercher des pulsars millisecondes dans les 
donnees LAT. De plus, les spectres finement resolus en phase vont a present nous permettre 
de definir plus en details les modeles d'emission, tout en etudiant les profils pulses des autres 
longueurs d'onde. Finalement, le debat sur I'association entre le pulsar PSR J0205+6449 et 
la supernova historique nous rappelle qu'il est essentiel d'etudier les associations, d'une part 
pour revolution des systemes stellaires, et d'autre part pour evaluer les distances de notre 
Galaxie. 
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Annexe A 



Premiere annexe 



Cette annexe presente un echantillon des courbes de lumiere du premier catalogue de 



pulsars du telescope Fermi (Abdo et al 2009-psrcat) 



O PSR J0007+7303 decouvert en y dans le reste de supernova CTAl (Abdo et al 2008) 



le pulsar millisecond PSR J0030+0451 (Abdo et al 2009b) 



© le pulsar radio PSR J0659+1414 associe au reste de supernova Monogem ring 



Counts 



Counts 




Figure A.l - Courbes de lumiere de PSR J0007+ 7303 (? = 316ms;. 
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Resume : Le Large Area Telescope a bord du satellite Fermi, lance le 11 juin 2008, est 
un telescope spatial observant I'univers des hautes energies. L'instrument couvre I'intervalle 
en energie de 20 MeV a 300 GeV avec une sensibilite nettement amelioree et la capacite de 
localiser des sources ponctuelles. II detecte les photons y par leur conversion en paire electron- 
positron, et mesure leur direction et leur energie grace a un trajectographe et un calorimetre. 

Cette these presente les courbes de lumieres et les mesures spectrales resolues en phase des 
pulsars radio et gamma detectes par le LAT. La mesure des parametres spectraux (flux, indice 
spectral, et energie de coupure) depend des fonctions de reponse de l'instrument (IRFs). Une 
methode developpee pour la validation en orbite de la surface efiicace est presentee en utilisant 
le pulsar de Vela. Les efScacites des coupures entre les donnees du LAT et les donnees simulees 
sont comparees a chaque niveau de la rejection du fond. Les resultats de cette analyse sont 
propages vers les IRFs pour evaluer les systematiques des mesures spectrales. 

La derniere partie de cette these presente les decouvertes de nouveaux pulsars y individuels 
tels que PSR J0205+6449, J2229+6114, et J1048-5832 a partir des donnees du LAT et des 
ephemerides radio et X. Des analyses temporelles et spectrales sont investies dans le but de 
contraindre les modeles d'emission gamma. Finalement, nous discutons les proprietes d'une 
large population de pulsars gamma detectes par le LAT, incluant les pulsars normaux et les 
pulsars milliseconde. 

Mots-cles : Pulsars - Fermi Large Area Telescope - Surface Efficace 

Abstract : The Large Area Telescope (LAT) on Fermi, launched on 2008 June 11, is a space 
telescope to explore the high energy y-ray universe. The instrument covers the energy range 
from 20 MeV to 300 GeV with greatly improved sensitivity and ability to localize y-ray point 
sources. It detects y-rays through conversion to electron-positron pairs and measurement of 
their direction in a tracker and their energy in a calorimeter. 

This thesis presents the y-ray light curves and the phase-resolved spectral measurements 
of radio-loud gamma-ray pulsars detected by the LAT. The measurement of pulsar spectral 
parameters (i.e. integrated flux, spectral index, and energy cut-off) depends on the instrument 
response functions (IRFs). A method developed for the on-orbit validation of the effective 
area is presented using the Vela pulsar. The cut efficiencies between the real data and the 
simulated data are compared at each stage of the background rejection. The results are then 
propagated to the IRFs, allowing the systematic uncertainties of the spectral parameters to 
be estimated. 

The last part of this thesis presents the discoveries, using both the LAT observations 
and the radio and X ephemeredes, of new individual y-ray pulsars such as PSR J0205-I-6449, 
and the Vela-like pulsars J2229-I-6114 and J1048— 5832. Timing and spectral analysis are 
investigated in order to constrain the y-ray emission model. In addition, we discuss the 
properties of a large population of y-ray pulsars detected by the LAT, including normal 
pulsars, and millisecond pulsars. 

Keywords : Pulsars - Fermi Large Area Telescope - Effective Area 



